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Über die Fortpflanzung kleiner Störungen in flüssigen Medien, 
betrachtet als Ausbreitungsvorgang von Unstetigkeiten in den 
Lösungen der Bewegungsdifferentialgleichung”) 


Von Johannes Nitsche in Berlin 


Die Fortpflanzung einer kleinen punkthaften Störung in zwei längs ebener Trrennfläche aneinandergrenzen- 
den Flüssigkeiten wird nach der Methode der Unstetigkeiten als Ausbreitungsvorgang kleiner Druck- und Ge-. 
schwindigkeitssprünge behandelt. Die Verteilung des Drucksprunges auf der Stoßfront wird angegeben und 
genauer diskutiert für den Fall, daß die Störung von einem Punkt der Trennfläche ausgeht. Die Verhältnisse 
bei auftretender Totalreflexion werden an einem Beispiel betrachtet. 


The propagation of a small point perturbation in two liquids, adjacent along a plane, is treated as pro- 
pagation process of small changes of pressure and velocity, according to the method of discontinuities. The 
distribution ofthe jump of pressure at the shock frontis given, and the case is discussed in detail, that the per- 
turbation emerges from the separating plane. The conditions of the occurence of total reflection are considered 
in an example. 


La propagation d’une petite perturbation en forme de point, proc&dant dans deux liquides adjacents le long 
d’un plan separant, est traitee comme procede de propagation de petits changes de pression et de velocite, suivant 
la methode des discontinuites. La distribution du choc de pression sur le front de pression est indiquee et dis- 
cutee en detail pour le cas que les perturbations sortent d’un point du plan separant. Les conditions se pro- 
duisant & l’occurence d’une reflexion totale sont considerees & parlir a un cexemple. 


PacnpocrpaHeHue MaI0OTO TOYeYHOTO BO3MYINCHHA B ABYX CMEKHEIX 3KUNKOCTAX, Pasrpa- 
HHYEHHBIX ILIOCKOH IOBEPXHOCTBEO Pasıena, TPAKTYeTCA MO MeTony PAsPEIBOB, KaK IIPOMecE 
PacupocTpaHeHHnA HeÖONbIIHX CKAYKOB NABTeHHA U CKOPocTu. JlaeTca pacııpeneneHue CKA4K& 
AaBJICHHA IIO HAIOPHOMY PPOHTy H AeTalbHO OÖCY:ENdeTcH B TOM CIY4&e, KOTNA BO3MYINeHHE 
HCXOAUT M3 ONHOK TOYEH NOBePXHOCTH pasıeıa, YCAIOBHA IPH BOSHHKHOBEHHH IOJIHOTO OTPA- 
’KeHUA PACCMOTpeHEI HA IIPHMepe. 


In der vorliegenden Arbeit soll die Ausbreitung einer kleinen punktförmigen Störung in 
zwei längs einer ebenen Trennfläche aneinandergrenzenden Flüssigkeiten betrachtet werden. 
Insbesondere interessiert der Fall, daß die Störung von einem Punkt in unmittelbarer Nähe dieser 
Trennfläche ausgeht. Dabei tritt das Phänomen der v. Schmidtschen Kopfwelle auf (O. v. 
Schmidt [11]!)). Schon das Huyghens sche Prinzip läßt den Verlauf der Wellenfronten 
verstehen, wie sie auf den v. Schmidtschen Photographien erkennbar sind. Für die skiz- 
zierten Verhältnisse soll die akustische Näherung ausreichend sein. Das bedeutet, daß allein den 
Störkomponenten quadratischen Glieder gegenüber den linearen vernachlässigt werden. Das 
Problem ist (für periodische oder punktförmige Störungen) in der Literatur schon mehrfach 
behandelt worden?). In jenen Arbeiten wird für das gesamte Störungsfeld eine Integraldarstellung 
hergeleitet und diskutiert. Dagegen wird hier das Problem als Ausbreitungsvorgang von Stößen 
behandelt, d.h. von Unstetigkeiten in den ersten Ableitungen des Geschwindigkeitspotentials. 
Diese Methode ist schon mehrfach auf andere, einfachere Fragestellungen angewandt 3) worden. 
Unter allgemeinen Voraussetzungen (physikalischen Erhaltungssätzen) können solche Unstetig- 
keiten nur beim Überschreiten einer charakteristischen Fläche auftreten, und die Sprungintensität 
(z. B. der Drucksprung) genügt dann längs der Bicharakteristiken einer gewöhnlichen linearen 
Differentialgleichung. Die Methode liefert unmittelbar die physikalisch hervorstechenden Eigen- 
schaften; allerdings ergibt sie nicht das gesamte Störungsfeld. Das Ziel der Untersuchung besteht 
vielmehr darin, unter Verzicht auf die Kenntnis der Gesamtlösung, die Intensitätsverteilung der 
Störung längs der Stoßfronten zu einem beliebigen Zeitpunkt anzugeben. Zu diesem Zwecke 
muß die erwähnte Gleichung für die Änderung der Sprungintensität jeweils bis zur Stoßfront im 
vorgelegten Zeitpunkt integriert werden. An der Grenzfläche der beiden Flüssigkeiten sind dabei 


*) Herrn Professor Tollmien binich für die Anregung zu dieser Arbeit und die mir gebotene Möglich- 
keit zu ihrer Durchführung im Max-Planck-Institut für Strömungsforschung in Göttingen sehr dankbar. 

1) Dis eckigen Klammern beziehen sich auf das Literaturverzeichnis am Schluß der Arbeit. 

2) Vgl.H. Ott [8], F Sauter [9], A. Schoch [12]. 

3) Vgl. E.B. Christoffel fl, 2], er behandelt auch die Reflexion an einer festen Wand. Weiterhin 
seienerwähnt J. Hadamard [6],G. Zemple&n [14],A. E.H. Love [7], R. Courant und D.Hil- 
bert [3]. 

25 


F ne > Z. angew. Math. Mech, 
362 Nitsche, Über die Fortpflanzung kleiner Störungen in flüssigen Medien 34.34 Nr. 10/11 Okt./Nov. 1954 


die Übergangsrelationen für den einfallenden, reflektierten und durchgehenden Stoß in geeigneter 
Weise abzuleiten und zu verwenden. Das Brechungsgesetz, das ja von geometrischer Natur ist, 
folgt schon aus dem auch hier anwendbaren Huyghensschen Prinzip. 

Es zeigt sich, daß für den Grenzfall, bei welchem die punktförmige Störung von der Trenn- 
fläche selbst ausgeht, diev. Schmidt sche Kopfwelle überhaupt nicht auftritt. Die Verteilung 
des Drucksprunges auf der Stoßfrontim schalldichteren Medium ist dann unstetig, man vergleiche 
Bild 9 am Schluß der Arbeit. Dieses Resultat steht nicht im Widerspruch zu den Experimenten; 
denn einerseits besitzt das Störungszentrum in Wirklichkeit stets eine endliche, wenn auch evtl. 
sehr kleine Ausdehnung, und andererseits sind die in den Experimenten ausgelösten Störungen 
möglicherweise so stark, daß die Rechnung mit den linearisierten Gleichungen nicht statthaft ist. 
(V. Schmidt arbeitet mit Knallfunken, vgl. auch die Bemerkung bei W. Schaaffs [10], 
S. 134). 

= Fall, daß die Störung von einem Punkt der schalldichteren Flüssigkeit ausgeht, ist 
nicht behandelt. Es scheint, als ob die dann für gewisse Winkelbereiche auftretende Total- 
reflexion, bei derjaschon das Huyghens sche Prinzip versagt, von der Methode der Unstetig- 
keiten nicht erfaßbar ist: scharfe Sprünge der einfallenden Welle können bei Reflexion und Bre- 
chung stark verwischt werden, andererseits treten Unendlichkeitsstellen der Lösungen. auf. Diese 
Verhältnisse werden am Schluß der Arbeit an einem einfachen Beispiel skizziert, man vergleiche 
auch eine in diese Richtung gehende Arbeit von F. A. Fischer [4]. 


Der ebene Fall 
I 


Das Problem ist an sich axialsymmetrisch. Einfachheitshalber soll aber der ebene Fall 
ausführlich diskutiert werden; der axialsymmetrische Fall ist dann ganz analog zu behandeln. 
Es sei erlaubt, die Methode der Unstetigkeiten für unsere Zwecke kurz darzustellen. Für die 
einführenden Bemerkungen wird zunächst eine homogene Flüssigkeit ohne Einfluß von Reibung 
und Schwere ins Auge gefaßt. p und o bedeuten Druck und Dichte, ? und 9 deren Mittelwerte. 
Die Differenzen p — ?, o— 0 sollen klein sein. c ist die Ausbreitungsgeschwindigkeit kleiner 
(stetiger) Störungen. Die Verhältnisse in der Flüssigkeit werden von dem Geschwindigkeits- 
potential 9 = 9(z, yY, t) beherrscht, aus dem sich auch Druck und Dichte gemäß 


_5_;%9 RL 
berechnen. (x, y,t) genügt der Wellengleichung 
2 op 100 
Ge Zu Mar aus ng vibaaandl2E 


Nun soll die Ausbreitung eines Drucksprunges, d. h. zugleich die Ausbreitung von Unstetig- 
keiten in den Geschwindigkeitskomponenten u — 9,, 9—= 9, der Flüssigkeitsteilchen, verfolgt 
werden. Die mit der Zeit t fortschreitende (als zweimal stetig differenzierbar vorausgesetzte) 
Stoßfront läßt sich in der Form f(z, y,t) = 0 darstellen. Zu einem Zeitpunkt i—= i, herrsche 
außerhalb der Stoßfront f(x, Y, t,) = 0 insbesondere Ruhe; innerhalb der Stoßfront in ihrer 
Nachbarschaft sollen sich aber endliche (kleine) Geschwindigkeiten vorfinden. Die (Normal-) 
Geschwindigkeit der Stoßfront werde vorläufig mit a bezeichnet. Diese Annahme zieht die 
partielle Differentialgleichung 


1 
k+h 37 9 neuen ee ee 
nach sich. Bei der Festsetzung f, — — a Yf} + fy gibt der Vektor {fa» fy} die Fortschreitungs- 
. . AT, 
richtung der Stoßfront an. Die zugehörige „äußere“ Seite werde durch den Index »— , die 


andere „innere“ Seite durch den Index „+“ gekennzeichnet. Für eine Größe n(&, Y, t), die sich 
beim Überschreiten der Stoßfront sprunghaft ändert, sollen diese Indizes die Grenzwerte bei 


Annäherung von der zugehörigen Seite her bedeuten. Wie üblich wird noch eine ‚eckige Klammer 
als Symbol der „Sprungbildung“ eingeführt: j 


Meter 

Die Gleichung /(z, y,t)—= 0 kann auch als Gleichung einer Fläche (wir sagen hier: der Wellen- 
fläche) im dreidimensionalen z, %, i-Raum mit dem Normalenvektor I fı} aufgefaßt werden 
Es handelt sich zunächst darum, die Relationen aufzustellen, denen die Sprünge in Druck 

und Geschwindigkeit der Flüssigkeitsteilchen beim Überschreiten der Stoßfront genügen müssen 
Ihre Herleitung geschieht ganz analog zu der beim Fall eines Verdichtungsstoßes in einer kom- 


ad amd Bd a nd Bl m ln nn a u 5 


Kau- Zaza> ua 


a um 
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pressiblen Flüssigkeit, nur werden jetzt allein die in den Störgrößen linearen Glieder berück- 
sichtigt. Man fordert, daß die physikalischen Erhaltungssätze (hier Kontinuitätsgleichung und 
Impulssatz) auch trotz der linienhaften Unstetigkeiten ihre Gültigkeit behalten ®). Die Ge- 
schwindigkeit der Flüssigkeitsteilchen in Richtung der Stoßfrontnormalen heißt 


1 ‚ 0p 0% 
A I A 
Emo ta 
Die Kontinuitätsgleichung führt (nach Multiplikation mit Yf3 + f}) auf die Bedingung 
pp ri 0 0 al 
I (1.3 +hay—alit R) -/e pe +hay—al+ 3) DZ ©. 
Aus dem Impulssatz folgt (nach Multiplikation mit (f2 + f3)): 
de Lou ö Pasta? | 
lehrer n Erna] +4 mr 
0 0 0) 0] ac N ME ‚ 
eltern lerne VER) +tme 0, 
Bezüglich der Größenordnung der einzelnen Glieder ist folgendes zu sagen: Wir betrachten eine 


wirklich fortschreitende Stoßfront. Erfahrungsgemäß besitzt ihre Geschwindigkeit a die gleiche 
Größenordnung wie c. Aus der Betrachtung der Dichten für kinetische und innere Energie: 


T-5 +9), Ve A (6) 


folgt, daß o,, 9,» = > = 9; die gleiche Größenordnung haben. Die Bedingungen (5) und (5’) 


führen daher, wenn nur die linearen Glieder berücksichtigt werden 5), auf 


keit) aikthlel:--: Her. (6), 
Beth +RlEl ::.:.:... er (0 

Wegen (1) gehen diese Gleichungen über in 
Rn Te 
Adern et: 


[p;] muß von Null verschieden sein, weil sonst auch [Y,] und [9,] verschwänden, ein Fall, der aus- 
geschlossen wird. Somit ergibt sich 


ae er Se - Brandner (8), 


d.h. in Worten: die Störung breitet sich mit der konstanten Geschwindigkeit c aus. Die Wellen- 
fläche f(x, %, t) = 0 muß also eine Charakteristik der Wellengleichung sein. Die Stoßfront durch- 
läuft im Verlaufe der Zeit alle Kurven einer Parallelkurvenschar; ihre gemeinsamen Normalen 
entstehen durch Projektion der (geradlinigen) Bicharakteristiken der Wellenfläche f(x, y, t) = 0 
des x, y, t-Raumes auf die x, y-Ebene. 

Der auf die zur Stoßfront tangentiale Richtung angewendete Impulssatz liefert als Resultat 
die Stetigkeit der Tangentialgeschwindigkeit, d.h. 


DerBlslekelsean , 2... Ss. ON 


In Verbindung mit (3), (7), (8) lassen sich die Sprungrelationen somit in der folgenden symme- 
trischen Schreibweise zusammenfassen: 


BEE rel. ed. ee 0), 


e,)=0'fs, 9,1 8:15 [91 9 ft a (11). 
Der Proportionalitätsfaktor © ist dabei nun als Maß der Unstetigkeit aufzufassen. Er hängt von 
der Normierung des Vektors {f,, fy, fı} ab. 


4) Vgl.z.B.E.B. Christoffel [1], S.83—89; F.Schubert[l3]. 
5) Man vergleiche die Diskussion bei E. B. Christoffel [1], S. 86—88. 
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Es soll nun noch die von der Längeneinheit der Stoßfront übertragene Leistung L berechnet 
werden. Man betrachtet dazu ein kleines „flüssiges“ Trapez ds - c dt, über welches die Stoßfront 
hinwegstreicht (vgl. Bild 1). Nach (6) erhält man dann 


I 
Miytra0 1 [2 Bla + ld + za [a] ds ode 
oder mit den Formeln (11): | w 
RE u ar (12). 
EN 6 
ds [6 
en) Wir werden später die Normierung % +, =1, fi =—e annehmen. 
“4 F An sich ist es inkonsequent, nach den vorausgegangenen Vernach- 
Bild1 lässigungen einen in den Unstetigkeiten quadratischen Ausdruck zu 


bilden. Das geschieht nur der Einfachheit halber: die Lösung der Trans- 
portgleichung für die Sprunggröße legt die Benutzung von ©? nahe. Für den Ausdruck 


(ds, — Bogenelement der Stoßfront) gilt aus physikalischen Gründen ein Erhaltungssatz, wie 

sich später auch auf mathematischem Wege ergeben wird (vgl. Bild 2, n hängt mit?durchn — et 
zusammen, siehe weiter unten). 

Man kann die im Vorhergehenden abgeleiteten Bedingungen (11) für das Geschwindigkeits- 

potential @ formal folgendermaßen aussprechen ®): 9—=g(x, y,t) wird als stetige Lösung der 

Wellengleichung (2) gesucht, welche beiderseits der Fläche f(z, y, t) = 0 

— jeweils mit Einschluß dieser Fläche — zweimal stetig differenzierbar 


ist, deren erste Ableitungen sich jedoch bei ihrem Überschreiten sprung- 
\\ A haft ändern. Dann bleiben die beiden ‚inneren‘ Ableitungen 
gr 
N m URAN ie 
\ Vi er lang V> hy ha 


stetig. Jede von diesen unabhängige ‚„herausführende‘‘ Ableitung 
macht jedoch einen Sprung. Die Charakteristikenbedingung (10) bzw. 
Bild2 (8) kann auch mit Hilfe des Hamiltonschen Prinzips gewonnen 
werden ?). 


u 
Gl. (10) ist eine partielle Differentialgleichung für die charakteristische Fläche; ihre Bi- 
charakteristiken sind die Strahlen mit den Richtungsvektoren 12 I» — & 1) . Als Ausbreitungs- 
C 


gesetz für © längs dieser Strahlen ergibt sich bekanntlich eine gewöhnliche Differentialgleichung, 
welche auch in der Form 


() 
NUN EEE li 


geschrieben werden kann.S). Dabei ist die charakteristische Fläche f(x, y,t) = 0 in eine Flächen- 
schar f(z, y, t) = const eingebettet gedacht. 


Die durch eine vorgegebene Anfangskurve x — &(£), Ne Y(E), t — 0 verlaufende Lösungs- 
fläche von (10) heißt in Parameterdarstellung 


Er 
y-yEn)=IO+n:sina(d) 


} | ME SR 
= n° ri 
% = a(£) muß der Streifenbedingung 
rd 
PRauNTTE BRIAN, ER a SER SER BER EN FE 2 PT ne 3! 


°) Vgl. Courant-Hilbert [3], S. 356 ff. 
?) Vgl. G. Zemplen [15]. 
®) Vgl. Courant-Hilbert [3], S. 362—364. 
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genügen (d. h. Strahlen senkrecht zur Stoßfront). Es werden nur reguläre Teile der Wellenfläche 


betrachtet, d. h. solche, in denen die Funktionaldeterminante 


%(z,%Y) ) : [e) 
2 NEN) — % SINN — YEeCOSa—NOEı. . 2... 17 
N (En) & € 1%e ( ) 
von Null verschieden ist. Die in der &, y-Ebene fortschreitende Stoßfront erhält man, indem 
man in (15) n = ct einsetzt. Das Differential ihrer Bogenlänge heißt 


de, = Va+typE=AEndE..... can. ... (18). 


Die Einbättung der durch (15) gegebenen charakteristischen Fläche I= 0 in eine Flächenschar 
f= eonst kann man z. B. so vornehmen, daß man unter Einführung eines weiteren Parameters 


1 
(in Yılurch =, zn + & ersetzt. Bezeichnet man die Auflösung nach & und n der ersten beiden 


Formeln von (15) dufeh &=£&(%, y),n=n(, y), so kann man die Flächenschar f(x , y , = const 
euch in der Form 


1 
fa, y,t) = t— = YES NEE Re . (19) 

darstellen. Für die Normalenvektoren ergibt sich 
1 Tr 
Ne: a ve, - 2 ; | a a EEE (20). 


Setzt man diese Werte für f,, fy, fr in die Transportgleichung (14) ein, so findet man die Differen- 
tialgleichung 


OT 
Ye, + 7 ee (21) 
Ihre Lösung heißt 
A(&, no) 
OXE,n)= IHEnM) —-——°) : 20. . sonerz (22) 
Wegen (18) erweist sich nun auch mathematisch die Gültigkeit des schon eingangs erwähnten 
Erhaltungssatzes 
uw Kin)es= Men) ER Mar ehr ki Ma (23), 
bzw. in integraler Form 
& 
S L(&, n) ds, = const (unabhängigvonn) . ......- (24), 
& 


vgl. Bild 2. 
Il 

Nunmehr soll die Ausbreitung von Stößen in zwei verschiedenen Flüssigkeiten (1) und (2) 
untersucht werden, welche längs einer Diskontinuitätslinie ® aneinandergrenzen. % falle im 
Ruhezustand mit der y-Achse zusammen. Es gelte ec,> c,. Die Störung gehe zur Zeit = 0 
vom Punkt = —k<0,y=0 aus. Im Folgenden tritt der Ansatz einer „punkthaften‘“ 
Unstetigkeit nicht explizit auf. Man braucht sich nur ein Anfangswertproblem vorzustellen, bei 
dem die Ausgangslage der Stoßfront ein sehr kleiner Kreis ist. Sein Radius geht, wie sich zeigen 
wird, in die Endformeln nicht ein. Die Unstetigkeiten in den Anfangsdaten können bereits so 
angenommen werden, daß nur eine Sprungausbreitung nach außen stattfindet, d.h. daß von den 
an sich möglichen zwei Stoßfronten die anfänglich kontrahierende nicht auftritt. Man überzeugt 
sich, daß dieser Fall beispielsweise bei den folgenden Anfangsdaten vorliegt (e< 1): 


— [@+k®2+y] für +2 + yP<e 


0 für ++? >, 
1; 
nz für % En k 2 . 2 RS e? 
9, Y 0) — 0] ( y 
| 0 für «+ KH yP? > 
Der Überdruck im Störungszentrum beträgt hier P. 


In der Flüssigkeit (1) existiert ein Geschwindigkeitspotential 9,(@, %, t), in der Flüssigkeit 
(2) entsprechend ein Geschwindigkeitspotential 9,(&, y, t). Aus der Forderung, daß der Druck 


9) Eine analoge Formelfindetsich z. B.beiE.B. Christoffel [11],8.100; F.G. Friedlander[5] 


2 
o(z, y, 0) = 280, C, 
’ 185). 
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und die Normalkomponente der Geschwindigkeit beim Durchgang durch # stetig bleiben müssen, 
folgen für die Grenzwerte bei Annäherung an # die beiden Bedingungen 


ee = BB. 
=) BE a oa en 


ee), er ) Tall. als > == m (Zi em). 


Um diese Übergangsbedingungen (26) und (27) an der Diskon- 
tinuitätslinie 9 zu befriedigen, muß man neben der einfallenden (e) 
in jeder der beiden Flüssigkeiten noch eine weitere Stoßfront an- 
setzen: die reflektierte (r) in der Flüssigkeit (1) und die gebro- 
chene (g) in der Flüssigkeit (2). Es handelt sich um die Bestim- 
mung der zugehörigen Sprunggrößen ©, und ©, aus ©,. Für die 
Normalenvektoren {fz, fy, fı} und {f2, f,, f7} hat man aus geometri- 
Bild 8 schen Gründen das Reflexions- und Brechungsgesetz 


a a ” : (0890, > 0) Pa . (28). 


SIN% .Cı 0 


K=Z=NT— 0; 


Mit der schon in Abschnitt I erwähnten Normierung, läßt sich also schreiben: 


= COS o,, sine, E=—6, 
Bere snag) era HIT EIER (29). 
2 =.6003 0% vesın az De or 


Bild 4 


Bei Beachtung der Vorzeichenfestsetzung ergibt sich an Hand des Bildes 4 und der Übergangs- 
relationen (26) und (27): 


>), ,=-Ro+re+:-. a a er 


er ee) = Re re 


en (93) -( u 99,\”) e 
2), 1), a 


nr Da & 7 
ah rn 


2 2 


und 


Die durch Punkte angedeuteten Glieder verschwinden aus Stetigkeitsgründen, wenn die Punkte 
A, A’, B, B’ nach dem Punkt P zusammenrücken. Dann findet man die beiden Bedingungen 


BO=BO HN BEN -ERAHEHRN 0.2... 
oder, unter Berücksichtigung von (29): 

9, 608 &, = (9, — 9,) cos o,, 9, EHI) aa. (31). 
Die Auflösung heißt 


0, = R(&,) : 9; ; 9, Aa (32) 


ie dr Wellenflächen erhält man die Parameterdustelungen 
a=—k Frau ı nr 


für ncos,<k 


we 


für n cos w>k 


r 2 j 4 c AR 
(9); ya e(1 2) 180, +n2sino, für 2,508 du = EAN (34,). 
T — N . Vjeg 
Dabei gilt 

080, = -yı-2 —, sin? 0 (35). 

Weiter findet man, wenn man den Parameter & in ds) durch FR ersetzt: 
49=—n, A nn SAD -—— —— 1 =: +no,cos x . (36). 

1 


2. cos? &, 
1—5 sin? , 
e e@ 
7 


Wir interessieren uns für die Verhältnisse zu einer Zeit =r. Die zugehörigen Stoßfronten 
ergeben sich, indem bei (e) und (r) für n = c, r bei (g) für n = cz; r eingesetzt wird. 

Auf den so mit Hilfe des Parameters & = «, dargestellten Stoßfronten ist nun die Ver- 
teilung von L zu berechnen. Dabei wird jetzt angenommen, daß ©, für sehr kleine (positive) 
Zeiten längs der ganzen Ausgangsstoßfront gleichmäßig verteilt ist. Es sei etwa zur Zeit 


= 0<kle 
: Del) TAN wen. Aukraknan Aut De (37) 


unabhängig von a. Mit den in (25) auftretenden Konstanten e und P hängen ö und ©, folgender- 


maßen zusammen: 
€ P 
Ö =, Pe 

cı Bi cı 
(Der krärücki im Störungszentrum beträgt bei dem Beispiel (25) also P=9, c, 9,). Nach (22) 
findet man 


PPFRAPP \(e): 
Ok) = 2-0, 
also 
re ee: 
2PPU TFT): 
Hung "9 Cosa Mesa 9 0 Cosa) mcos« 
A a s„cıöcosa k _ a6 pe 
— R?(o) Gere ann MT, R?(a) : 9 ; 
also 


Man = ERW... . Bor dan Dil „1 35440): 


\ 10) Dassind die bekannten Übergangsbedingungen; vgl. z. B. (für periodische Störungen) A.Schoch[12], 
Ss. 141. 


angew. Math. Mech. 
368 Nitsche, Über die Fortpflanzung kleiner Störungen in flüssigen Medien pa4.34 Nr. zo Okt./Nov. 1954 
20,0 > > 2 
Ele. NEST OT ES En 


3.2Knalg): 1 
[23 k ) anf Ca ) 
) u ee O5 
5 h 6, k 4% (« 01 0050) _ G24).08 (x a \7 61 0080) 
o(& n) == 9% 3 6, Cc0S / AYKo, n) I: e\ cos“ AP(c, n) 
’ 40 (« = ) 
c,6c0Sa C; & 
=) HT Ta 
Nun gilt Fer 
aus, % k )-— Mae on Wale sin? , 
6, Cosa cos? & ed 
so daß sich 


2 
1— a sin? & 
9,(%; 63 7) = I @ (&) ° EPOPYRFGIGE BEIN PET. el er Bi N 
2 ee Ze gi 
ae il! 2) 
ergibt. Führt man noch die ee ein: 
e=6 Nadel een he RR: 
1 

so erhält man schließlich 


R k 
I, = Se für cosa<s mn 
k 
ep | für cosa> — 
5 a : 
” 1 > sin? « 
Die,n)= A Ga) ———ı  —— für coso> ss 
De 2%C, 70,6 Eos + k ee DH 
q 6, TCOS“& d 


wobei R>(o) und @°(&) aus (33) zu entnehmen sind. 
Mit den Differentialen ds®, ds” und ds‘ der Bogenlänge auf den drei Stoßfronten findet 
man durch Einsetzen nach einigen Umformungen nun auch den Erhaltungssatz 


[LICH [LICHT URS, (44). 
(e) (r) (9) 


IV 
Wenn % sehr klein ist, d. h. wenn die Störung von einem Punkt nahe der Grenzlinie 9 aus- 
geht, so verlaufen einfallende und reflektierte Stoßfront im Medium (1) dicht beieinander und 
lassen sich im Grenzfall k—-0 als einheitliche Stoßfront auffassen. An der Grenzlinie % kompen- 
sieren sich ihre Drucksprünge dabei gerade. Schon hieraus kann man schließen, daß für k—0 
dieSchmidtsche Kopfwelle in der Nachbarschaft von # verschwindet. 
Zur genaueren Diskussion der Sprungintensität L,(a,r) in ihrer Abhängigkeit von « 
(bei festem 7) wird der Parameter & = cos «a eingeführt. Die Sprungintensität wird wieder mit 
L,(£, t) bezeichnet. Nach (43,) findet man dann 


LE)= 3,5," 4% 010, A@), (KzE<i)...2..) 
7 Cı 
wobei 
& &2 1 m2 
AN=— ger N ... (46 
(ö,& +0, Ym? + &)? E+mS a 
mit den Abkürzungen 
m"®—=n’—l1, en. 


gesetzt ist. Die gebrochene Stoßfront besitzt nach (34,) die Darstellung 


el, rl) +n-Errze aeecy Ws. 
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Die Ableitung dA/dE von A stimmt bis auf einen negativen Faktor mit dem Ausdruck 
BO-E—ZER—ZmES- ZURTMFER......() 


überein. Daraus erkennt man vorläufig die Gültigkeit der Feigertlen Beziehungen 


Aa) Kr 0 für große K(<I) 
20: für kleine „Kay Imınia Bil Ed I 007 ee (50) 
AR) >0 


L, nimmt also in jedem Falle nach der Diskontinuitätslinie zu (d.h. für &—K) ab. Vielleicht ist 
es wünschenswert, ein Zahlenbeispiel vor Augen zu haben. Ein solches ist für den axialsymme- 
trischen Fall durchgerechnet; in den Bildern 6 bis 9 sind die Stoßfronten 
und die zugehörigen Intensitätsverteilungen A(£) für zwei Werte von K 
dargestellt. 

Besonders interessant ist der Fall, daß die Störung praktisch in der 
Diskontinuitätslinie beginnt, d.h. daß % sehr klein ist. Man darf allerdings 
in den Formeln nicht einfach k = 0 setzen, weil dann die Ausdrücke, die für 
a-Werte mit k/c,rcosa=1 (d.h. für &-Werte mit &=K) ins Gewicht 
fallen, gestrichen würden. Vielmehr wird es notwendig sein, die 
interessierenden Größen nach geeigneten Potenzen von K zu entwickeln. 
ı wird dabei als konstant angesehen. 


Zunächst zur Form der Stoßfront (48): Für WerteE>K hat sie die Bild 5 
Darstellung 
a = 2 ym 24, y-fıe Regen Be De: 
sie ist folglich ein Kreis vom Radius c, r. Für &-Werte in der Größenordnung von K findet man 
dagegen 
K\ m Kar K 
= 021-5) a Fu =. l1- 2). +22} es Zaf (52), 


also die Darstellung einer Geraden mit dem we ı (vgl. Bild 5) 


il een Yd—a —E=tgy. 


Man erkennt deutlich den von den schönen v. Sc h midtschen Aufnahmen her bekannten 
Verlauf der Stoßfront. Die geradlinigen Teile sind die Schmidt schen Kopfwellen. 

Nun soll die Sprungintensität Z, diskutiert werden. Aus den Beziehungen (50) ergibt sich, 
daß Z, zunächst wächst, wenn man von der Mitte der Stoßfront aus nach deren Seiten zu fort- 
schreitet, dann aber wieder abnimmt. Genauere Auskunft liefert die Funktion B(£) aus (49). 


Wegen B(1)>0, B(K)< 0 hat die Gl. B(&) = 0 mindestens eine Wurzel & — E im Intervall 
K<E£öx<l. Dabei gelten für jede Wurzel & = £ jedenfalls die Ungleichungen 


u _ 1/4 
a me) an << mat +24, ma) En... 68 
2 0, 2 03 2 


Aus der Betrachtung von B’(£) erkennt man, daß es bei hinreichend kleinen Werten von K 
auch nur eine derartige Wurzel £—= £& gibt. Setzt man eine Puiseux-Entwicklung 


Bee EB SRiN ee), .,.. . do ad > o2r(54) 


an, so findet man 


1/3 0, = 
Zi m> Kıa 1 sya&» 02 mr... 
2 
& ge 


0 
hörige Punkt z, % berechnet sich zu 


les 3% „ mRUR 4.) 


u Tat 1/3 De | 
= WO = ar. Val K!r 


. (56). 


——— 


ne = 
1 
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Als Grenzlage für K— 0 ergibt sich also der Punkt 


28 1 5 | 
er a ger (57). 
2 . . . 
Das ist gerade der Punkt P,, in dem die Tangente durch P,: (0, c, ) den Kreis mit dem Radius 
c, um den Ursprung (0, 0) berührt. | j 
Für sehr kleine X kann man an Hand des Bildes 9 also folgendes sagen: Die Sprung- 
intensität wächst von P, bis P, und nimmt von P, bis P, ab. Zur Übersicht über die Größen- 


ordnungsverhältnisse mögen die folgenden Angaben dienen: r 


1 1 N 1 1 
ru 3 ur .—. BCRUNE A(K —— '— £Ä?rt:.- =0 (58). 
Az f ı a, te) legen re 
pe 
n 01 R 
Für das vielleicht einer Messung zugängliche Verhältnis Va /Y A(l) ergibt sich 
Drucksprung bei P, YAad) 1.0: 
a Er ae a I EEE IE EIER FINE DIN: 
Drucksprung bei P, Y A() n 0, 


Darüber hinaus erkennt man nun aber genauer noch folgendes: Durch Auflösung der 
Formel (48,) ergibt sich der zu einem Punkt der Schmidt schen Kopfwelle mit der Abszisse x 
gehörige Parameterwert in der Form 


1 m 
une hm (0 <a<ar:”) EEE (60). 
LE.‘ zZ 
M 5 T 
Die Sprungintensität an dieser Stelle bestimmt sich aus 


ME) = zum ea, s 


Für den Grenzfall X—0, bei welchem die punktförmige Störung von einem Punkt der Grenz- 
linie ausgeht, treten die Schmidtschen Kopfwellen somit überhaupt nicht auf. Die Vertei- 
lung des Drucksprunges auf der gebrochenen Stoßfront ist unstetig, vgl. Bild 9. Wie schon in 
der Einleitung erwähnt, muß dieses Resultat nicht als Widerspruch zu den Experimenten auf- 
gefaßt werden. 
Der axialsymmetrische Fall 
Die für den ebenen Fall durchgeführten Betrachtungen lassen sich ohne Schwierigkeiten 
auf den praktisch wichtigeren axialsymmetrischen Fall übertragen. Vermöge 
, Y=rcos7%, ee Peg . (62) 
werden die Koordinaten x, r, y eingeführt; aber alle auftretenden Größen sollen unabhängig von 


y sein. Die partielle Differentialgleichung für eine charakteristische Fläche heißt auf Grund 
einer zur früheren ganz analogen Überlegung 


| 
Han )=05 [rasjee 7 TEEN ELITE (63). 
Dagegen tritt in der Differentialgleichung für @ jetzt ein Glied mit einer ersten Ableitung auf: 
1 1 
Past Pr Tape | ER ee en (64). 


Dieses muß bei der Aufstellung der Transportgleichung für die Sprungintensität berücksichtigt 
werden. Wir geben sogleich das Resultat an. Die Gleichungen für die Stoßfronten können aus 
(34) übernommen werden, wenn dort y durch r ersetzt wird. Für die Verteilung der Leistungs- 
dichten ergibt sich 

€ 


L,(&, T) = And für cos«a < ee 
€ N k 
L,(&, 7) = And R*(&) für cos«a > et 
2 
n 1— sin? a (65) 
2 . . . 
Dina _ ‚SQ ge 
0%, 7) aan G°(6) — 


2 2 2 

2 cos? + & (1 2) = a £ ie =) 

2 C1TCOS & c)la c,Tcosa di 
- 


ei k 
für cos > —— 
ec T 


” 


sodaß“ 


Pt 


zu setzen ist. 


7o 


TOR 177) SEEN l rar 2 
= En ar ENT alte ea a ! 
Da: auf S.369 erwähnte und den Bildern 6 bis 9 zugrunde 


n Verhältnisse durchgeführt: I "Eur, hob maih 
Medium (1): NaCl-Lösg., 6, —=1,l5gfem’, ec, — 1650 mjsec; 
Medium (2): Xylol), &,=0,86gjem’, c, — 1165 m/sec, 


und: m2— 1,005 " 


Ale) für K-01 


Bild 6. Die Stoßfronten für X = 0,1 Bild 7. Die Intensitätsverteilung auf der gebrochenen 
z Stoßfront für ‚1 


ee | Ale) für K-00 
k = 000033 “= E00 
Bild 8. Die Stoßfronten für X = 0,0001 Bild 9. Die Intensitätsverteilung auf der gebrochenen 
= er ; Stoßfront für X = 0,0001 
Wiederum erhält man insbesondere 

1 1 ei! 

Al)=-— — — —'o 

= een di BE 
n 0ı 
A 1 5 

Ai 17 EUR VEN LOS) 
01 
1 1 


A ee VREHMNZEN, 
a Dee 
11) Nach einer Bemerkung von W. Schaaffs [10], S. 135, ist die hier von O. v. Schmidt über- 


nommene Geschwindigkeit c,für Xylolnicht zutreffend. Dasist aber für die hier nur derÜbersicht wegen gegebenen 
zahlenmäßigen Resultate ohne Belang. 
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also 


Drucksprung bei P, Yad ” ( 1 1% 
Drucksprung bei PR, yAd) 0 


4 te 2.) 


Beispiel: Die Totalreflexion einer ebenen Stoßfront 


Im Falle der Totalreflexion, bei dem schon das Huy ghen s sche Prinzip nicht anwendbar 
ist, versagt die Methode der Unstetigkeiten. Die gebrochene Stoßfront, welche man als Lösung 
der Eikonalgleichung zu bestimmen hätte, würde nicht mehr reell sein. Bei der Reflexion einer 
Unstetigkeit treten jetzt Singularitäten auf. Auch das war bei den bislang zugrundegelegten 
Annahmen ausgeschlossen. Wir wollen die zu erwartenden Verhältnisse an einem einfachen 
Beispiel diskutieren !2). Dazu betrachten wir die Reflexion und Brechung einer stufenförmigen 

ebenen Druckwelle an der ebenen Grenzfläche zweier Gase mit den Dich- 
A ten o, und o, und den Schallgeschwindigkeiten c, und c, (> c,). Alle auf- 
tretenden Größen werden von derz-Koordinate unabhängig angenommen. 
Das Druckfeld der einfallenden Störung im Gas (1) wird mit p, bezeichnet, 
entsprechend sind p, und p, aufzufassen. Für p, nehmen wir, unseren Ab- 

* sichten entsprechend, den folgenden stufenförmigen Verlauf an: 


27 a 1? für ae, | 
1 
, Pr= 4 ee EN) 
ee | ee er. _ı a 
1 


(vgl. Bild 11). Es handelt sich um die Berechnung von 9, und 7,. 
Der Kürze halber sollen sogleich die Resultate angegeben werden. 9 sei größer als der 
Grenzwinkel der Totalreflexion (d.h. sin®> c,/c,). Werden die Abkürzungen eingeführt: 


— 
rad, — V2sing9— 1, 
02 €8 c (71) 
_ e0s?p9—r?W? .2rW cosp o_ _2e0s?p [ i 
cos? + r2 W?’  c0s?p + r2 W?’ cp +rW: ] 
so erhält man für p, und 7,: 
a4 DI Au ) \ 
ee ARE für || < a ” 
„= 
B a+E a EN 
Ar eh Po 
und 
= = en (a 9%7 
en C C2 | Ca ) B aw“ 
9 = ze. we-9N Tr Artg m ert9 37108 RT ER (73) 
K ICH 2 w3 (7-8) 4 
mit Ee=t— ysinp h 
cı 


Der Druck der reflektierten Störung wird über den Geraden 


un, x c0osp : ySny Ines 
1 

unendlich. Diese beiden Geraden schneiden in jedem Augenblick die y-Achse gerade in den 
Punkten, in denen sich die Unstetigkeit des einfallenden Stoßes befindet. Der Einfluß der 
Störung im zweiten Medium ist jederzeit an allen Stellen des zweiten Mediums vorhanden, wenn 
er auch für große x, y-Werte stark abklingt. 

Der Anschaulichkeit halber sollen die Resultate für ein Zahlenbeispiel zeichnerisch dar- 
gestellt werden. Wir setzen etwa «— 1 und wählen das System Luft (1)-Wasserstoff (2). Bei 
760 mm Hg und 20° C hat man: c, — 343 m/sec, c, = 1305 m/sec und o, = 1,205 - 10-3 g/em?, 


er 12) Man vgl. die Arbeit von F. A. Fischer [4], wo die Totalreflexion einer „Öö-förmigen‘“ Störung behan- 
elt ist. 


Z. angew. Math. Mech. oe, 1 48 ee a RES . e 
Bd. 34 Nr. 10/11 Okt. Nov. 195& Nitsche, Über die Fortpflanzung kleiner Störungen in flüssigen Medien 373 


9 = 8,38 L 10” g/em?. Der Winkel der Totalreflexion beträgt 15,25°. Für wird 20° gewählt. 
Bild 11 zeigt den Verlauf der Druckstörung p,, Bild 12 den von 9, und Bild 13 den von p,, alle 
berechnet zur Zeit = 0. Diese Druckberge muß man sich nun noch mit der Geschwindigkeit 


c,/sin p in der positiven y-Richtung fortschreitend denken. Den Gesamtdruckverlauf würde man 
durch Überlagerung von p,, p, und p, erhalten. 


Was die Frage der Energieverteilung betrifft, so 
y kann man zeigen, daß, sofern man den gesamten Zeitraum 


— 00< t<.oo ins Auge faßt, keine Energie in das zweite 
Medium übertritt. 


Schlußbemerkung: Die Methode der Unstetigkeiten, 
angewendet auf Refraktionsprobleme, versagt in den 
Fällen, in denen die geometrische Durchführung, d.h. 
das Huyghenssche Prinzip versagt, in denen also 
eine (oder mehrere) der zur Erfüllung der Grenzbedingun- 
gen notwendigen Stoßfronten nicht mehr reell ist 23). 


/ 
7 
ZZ 


04 


Bild 12. Die reflektierte Druckstörung Bild 13. Die gebrochene Druckstörung 
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13) Die offenbar einzige diesbezügliche Bemerkung in der Literatur macht E. B. Christoffel [2], 
S. 243, der darauf hinweist, daß man in diesen Fällen wahrscheinlich auf die Ausgangsdifferentialgleichungen 
zurückgehen muß. 


Eingegangen am 13. April 1953. 
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Der höchstmögliche Auftrieb von Tragflügeln 
Von Karl Nickel in Sierras de Cördoba 


Es wird untersucht, um welchen Betrag der Höchstauftrieb eines Tragflügels mindestens absinkt, wenn der 
Flügel außer Auftrieb noch aerodynamische Momente erzeugen soll. Beispiele. Mit den entwickelten Hüfsmitteln 
lassen sich noch andere Aufgaben behandeln. 


It is investigated, how much the maximum lift force of a wing is at least diminished, if the wing effects aerody- 
namic moments, besides thelift. Examples are given. The results are applicable to other questions. " 


Il est examine, de combien le maximum de la force de montee d’une aile portante est diminue au minimum, 
si, outre la moni£e, elle doit produire encore des moments aerodynamiques. Des exemples sent donnes. Les 
moyens developpes sont applicables encore a d’autres problemes. / 


Mecsenyerca, HACKOABKO MAKCHMANBHAA TIONbBEMHAA CHA KpEla IIO MeHbIIeÄ Mepe IIOHH- 
3HTCH, ECJIH OHO, IIOMHMO NONBEMHOÄ CHJIEI, MOLKHO EINE IIPOUSBONHTb AUHAMHYECKHE MOMEHTEI. 
Ilynmepsi. PesysbTaTtt IIPHMeHHMEI U K APyTuM s3anayaM. 


1. Einleitung 

Über den Höchstauftrieb von Tragflügeln sind schon viele experimentelle und theore- 
tische Untersuchungen bekannt geworden. Dabei wird gewöhnlich die Aufgabe behandelt, zu 
einem gegebenen Tragflügel den Maximalauftrieb zu bestimmen und zu untersuchen, wie sich 
Flügelanbauten, Klappen, Absaugung und ähnliche Änderungen auf ihn auswirken. Da die 
rechnerische Behandlung dieses Problems jedoch wegen der großen auftretenden Schwierig- 
keiten über einige Ansätze noch nicht hinausgekommen ist, kann auf das Experiment vermut- 
lich noch auf lange hinaus nicht verzichtet werden. 

Die nachfolgenden theoretischen Untersuchungen verfolgen im Gegensatz dazu nicht das 
Ziel, den Zusammenhängen etwa zwischen Flügelgestalt und Maximalauftrieb nachzugehen, 
sondern behandeln eine andere Frage aus der Flugmechanik. Es ist dies die Aufgabe, das 
Absinken des Maximalauftriebs festzustellen, wenn der Flügel nicht nur Auftrieb erzeugen 
soll, sondern ihm dazu noch weitere Bedingungen auferlegt werden; wie z. B., wenn ihm ein 
bestimmtes Rollmoment oder (bei Pfeilflügeln) ein vorgegebenes Längsmoment aufgezwungen 
wird. Da man im allgemeinen ein Interesse daran hat, diesen Auftriebsverlust möglichst klein 
zu machen, soll darüber hinaus an dieser Stelle gleich das Extremalproblem behandelt werden, 
den Auftrieb unter solchen Nebenbedingungen möglichst groß zu machen, oder, anders aus- 
gedrückt, die eintretende Auftriebseinbuße zum Minimum werden zu lassen. 

Selbstverständlich kann eine solche Aufgabe beim heutigen Stande der Forschung nicht 
ohne gewisse Annahmen gelöst werden, die unter Umständen nicht immer, oder sogar nur 
näherungsweise erfüllt sind. Um nun für das folgende eine möglichst einfache Grundlage zu 
haben, soll angenommen werden, daß für jeden einzelnen Flügelschnitt sein Höchst- 
(sowie Geringst-Jauftrieb bekannt sei und sich im Verlaufe der weiteren 
Untersuchungen nicht mehr ändern möge!) Diese Grenzauftriebe der 
einzelnen Profile brauchen also nicht aus ‚„‚zweidimensionalen‘‘ Experimenten 2) bekannt zu 
sein, sondern sie können an einem zu untersuchenden Tragflügel gemessen sein und daher 
durchaus auch Induktionseinflüsse der anderen Flügelteile enthalten. Die einzige Annahme 
ist, daß sie sich nicht ändern sollen. 

Diese Voraussetzung, daß die Grenzauftriebe eines Flügelschnitts bei gegebener Flügel- 
gestalt fest sein sollen, hat jedoch sicher keine strenge Gültigkeit. Denn schon allein eine 
Variation der Auftriebsverteilung — wie sie bei unverändertem Flügel etwa durch Erhöhung 
der Anblasgeschwindigkeit und damit der Reynoldsschen Zahl verursacht werden kann — 
beeinflußt durch die veränderte Druckverteilung auch die Grenzschichtvorgänge und damit 
den örtlichen Höchstauftrieb. Da sich jedoch heute wohl kaum eine bessere Annahme treffen 
läßt und von der vorliegenden in manchen Fällen, so z. B. für sehr schlanke Flügel, genügend 
genaue Ergebnisse erwartet werden können, soll die weitere Rechnung darauf aufgebaut werden. 

Da außer der obigen Voraussetzung weiter keine wesentlichen Annahmen mehr getroffen 
werden, ist man auch völlig frei von allen Tragflügeltheorien geworden. Das ist insofern ein 
Vorteil, als damit die nachfolgenden Untersuchungen ganz allgemein für beliebige Tragflügel 


1) Diese Annahme folgt aus dem Grundgedanken der Prandtlschen Tragflügeltheorie, wo die Eigen- 
schaften aller Tragflügelprofile als bekannt angesehen werden, und wo z. B. der Auftriebsanstieg als unabhängig 
vom Einbau in den Tragflügel vorausgesetzt wird. Während man sich aber bei dieser Theorie i.a. nur auf den 
linearen Teil der Profilpolare beschränkt, wird in obiger Annahme ausdrücklich darüber hinausgegangen. Die aus- 
gezeichnete Übereinstimmung der Prandtlschen Tragflügeltheorie mit dem Versuch braucht deshalb für die 
folgenden Überlegungen nicht einzutreten. 

2) Es ist wohl zu beachten, daß die Auftriebscharakteristik eines „Profils“ stets an einem endlichen Flügel 
gemessen ist und bei hohen Auftriebswerten deshalb mit der „wahren“ Auftriebskurve des Profils in zweidimen- 
sionaler Strömung nur sehr bedingt übereinzustimmen braucht. 
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gelten, ob nun deren Auftriebsverteilungen nach Theorien der ‚‚tragenden Linie‘ oder ‚‚tragen- 
den Fläche‘ oder anderen berechnet werden. Allerdings muß dafür der Nachteil in Kauf 
genommen werden (vgl. S. 377), daß es im allgemeinen nicht möglich ist, die Lösung eines 
derart voraussetzungsarm gestellten Problems sich physikalisch realisiert zu denken. Man 
kann sie nur noch als eine Grenze des Erreichbaren betrachten, die i. a. selbst gerade nicht 
mehr verwirklicht werden kann. 

Die weitere Gliederung der vorliegenden Arbeit ist folgende: In 2. wird die Aufgaben- 
stellung an vier Beispielen ausführlich erläutert, in.3. wird das Problem mathematisch gestellt 
und gelöst. Die gewonnenen Ergebnisse werden dann in 4. nocheinmal zusammengefaßt und 
die getroffenen Voraussetzungen kritisch betrachtet. Um die Ergebnisse lebendiger werden 
zu lassen, werden dann schließlich in 5. zwei der Beispiele von Punkt 2. unter speziellen An- 
nahmen bis zur Lösung durchgerechnet. 


2. Aufgabenstellung 


Ein Tragflügel erstrecke sich von y = «a bis y = b und habe die örtliche Auftriebsdichte 
ca(y)t(y) (c.(y) = örtlicher Auftriebsbeiwert, durch Schränkung, Profiländerung und andere 
Maßnahmen in weiten Grenzen veränderlich; i(y) = Flügeltiefe, wird in den folgenden Über- 
legungen als fest gegeben betrachtet). Maximal- und Minimalauftriebsbeiwerte („mas Y) und 
Camin (Y) seien bekannt und im folgenden fest, es gelte also stets 


Camin (Y)S(Y)S Camaz (4) - u nn 5 5 AH. 


d 
Dann kann der Gesamtauftrieb des Tragflügels A = [ ©.(y) t(y) dy (der unwesentliche 
a 


konstante Faktor - Staudruck - wird hier und auch bei den folgenden Integralen stets vernach- 
lässigt) für alle möglichen c,(y)-Verteilungen offenbar nur zwischen dem dadurch gegebenen 
Maximal- und Minimalauftrieb variieren: 


b b b 
Amin — / Camin (Y) t(y) dy< / (y) t(y) dys / Camar(y)t(y) dy—= Amas 


und der Maximalauftrieb A,„s, wird nur dann erhalten, wenn für alle yaus asysbgilt: 
Ca(Y) = Camaz(Y). (Diese Bedingung läßt sich u. U. durch passende Flügelschränkung errei- 
chen.) Ob jedoch dieser Maximalauftrieb A„«,; überhaupt erreicht werden kann, hängt nicht 
nur von aerodynamischen Gegebenheiten sondern auch rein mathematisch davon ab, ob die 
C„(y)-Verteilung nicht vielleicht noch durch andere Bedingungen beeinflußt wird. Einige Bei- 
spiele ?) sollen das zeigen: 


a) Die cyuma» (y)-Verteilung eines Tragflügels sei unsymmetrisch zur Flügelmitte 
(z. B. durch einseitigen Motorausfall, Eisansatz o. ä.), und das Rollmoment dieser Verteilung 


b 


b 
| eumae 0) [vH] ar 


möge nicht verschwinden. In diesem Falle ist also der höchsterreichbare Auftrieb Ayaz 
der nach oben für c„(Y) = Camas (Y) erreicht wird, mit der Annahme eines verschwindenden 
Rollmoments — wie es für einen stationären Geradeausflug gefordert werden muß — nicht 
verträglich. Man wird also eine c„(y)-Verteilung zu suchen haben, für die das Rollmoment 
verschwindet, die überall unterhalb von c,mas (Y) bleibt, und die gleichzeitig den Höchstauf- 
trieb möglichst wenig herabsetzt. 

In gewissem Sinne eine Erweiterung dieser Fragestellung ist die Aufgabe 


b) Ein Tragflügel sei einseitig durch ein Gewicht belastet, so daß der Flügelschwerpunkt, 
und damit auch der Angriffspunkt der Luftkraft im Geradeausflug, exzentrisch am Hebelarm 


RORON z 2 2 ay 


- b 
J«W t(y) dy 


3) Man vergleiche damit K. Nickel, Lösung eines Minimumproblems aus der Tragflügeltheorie, 
Z. angew. Math. Mech. Bd. 31 (1951) S. 72—77, wo ein anderes Extremalproblem mit ähnlichen Nebenbedingun- 
gen behandelt wird. 
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sitzen möge. Dann wird wieder — analog zu a) — im allgemeinen für den Angriffspunkt der 
Luftkraft bei c,(y)= Cama, (y) gelten 


b 
Re) (Y) lv nn = as 
gi = Fe, 
N Camaz (Y) t(y) dy 


und man hat das Extremalproblem, den Auftrieb A möglichst groß zu machen, wobei aber die 
Bedingung (2) eingehalten werden muß und gleichzeitig noch c„(y) durch die Ungleichung (1) 
eingeschränkt ist. PR 
c) Bei einem Pfeilflügel nach Bild 1, bei dem die Angriffspunkte der örtlichen Luft- 
kräfte von der als Bezugsgeraden gewählten y-Achse den Abstand A(y) haben, ändert sich bei 
Veränderung der c„(Y)- Verteilung offen- 
= 6 y sichtlich auch das Längsmoment 


Sen anWay - - - @) 


Um bei einem Schwanzflugzeug nicht 

zu starke Höhenleitwerksbelastungen 

zu bekommen, wird man für das Längs- 

moment Grenzen angeben müssen, inner- 

halb deren eine Änderung noch tragbar 

Bild 1. Pfeilflügel ist. Noch enger ist man bei schwanz- 

= losen Flugzeugen gebunden, bei denen 

mit der Schwerpunktlage im stationären Geradeausflug auch automatisch die Lage des 


Druckmittelpunkts ; 


[ca(y) t(y) h(y) dy 
Bm ee ET 
KO) t(y) dy 


festliegt. Zwar hängt die Linie A (y) noch von den Profildaten ab und kann daher (etwa durch 
Profiländerung oder Flügelklappen) noch etwas verschoben werden. Doch handelt es sich bei 
den üblichen Tragflügeln meist um langgestreckte Gebilde, bei denen diese Veränderung im 
Verhältnis zur Größe von h(y) wenig ausmacht, so daß man sie also etwa als Korrekturgröße 
auffassen kann und im übrigen mit festem Hebelarm A (y) rechnet. 

Auch hier hat man, wie bei a) und b) eine Nebenbedingung, die für c,(%Y) = Cgmaz (Y) 
im allgemeinen nicht erfüllt sein wird. Man wird auch hier wieder nach derjenigen c„(%)- 
Verteilung fragen, die den Höchstauftrieb möglichst wenig absinken läßt, während anderer- 
seits das Längsmoment oder die Druckpunktlage einen vorgeschriebenen Wert behält. Nach 
Lö:ung des Problems kann man sich dann die Änderung des Auftriebs z. B. über der Druck- 
punktlage auftragen und damit die ‚„‚günstigste‘“ Schwerpunktlage feststellen. 


+5 


Wählt man insbesondere die Funktion h(y) = y— ee ‚ so erweisen sich a) und b), 
mathematisch gesehen, als Sonderfälle von c). 2 


d) In den vorhergehenden Beispielen war ein in c,(y) lineares Integral, der Auftrieb, 
unter einer gewissen Nebenbedingung zum Extremwert zu machen. Dieses Integral braucht 
aber nicht unbedingt den Auftrieb zu bedeuten. Z. B. kann man auch in c) das Längsmoment 
(3) zum Maximum zu machen suchen und dabei den Auftrieb A als gegeben vorschreiben. 
Diese Frage tritt z. B. auf, wenn man für ein schwanzloses Flugzeug diejenige Auftriebsvertei- 
lung sucht, die eine möglichst große Höhenruderwirkung besitzt. 

Offenbar stimmt auch diese Problemstellung mathematisch völlig überein mit den vor- 
hergehenden, denn ersetzt man in c) die Tiefe t(y) durch eine fiktive Flügeltiefe t(y) h( Yy) und 
die Funktion h(y) entsprechend durch 1/h(y), so geht c) in die obige Frage über. | 

Im vorangehenden war c,() immer durch die Ungleichung 


Cmin(Y)S C(Y)S Camaz(Y) - eine Arne Feie Eee Star BL) 


eingegrenzt. Für die mathematische Fragestellung der Ziffer 3., die alle diese Beispiele zusam- 


menfaßt, ist es nun natürlich belanglos, daß c„(y) in (1) gerade durch seine physikalischen 
Grenzwerte eingeschränkt wird. Ändert man (1) um in 


N) a(Y)< 0wr(Y) ee) ER in sinn: (la), 


Tr 


N 
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wobei c,,(%) und c„,(4) willkürliche c,(y)-Verteilungen sind, so bleiben sämtliche zu ziehenden 
mathematischen Folgerungen und damit auch alle Ergebnisse erhalten, wenn man mm und 
Camas beziehentlich durch «,, und c„, ersetzt, da in den folgenden Überlegungen von der 
physikalischen Bedeutung von Camaz BZW. C,mim Kein Gebrauch gemacht wird. (Es möge 
angemerkt werden, daß mit dieser Erweiterung der Fragestellung u. U. die in der Einleitung 
erwähnten Bedenken entfallen können, die mit der Annahme (1) verknüpft sind.) 


. Z.B. besitzen manche Flügelprofile einen gewissen „‚Laminarbereich“ in c,, innerhalb 
dessen sie besonders niedrige Widerstandsbeiwerte aufweisen. Will man dann den Gesamtauf- 
trieb des Tragflügels bei kleinem Widerstand groß machen, so kann man in den vorhergehenden 
Beispielen die Forderung (1) durch (la) ersetzen, worin &, (Y) und &; (Y) die Grenzen dieses 
Laminarbereichs sind. 


Unter Vorwegnahme eines Teils der Ergebnisse der Abschnitte 3. und 4. soll hier schon die 
Lösung der Beispiele a) bis d) angegeben werden: Die günstigste c, (y)-Verteilung, die also den 
Auftrieb A zum Maximum macht unter gleichzeitiger Erfüllung der Beschränkung (1) und 
einer der weiteren Nebenbedingungen, setzt sich stückweise aus &min(Y) und C,mas(Y) 
zusammen (siehe Bild 2). Über- 
raschend ist, daß also die Lösung nur 
insofern von der Integral-Neben- 
bedingung abhängt (z.B. von der 
Funktion h(y) in Beispiel c)), als 
durch diese die Sprungstelle im 
Ca (y)-Verlauf festgelegt wird, daß sie 
in ihrem Charakter jedoch von die- 
sen Nebenbedingungen gar nicht be- 
einflußt wird. Ebenfalls nicht von 
vornherein zu erwarten ist die Tat- hi 
sache, daß diese Lösung für beide 


M Ca(y) 
a max Y) als 


on _ 


ne nn 


—— 
——— 
— 


Arten von Nebenbedingungen gilt, 
wie sie durch (3) und (4) vertreten 
werden: konstantes Längsmoment 


Ca min (J) 


Bild 2. Günstigste [7 (y)-Verteilung im Beispiel 2 a). (Unsymmetrische c, mar. 
Verteilung entspricht etwa einseitig blockiertem Vorflügel.) 


und konstanten Druckpunktsabstand. 


Da damit die c,(y) t(y)-Verteilung ebenfalls Sprünge bekommen und weiter an den 
Flügelenden meistens auch nicht verschwinden wird, ist die so gefundene Funktion c,(Y) t(y) 
im allgemeinen keine mögliche Auftriebsverteilung des Tragflügels !}Dies ist jedoch meist 
nur ein scheinbarer Nachteil, da man i.a. auch Funktionen c,(y)t(y) mit Sprüngen beliebig 
genau durch mögliche Auftriebsverteilungen (die also u. a. stetig sind) wird annähern können. 
Der Hauptzweck der vorliegenden Abhandlung ist ja auch nicht, die vorteilhafteste c,(Y)- 
Verteilung zu gewinnen, sondern den dazu gehörigen höchstmöglichen Auftrieb des Trag- 
flügels zu erhalten. In diesem hat man dann eine obere Schranke für den unter den gegebenen 
Nebenbedingungen möglichen Auftrieb gefunden und weiß weiter, daß man ihm durch passende 
Auftriebsverteilungen beliebig nahe kommen kann, daß er also sogar eine obere Grenze 
darstellt. 


Damit ist ein Sachverhalt gekennzeichnet, wie man ihn bei Extremalproblemen auch 
sonst nicht selten antrifft, daß nämlich alle zulässigen Auftriebsverteilungen einen Auftrieb 
besitzen, der einem Grenzauftrieb beliebig nahe kommt. Die zu diesem Grenzauftrieb gehörige 
Auftriebsverteilung selbst ist jedoch nicht mehr zulässig. Durch diese Überlegungen erkennt 
man also jetzt, daß es nicht nur unnötig ist, eine weitere Bedingung hinzuzunehmen, durch die 
alle physikalisch möglichen Auftriebsverteilungen ausgewählt werden, sondern daß unter solch 
einer Forderung die formulierten Extremalprobleme gar keine Lösung mehr besitzen würden. 
Es hat sich also nachträglich herausgestellt, daß schon die geringen getroffenen Voraussetzun- 
gen der Seite 374 die Grenze des Zulässigen darstellen, wenn die an den Beispielen erläuterte 
Aufgabe überhaupt eine Lösung haben soll. 


3. Mathematische Problemstellung 


Zur Vereinfachung der Schreibweise werde y(y) = c,(y) t(y) und c(Y) = & min(Y) t(Y) , 
sowie ((Yy) = Czmas(Yy) t(y) gesetzt. Da nur positive Flügeltiefen i(y) einen physikalischen 


Sinn haben, verwandelt sich damit die Ungleichung (1) in die entsprechende 
c(y)sry)=C(y) (ib). 
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Man trifft nun folgende 
Voraussetzungen: Die Funktionen c(y), C(y) und A(y) seien auf der Strecke 
a<y<b definiert und dort stückweise stetig mit c(y)<C(y). Darüber hinaus sei h(y) 


noch stückweise differenzierbar mit rn. s O(h(y) monoton wachsend oder fallend). 


Es sollen die beiden folgenden Extremalprobleme gelöst werden: 


I. Die reelle Zahl M sei beliebig gegeben. Dann soll y(y) als stückweise stetige Funktion 
ina<y<bso bestimmt werden, daß 


b 
A) = Sy) dy = Max sit amrnbeo ah alu 
wird unter den beiden Nebenbedingungen 
b 
AL a Ed en EEE IRRE 
und 
NWEYW)SCWane 2 Zr BE a 
II. Die reelle Zahl H sei beliebig gegeben und zusätzlich zu oben sei noch zur Verein- 
fachung gesetzt c(y) = — KC(y), wo K eine beliebige positive Konstante bedeuten möge. 
Dann soll y(y) als stückweise stetige Funktion ina<y<bso bestimmt werden, daß 
b 
Air) = | v(WW) dy — Max. Neger, 2 EN 


wird unter den beiden Nebenbedingungen 


d 
h(y) dy 
m Sr@Wht 
Be — 4 enger in 
BRAOLE 
und 
au) el)... 2 Ne ne ee 

Es wird zuerst die Maximumaufgabe I. gelöst und zwar sind drei Fälle zu beachten, je 
nachdem A(y)>0, h(y)<0, oder A(y)S0 in a<y<b ist. Anschließend wird dann 
nachgewiesen, daß unter der Einschränkung c(y)—= — KC(y) die Gesamtheit der Lösungen von 
II. mit der von I. übereinstimmt, oder, genauer ausgedrückt, daß es zu jedem H’ gerade ein 
M’ gibt derart, daß die Lösung von I. für M’ gerade auch Lösung von II. für H’ist. Es wird 
dabei gezeigt, daß I. und II. immer genau eine eindeutig bestimmte Lösung haben, falls nur M 
aus einem gewissen Bereich gewählt wird. Für H ist dagegen jeder beliebige Wert zulässig. 


I. 


a) Seih(y)>20 und dh(y)/dy>0 in a<y<sb. 
Behauptung: Dann hat I. für alle M mit 


b b 
M, =[e(y)h(y) dys M<.[,O(y) R(y)dule U ee 
genau eine stückweise stetige Lösung y, und diese ist gegeben durch 
_ fly) für asy<p 
OR für Nil 3 en 


Dabei ist 7 durch die Nebenbedingung (6) M(y,)—= M im Intervall a<&n<b eindeutig 
bestimmt. 

‚Beweis: Offenbar wird die Nebenbedingung (1b) durch y, aus (9a) immer erfüllt, 
gleichgültig, wie 7 gewählt wird. Weiter ist sicher immer M, < M(y,) < M,, wobei die Grenzen 


für n— a bzw. n = b tatsächlich angenommen werden. Also bedeutet (8a) keine Einschrän- 
kung der Lösungsmöglichkeit. 


Es soll nun n durch (6) festgelegt werden. Es gilt 


n b 
M(y,) m C(y) h(y) dy+ f c(y) h(y) dy 
und also für n>a { 


dM 
an —kmlCmM)—enl>o. 


| 
. 
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Da nun noch M(y,)—= M,, M(y,) = M, ist, gilt also: Wächst n von « bis b, so wächst M(y,) 
monoton von M, bis M,. Zu jedem M mit M,< M< M, gibt es also genau ein n derart, daß 
M (Yn) —=ıM ist. 

Um nun die Maximumseigenschaft von y, nachzuweisen, setzt man f(y) = y,(%) — y(Y)» 
worin die Vergleichsfunktion y(y) stückweise stetig sein soll und die Nebenbedingungen (6) 
und (1b) erfüllen möge. Nach (1b) ist dann 


el at. Rz | 
IN <o ee es: (10a), 


und aus (6) folgt einerseits durch partielle Integration, andererseits durch Einselzen von f(y) 


b b b 
un=ie) [san TR Irwaay=0...... - (1a) 


Sei nun speziell 7(y) solch eine Vergleichsfunktion mit A(#)> A(y,), oder A(f)<0 
= d_ 
(f= 7» — 9). Dann ist nach (10a) auch [f(t)dt< 0, die beiden Summanden in (11a) haben 


v a 
also beide dasselbe Vorzeichen (<0). Also ist (11a) nur für f(y) = 0 erfüllbar, d.h. aus der 
Annahme A(y) > A(y,) folgt y=y,. Also ist y„(y) die einzige Lösung des Maximum- 
problems 1. 


b) Sei A(y) <O und dh(y)Ydy>0 in asy<sb. 
Behauptung: Wie bei.a), jedoch sind (8a) und (9a) zu ersetzen durch 


b 
M, = Ch) dy< M< fehl) We MM, atge 3 ae (8b) 
und 
eh EHÜLTASTHENE 
7a) se rreryhahrnden wirhigch ne (9). 


Beweis: Ganz entsprechend zu a), nur ergibt sich hier statt (10a) 


swlZ, = en a (10b) 
und statt (11a) ? 1 
y 
Mf)=h0b) | f)ay — | nn j dr ee (11b). 


ec) Sei AK) —=0 (a<d<b) und dh(y/dy>0 in asy<sb. (Wechselt h(y) bei 
ö das Vorzeichen, ohne daß h(9) = 0 ist, so kann man trotzdem ohne Einschränkung der 
Allgemeinheit k (9) — 0 setzen, da es ja auf den Wert von k(y) an der Stelle # für das betrach- 
tete Problem nicht ankommt.) 


Behauptung: Wie bei a), jedoch sind (8a) und (9a) zu ersetzen durch 


d b E d 
M=SCcWhmdy+SedhWdysMsewhwdy+[OoW)hWdy= M, (8) 
und 
C(y) für. asy<sn, RN 
A ET Me nn 
Yn\Y) = FeEn<D) - Seas ge). 
"1% c(y) für asysn, A 2, ist 
RE 22 Para 
Für n = ® ist y, doppelt bestimmt, entsprechend den beiden Fällen n„2# und ns #, doch 
wird dadurch wohl kaum eine begriffliche Schwierigkeit entstehen. 


Beweis: Ganz entsprechend zu a). (10a) bleibt im Falle „> % gültig und wird für 
n< ® durch (10b) ersetzt. (11a) ist etwas zu modifizieren und lautet jetzt 


DI b b 
un=|sonwar + | 2 |roaay tür 20 


bzw. (11e), 
M(f) = IKKOLXO dy— | za fl)didy für n< | 
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d) Sei dh(y)/dy<0 inas<y<b. Dann gelten ganz entsprechende Behauptungen 
und auch die Beweismethode kann unverändert übernommen werden, 
II. 


€) Behauptung : Zu jedem 4 mit = o<H<-® gibt es ein 9 derart, daß 
die Lösung y, von I. auch die Nebenbedingung (7) H (y,)—= H erfüllt und daß y, auch einzige 
Lösung von II. ist. ® 

Beweis: Sei A(yJ)=0 und dh(y)/dy>O in asy<sb. Aus der Voraussetzung 


c(y)=—KC(y) folgt c(yJ)<0<C(y). Dann ist mit y„(y) aus (9a) 
n b 
A= [CW)ay+jewäy, aso AuW=Jewdy<o, 


b i 
Am=| CWay>0 una Tr Cm) en)>0. 


Weiter gilt 


[own ay + Joa) ay Jen ay 
ee 7 ‚ also H(y.)= H(p) = —— 
JeW4W ie) e(y) dy Sewdy 


und 


dH C(n)— , b 
EN Fomıne— am ay+ fewWiR@— Alan >o. 
dn [A (yn)] a u) 
Nun ist A (y,) stetig in n, nimmt also den Wert Null genau einmal an. Sei dies an der Stelle n, 
dann findet man durch partielle Integration 
d 


b Y 

= dh 

1m | von — [HR | miwaray <o, 

‘a a a 

M(y„) verschwindet also an dieser Stelle nicht. Da H(y,) aber auch stetig in » ist, nimmt 

H (y„) also — wenn man etwa n = a ausschließt — jeden Wert aus (— ©, + oo) genau einmal 
an *). Also kann man zu vorgeschriebenem H genau ein n so finden, daß 4(y,) = H ist. 
Für k(y)<0,h(y)=0, sowie dh(y)/dy<O in a<y<b schließt man genau so. 

Man denkt sich nun zu gegebenem A die Funktion y„(y) nach oben so bestimmt, daß die 

Bedingung (7) erfüllt ist. Sei weiter Y(y) eine Vergleichsfunktion, die also stückweise stetig ist 

und die Bedingungen (7) und (1b) erfüllt. Ist nun A(y) > A(y,), so setzt man (4A(Y) #0 


wegen |4H(Y)| < ) g 
aD N RE 


IYyV> 


Offensichtlich erfüllt auch 7(y) noch die Nebenbedingungen (7) und (1b), wegen 


49) = ie AB) = Ay) 


ist aber auch 


M()= HA(y)= M(y,); 


d.h., y(y) ist auch eine Vergleichsfunktion für das Problem I. mit M — M(y,). Da jedoch I. 
die eindeutig bestimmte Lösung y„(y) besitzt, gilt also 7— y„, oder nach (12) ebenfalls 
Y= m. Also ist y„(y) auch Lösung der Aufgabe Il. 


4. Ergebnis 
Längs des Flügels von y— a bis y= bsseien die Flügeltiefe t(y) > 0°), sowie die beiden 
örtlichen Grenzauftriebsbeiwerte Comin (4) und- vCgmas (de, min Camas) als stückweise 
stetige Funktionen gegeben. Ferner sei der Hebelarm h(y) dort eine stückweise stetige, stück- 
weise stetig differenzierbare und monotone Funktion. 
4) Die Stelle H(y) = + oo hat offenbar keinen physikalischen Sinn, da ihr immer der Auftrieb A(y) = 0 


zugeordnet ist, man also nicht gleichzeitig noch A(y) = Max. verlangen kann. 
5) An den Flügelenden y = a und y = b kann t(y) = 0 zugelassen werden, 


En 


| 


siar iu link ü j t ’ AzEH Aa Fe a ma RT wadsn j seh 
Er ee Areal guy) vH dq>bRil we Female Au 

N e Hr ü ] b #a1.Äl bi [| 007 zafe via. = 
AR 2-26 . : y . . hi - 5 B . ae F Bee a 
unter diesen Bedingungen genau dann seinen Maximalwert an, wenn c„(y) stückweise gleich 


= S R be } A r 1 . Pi 4 er . a 
hränkt, so nimmt der Auftrieb WINE mb 30% Klug ob ind 
la Au) yı vn | in 1dam Ianin 59] a9 Ilona 
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Camin(Y) DZW. C4 mas (Y) ist. Und zwar gilt 


f = lee | art dh & R 

00 [mas für asy<n, AU were nd 
r RE - wenn oder = ea rar 
Camn(y) für n<ysb, k(n)<0, gm Inöi 
Gy) = dh (13). 


Ca min (Y) für a zsys Nn> | 
l 
| 


a) et y<b, 


SI 


j ) n 
Die Trennstelle nin (13) ist durch die Momentenbedingung (3) eindeutig (bis auf den Fall 7 = ®) 
festgelegt. 

Die Werte M, und M, sind die Grenzwerte des Moments, die unter der Bedingung (1) 
überhaupt erreichbar sind. Sie lauten, wenn h(y) bei ® das Vorzeichen wechselt 


J Camino u) ha dy für A(y)>0, 
Male uEk ey a > 
res b dh (14a). 
| J Camaz(y) t(y) h(y) dy +] Samin (9) t(y)h(y)dy für rs 0, 


8 b 
[ann DU Amdy + [ana (EM ho)dy tür Z<0 


(14b): M, erhält man aus M,, wenn man darin cam und (mas vertauscht. 
‚Ist zusätzlich zu den obigen Voraussetzungen noch Cymin (Y) = — KCamaz (Y) (K eine 


beliebige positive Konstante), so ist die c,(4)-Verteilung aus (13) auch noch Lösung des folgen- 


b 
den Problems: Man sucht den Auftrieb A = ji Ca(Y) t(y) dy zu einem Maximum zu machen, 


wenn der Hebelarm ZH der Luftkräfte 
b 
Jay) t(y) h(y) dy 


| N En a 
Hr / (y) t(y) dy 


beliebig vorgeschrieben ist und c„(y) wieder durch die Bedingung (1) eingeschränkt wird. 


Zusatz (ohne Beweis): Sucht man den Auftrieb A unter den obigen Bedingungen 
zum Minimum zu machen, so gilt derselbe Satz wie oben, nur sind in (13) &,mas und 
Camin miteinander zu vertauschen. 


Erweiterungen (ohne Beweis): Eine Abminderung der Voraussetzungen über die 
Funktionen £(%), Camin(Y) und Cymas (y) hat wenig Bedeutung, da man mit der stückweisen 
Stetigkeit dieser Funktionen praktisch immer auskommen wird. Dasselbe trifft auch für h (y) 
bezüglich der stückweise stetigen Differenzierbarkeit zu. 


6) Mit der Verabredung h(9) = 0, falls h(y) bei 9 das Vorzeichen wechselt. 
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a) Dagegen sind durchaus praktisch vorkommende Fälle denkbar, in denen h(y) nur 
noch stückweise monoton und zwar teils monoton wachsend, teils monoton fallend ist. Für 
diese Fälle läßt sich durch Gegenbeispiele zeigen, daß die obigen Behauptungen i. a. nicht mehr 
aufrecht erhalten werden können, daß eine stückweise aus C,maz(Y) und C4min(Yy) zusammen- 
gesetzte c„-Verteilung durch die Nebenbedingungen eindeutig bestimmt ist und die Aufgabe 
eindeutig löst. Die Monotonie von h(y) ist also für das Bestehen des formulierten Ergebnisses 
notwendig. 

b) Wird bei der Aufgabe II. auf die Bedingung &umin(Y) = — K Camaz(y) verzichtet, 
so braucht es jetzt nicht mehr eindeutig möglich zu sein, zu jedem H eine c„ (y)-Verteilung mit 
diesem Druckpunkthebelarm zu finden, wie man durch einfache Gegenbeispiele nachweist. 
Eine Äquivalenz von I. und II. ist damit nicht mehr möglich. 

c) Erweitert man schließlich die Aufgabe I. (oder II.) auf mehrere Nebenbedingungen der 
Gestalt (3) (oder (4)), wie etwa durch Kombination der Beispiele b) und c) des Abschnitts 2., so 
läßt sich wieder durch Gegenbeispiele zeigen: Es ist im allgemeinen nicht möglich, «,(%) 
als eine Funktion zu bestimmen, die stückweise gleich C4min (Y) Oder Ca mas (Y) ist derart, daß 
die sämtlichen Nebenbedingungen erfüllt sind; auch wenn diese nicht mit (1) in Widerspruch 
stehen. Falls eine solche Bestimmung trotzdem möglich ist, so braucht sie weder eindeutig 
möglich zu sein, noch braucht dann der Auftrieb A zu einem Maximum zu werden. 


5. Beispiele 
a) (vgl. 2b). Es soll der Maximalauftrieb eines Tragflügels bestimmt werden, bei dem 
durch unsymmetrische Belastung der Druckmittelpunkt um die Strecke e außerhalb der Mitte 


liegt. Es sei @,mar = — Camin = Konst. angenommen”), zur Vereinfachung der Rechnung 
werde C„maz(Y) = 1; Camin (y) = — 1 gesetzt, außerdem werde die Spannweite auf 2 nor- 
-7 0 +7 7 


2 
#1) 
T 
t(yy=tio) [1- (1-Z) ıyl] (-1=y=+7) 
BEE, 
2777 


Bild 3. Trapezflügel. 


miert (diese Vereinfachungen gehen nicht in das Ergebnis ein). Als Flügelgrundriß werde 
die Schar der Trapezflügel mit variabler Zuspitzung Z =(1)/t(0) verwendet (vgl. Bild 3), da- 
bei möge Z zwischen 0 (Dreiecksflügel) und 1 (Rechtecksflügel) verändert werden. 


Mit £(y) = t(0) [1 — (1 — Z) |yl] (-1s<ys+1) und 
RER A 

“NL für n<y<-+l 
wird der Auftrieb für —1 <n<0 


a 
A= [WıWdy=— nt) +A— Zn] 
und das Rollmoment 


Dr, (0) .. 
2 Jet) yay=-5B( )-2(1— ZU +M)), 


und also 
N Dre 131 —- )—-21—-Z(i1-+N) 
(M)=4I-- et 
3 In FI=Zr 
(Aus Symmetriegründen genügt es, sich auf e> 0 zu beschränken.) 


Es soll nun die Abnahme von A mit größer werdender Exzentrizität e des Druckmittel- 
punktes festgestellt werden, wobei A durch Division mit Ayar zu A/ Amaz normiert werden 


soll. Amaz wird offenbar für e—= 0, also n= — 1lerreicht. Also ist 
} £ A 2n +1 — Zn 
Anas.=t(0) (1 - Z) Munue s u e N 
Ser. 2 ‚ a A max 1+2Z ö 


?) Entspricht etwa gleichem, symmetrischen Flügelprofil längs der Spannweite, 


k 4 


Zen 
as e naad 
N 


wi sÄ 


Dee wie: di 
| Apr) Itw_dnaaie PETTTE 


En = = n»ES% \ 
‚MH „ 1n Aka } 27 H ey u; ün dene Er | 15 Bin u d r try “ is? Pr 
ei; AG bi; fl zuryt „_ = (I Amas)? karı Ar 1 asiao ö h UDRRER TR “13 al 7:3 


_ Dieses Ergebnis ist für Z= 0, 1/10, N Bet u 
1/5, 1/3, 1/2, Lin Bild 4. aufgezicm Amar Bl: 
net. Besonders interessant ist daran, 
daß A/Amas für e>0 nicht mit 

p- waagrechter Tangente in 1 einmün- 
j det, sondern unter der Neigung 


d.A/ Amaz 


- de  ke=0 


Liegt also z. B. der Auftriebsmittel- 
punkt um 10% der Halbspannweite 
exzentrisch zum Flügelmittelpunkt, 
so ist der Auftriebsabfall im günstig- 
sten Falle ebenfalls rund 10%. 


b) (vgl. 2c). Beidem Pfeilflügel 
in Bild1, der sich ebenfalls von 
a—=—1 bis b=-+]1 erstrecken 
möge, sei der Hebelarm Ah(y) der 
Luftkräfte proportional zu |y| (gerader 
Pfeilflügel)®). Es soll der Maximal- 
auftrieb in Abhängigkeit von der 
Druckpunktlage bestimmt werden, 
wenn wieder wie bei a) 2 


=—1 (für alle Z). 9 


07 


Camas(Y) = — Camn(y) =1 
angenommen wird und lineare Tiefen- 
verteilungen, d. h. Trapezflügel, be- 


trachtet werden. 00 91 02 03 0% 05 06 a7 
F In dem vorliegenden Falle ist 


monoton sei, nicht erfüllt. Wegen der 
Symmetrie des Flügels ist es jedoch statthaft, sich auf eine Flügelhälfte O<ys+1 zu 
beschränken, und dort ist h(y) monoton. itanf 


Mit {= 1) [1 — (1—Z)|y]] una 
—1. fü 1s y<—n 
ay)=ıtr1l fü nsystn 
5 —ı fr nme 


erhält man für den Auftrieb 
1 
4= IK) t(y)dy=t()An—2+1—- ZA— 2) 
und mit k(y) =|y| für das Längsmoment 


u (ae ini Mr +20 Da 2m. 


uw. 
\ 


aa a 


8) Es möge hier ausdrücklich angemerkt werden, daß die Linie der Druckmittelpunkte eines Flügels im allge- 
meinen Fallen ic h t gleich der {/4-Linie (Linie der Profilneutralpunkte) ist. Über die genaue Form des betrachteten 
Tragflügels wird daher bei obigem Beispiel nichts ausgesagt, sie ist mit h(y) und t(y) noch nic h t festgelegt. Im- 
merhin wird man für langgestreckte Flügel wenigstens näherungsweise die Druckmittelpunktlinie mit der t/4-Linie 
identifizieren können, 


-17 oO +7 Y 


VRR EE-- 2 
| E 


b) | | hiyeborlgtigp 


Bild 5. Gerader Pfeilflügel 


Der Wert A= Amaz wird offenbar angenommen, wenn c,(Y) = Cama(Y) =1, also 
n=-+1 ist. Diese c„(y)-Verteilung bedeutet aber eine bestimmte Druckpunktlage, die mit 


H= ( bezeichnet werden möge °). Es ist C 5477 d.h. C ist nicht fest, sondern 


ändert sich noch mit der Zuspitzung Z. Damit ist es zweckmäßig geworden, statt Z als Ver- 
änderliche die Größe C — H einzuführen, also den Abstand von dem Druckpunkt (, in dem 
der maximale Auftrieb erreicht werden kann (C — H> 0 ‚‚kopflastig“, 0 — H < 0 ‚„‚schwanz- 
lastig“‘). Weiter ersetzt man wieder wie oben A durch A/Amaz. Man findet wieder 


4 HR ZN 
EN. 5 1+Z h 
Drückt man nun hieraus wieder n durch A/Amas aus und setzt das Ergebnis dann in 
C — H(n) ein, so ergibt sich schließlich 


C— Hz0: 


ar) 1) 


2 A A 3/2 
= 3 |z: ( +1)— —1 firahzae ig 
/ Amaz Amsz 5 
11—- (AA) . 
ne fur, „7, 
4 Ale 
9) R.u. W.Horten, Stabilitätsbetrachtungen an Pfeilflügeln. Bericht 164/Vorbericht der Lilienthal-Ge- 

sellschaft für Luftfahrtforschung, P 37 (Geheim), S. 59—70, besonders 8. 63. Dort wird, bei etwas anderer Frage- 
stellung (Abkippverhalten), dieser „C-Punkt‘ wohl zum ersten Mal angegeben. Bedeutet h(y) den Hebelarm der 
{/4-Linie (vgl. Fußnote 8), so wird dieser Punkt neuerdings bei Stabilitätsuntersuchungen auch als „geometrischer 
Neutralpunkt‘““ bezeichnet, 


Ans =t0) (1 2), alen 


ta mas. arme 
EN 
Zi A Anne 
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Bild 6. Gerader Pfeilflügel nach Bild 5a: Abnahme des erreichbaren Höchstauftriebs mit dem Abstand des Druckpunkts Z vom C-Punkt. 
(Für € — H> Ofallen die Kurven für Z= 1/3, 1/5, 1/10 fast-mit den anderen Kurven zusammen und sind daher weggelassen.) 


Das Ergebnis ist für Z= 0, 1/10, 1/5, 1/3, 1/2, 1 in Bild 6. aufgezeichnet. Man beachte 
wieder die von Null verschiedene Neigung von A/ Amas bi O— H=0: - 
he 
AA 
lu wi 
Bgzr für OÖ—H<0. 


für C— H>0, 


Eingegangen am 8. Dezember 1952. 


Über Gleichungssysteme der Tschebyscheffschen 
= Approximation 


Von F. Wenzl in München*) 


Aufgaben der Tschebyscheffschen Approximation mit reellem Approximationsintervall lassen sich häufig 
auf Gleichungssysteme zurückführen, in denen auch die Stellen mawimaler Abweichung innerhalb des Approzi- 
mationsintervalls als Unbekannte auftreten. Die rechnerische Behandlung solcher Gleichungssysieme kann sich 
beträchtlich vereinfachen, wenn man diese Stellen nach dem Vorzeichen der zugehörigen Differenz zwischen Nähe- 
rungsfunktion und angenäherter Funktion in 2 Klassen einteilt und von jeder dieser Klassen die symmetri- 
schen Grundfunktionen aufsucht. 


Problems of the Chebyshev type, with real range of approximation, can often be reduced to systems of 
equalions. The unknowns of the systems are, among others, the points of the greatest difference within the 
range of approximation. The numerical treatment can be considerably simplified by classifying these points 
accordings to the sign of the difference beiween the given function and the approximative function, and by 
ihe determination of the elementary symmetric functions of the two classes. 


* Verfasser ist am 29. September 1953 verstorben. 
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imali i d’approximation 
&mes du type de l’approwimation de Tchebycheff avec un intervalle reel d’approx 
regt Äe Kira echt ädes eS d’&quations, dans lesquels aussi les points de plus is: a 
au dedans de l’intervalle d’approaimation figurent comme valeurs inconnues. Le maniemen Bi en 
tels systemes d’&quations peut öre simplifi& considerablement, si l’on partage ces Feng en Ei en 
suivant le signe de la difference entre la fonction donnee et la fonction approwimative, et delermine le 
symmeiriques Elömentaires de chacune de ces classes. 


3anayu anııporcnuManun YeöstmeBa c BEINCCTBEHHBIMHHTEPBA.IOM IPHÖNMKEHUA JACTO mer 
ÖbITb IIPHBENEHEI K CHCTEMAM YpaBHeHHÄ, B KOTOPEIX PurypHPpyIoT, B Bue ee BauTe 
M&KCHMA.IBHOTO OTKJIOHeHHA B IIPENENAX &IIPOKCHMHUPYIOIIerO UHTEPBA.IA. ae a 
O6pab0Tka TAKUX CHCTEM YPaBHeHHÜ CyIIecTBeHHO yUPoImaeTca, ecım Pas Er aTu toere 
110 3HAKY COOTBETCTBYIOINEH PABHOCTH Mekıy IPHONmKamımeh u IIPuÖNmKaeMmoNH DyHEI 
Ha BA Kllacca U YVCTAHOBUTL CHMMETPHYeCcKHe OCHOBHBI® PYHKIIHH KAKIOTO KIacca. 
Die Lösung Tschebyscheff£fscher Approximationsaufgaben mit reellem Approxi- 
mationsintervall — gegeben durch die Forderung, die Parameter a,in f(@;@, ;...q,)so zu 
‚bestimmen, daß | . =“ 
Max. |f| möglichst klein ... - ---..... (11) 
AsısB 
wird — führt häufig auf Bedingungen folgender Art: !) 3 


Zu der gegebenen Funktion f(# ;a, ;... a„) und den gegebenen Konstanten A und B 
sollen die Argumentwerte &, und A, sowie eine Konstante Z so bestimmt werden, daß 


Ami KERN ee 
ferner 
N&:A:.. A ezirr ame 1:2. A ee 
aa ee für kkeo2 ;.,.0 Sand, Sn ee 
lasdı 3... A für AS zB)... Haze 


Damit ist dann die gestellte Approximationsaufgabe auf das System der Gleichungen 
(1.2b) und (1.2c) mit den Nebenbedingungen {1.2a) und (1.2d) zurückgeführt. Praktisch ist 
allerdings mit einer solchen Zurückführung auf ein Gleichungssystem selbst bei rationalem f 
häufig noch nicht viel gewonnen, da man schon bei verhältnismäßig einfachen Aufgaben eine 
große Anzahl von Gleichungen höheren Grades bekommt. Im folgenden soll versucht werden, 
unter gewissen Einschränkungen für die gegebene Funktion f aus dem Gleichungssystem 
(1.2b), (1.2c) ein anderes System von Gleichungen niedrigeren Grades zu gewinnen. Der Grund- 
gedanke dafür ist, die Symmetrie-Eigenschaften des Gleichungssystems (1.2b); (1.2c) bezüg- 
lich der & (k=2;...n) auszunutzen, indem man zunächst nicht nach den £; selbst fragt, 
sondern einerseits nach den symmetrischen Grundfunktionen der &, mit ungeradem %, 
andererseits nach den symmetrischen Grundfunktionen der & mit geradem %& (&, und 
&n+1 bleiben dabei außer Betracht; sie sind durch die Forderung (1.2a) schon festgelegt). 

Wir beschränken uns im folgenden auf solche Fälle, in denen f rationalin ist, 
also darstellbar durch 


M 
SIR, 2... 00) 
u—=0 


= 


a; 
a Sl re Se 
u=0 


Gehören die Parameterwerte A, zu einer Lösung des Systems (1.2), so sei zur Abkürzung 


BT a 


Say Ar. A)eile e 
Sri dee 0 Re SR (1.4b) 
u=0 

und 
PS NE = Pla) rar Fa ine . (1.4c) 
u=0 


») Vgl.z.B.Tschebyscheff, Ges. Werke Bd. I, S. 271 (Petersburg 1899), Kirchberger, Über 
Tschebyscheffsche Annäherungsmethoden (Diss. Göttingen 1902) oder Cauer ‚ Theorie der linearen 
Wechselstromschaltungen, Bd. I, S. 548. 

?) Die Bedingungen (1.2a) bis (1.2d) besagen in Worten, daß für die speziellen Parameterwerte a, = Ay die 
Funktionen F(2) = f(x; A,);*** An) an (mindestens) n +1 aufeinanderfolgenden Stellen & des Approximations- 
intervalls / mit wechselndem Vorzeichen den Maximalbetrag |Z]| = Max |F| annehmen muß. 


(1.2c) folgt, die Existenz von fs vorausgesetzt, aus (1.2d). Für die rechnerische Behandlung ist esaber zweck- 
mäßig, (1.2c) eigens zu formulieren. 

®) Damit soll identisches Verschwinden des Koeffizienten von x 4 (und weiterer Koeffizienten) im Zähler 
oderim Nenner nicht ausgeschlossen sein; vorausgesetzt sei nur, daß nicht für alle (a) Rum = Su = 0. 
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gesetzt. Dann ist 


U 
| Erin 16.D bar tb. Boa re 

und nach (1.2b) 
UFLV/=0 an'den Stellen =’ & P wm 2 RD RERH). 


Von V (x) können wir annehmen, daß es im Intervall A < x < Bsein Vorzeichen nicht ändert, 
denn an einer in diesem Intervall gelegenen Nullstelle von V (2) müßte wegen (1.2d) auch U (x) 
verschwinden; die Polynome U und V hätten einen gemeinsamen Faktor und man könnte in 
(1.4d) oder, anders ausgedrückt, in (1.3) für den Fall der speziellen Parameterwerte a, = 4, 
durch diesen Faktor kürzen. Damit wäre aber die Aufgabe, die Bedingungen (1.2a) bis (1.2d) 
für die Funktion f aus (1.3) zu erfüllen, zurückgeführt auf eine einfachere Aufgabe, bei der die 
Polynome im Zähler und im Nenner von niedrigerem Grad als M sind ?). 

_ Aufgrund dieser Überlegungen können wir jetzt V (x) im Intervall Asxz<sB ohne 
Einschränkung der Allgemeinheit positiv voraussetzen. Aus (1.2d) und (1.4d) folgt dann nicht 


U 
nur — |Z] sy = |Z], sondern auch —|Z|-V— U<s0O<s|L|-V— U und damit z.B. 


auch einheitliches Vorzeichen von U—LV imganzen Intervall. Daher müssen diejenigen ö, im 
Innerndes Intervalls (4 ; B), an denen das Polynom U— LV verschwindet, (gerades k), 
notwendig doppelte Nullstellen sein. Entsprechendes gilt für ungerades k und U(s) + LV (2). 

Wir betrachten jetzt den Fall, daß » (also die Zahl der Parameter a,)eine ungerade 


; s 0 uk: 
Zahl ist. Dann sind nach (1.2a) je I im Innern des Intervalls (A ; B) gelegene Stellen 


g mit geradem bzw. mit ungeradem %k vorhanden und man erhält einerseits für ungerades 
Be 


n—$ n—3 
N): un a in; 
en 
andererseits für gerades Kke2a:..n-\ 
RT n=3 
I (@— &)=«' RE Fe N ei! 


2 

Das Polynom U + LV ist nach (1.4b) und (1.4c) höchstens vom Grad M. Nach (1.5) 
und den daran anschließenden Überlegungen enthält es die in (1.6a) zusammengefaßten 
Faktoren als doppelte Nullstellen und außerdem den Faktor & — £ı. Unter Berücksichtigung 
von (1.2a), wonach &, = 4, folgt daher 


n—1 n—3 


2 
U+LV/=(z2— 4) E u ut ng Frage 1) os 2 rn GR L>7E) 


2 


mit p= M—.n, wobei allerdings die Möglichkeit c,= 0 (und das Verschwinden weiterer Koefli- 
zienten) nicht von vornherein ausgeschlossen ist. 
Entsprechend gilt unter Berücksichtigung von (1.6b) 


n—1 n—3 2 
U-IV=(s— nl. run? re (dad, art...) . - (1.7b). 
= 


Durch Addition von (1.7a) und (1.7b) folgt 


n—1 2 n—l 2 
20-(6-4|2' au (War oreBi|e = ] (d,aP+...d,) (1.8a), 
2 


2 


durch Subtraktion 
n—1 \2 2 
2 2 
2 IV =(8— ale +... no) (Gar+...,)— el: +.. rn (d,@P+...d,) 
. : (1.8b). 

4) Gegeben sind in (1.3) die Funktionen A, und S,„ der Parameter a,. Unter der Voraussetzung, daß Zähler 
und Nenner einen gemeinsamen Faktor haben, kann man auch ohne Kenntnis der speziellen Systeme (4,), für die 
dies zutrifft, mit Hilfe des Euklidischen Algorithmus zunächst die Struktur dieses gemeinsamen Faktors, etw& 
2 —r(@ ;**-4y), ermitteln und dann Zähler und Nenner durch diesen Faktor dividieren (die formal zunächst noch 
auftretenden Reste sind dabei voraussetzungsgemäß für die betrachteten, wenn auch nicht numerisch bekannten 
Systeme a, gleich 0). 


BD CNTBRAT 
Bi 
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Die rechten Seiten dieser Identitäten in x sind nach der Bemerkung im Anschluß an 
(1.7a) formal (wenn die c, und die d, als Variable aufgefaßt werden) vom Grad M. Die linken 
Seiten sind nach (1.4b) und (1.4c) keinesfalls von höherem Grad. Wir erhalten daher durch 
Vergleich der Koeffizienten von & aus (1.8a) und (1.8b) ein System von genau 2- (M +1). 
Bestimmungsgleichungen für die n Unbekannten 4,, die nach (1.4b) und (1.4c) in U und F 
auftreten, für die n — 1 Unbekannten s, und t, sowie für die 2- (M — » + 1) Unbekannten.c, 
und d, und für die Unbekannte L. Dies sind insgesamt 2 M + 2 Unbekannte, d.h. ebenso- 
viele Unbekannte wie Gleichungen vorhanden sind. Die Anzahl 2 M + 2 dieser Gleichungen 
wird im allgemeinen größer sein als die Anzahl 2» der Gleichungen des Systems (1.2b) ; 
(1.26). Aber der Grad der einzelnen Gleichungen ist hier wenigstens unter der Voraus- 
setzung, daß die R,(a,) ganze rationale Funktionen von höchstens 3. Grad und die S,(@,) 
ganze rationale Funktionen von höchstens 2. Grad bedeuten, selbst höchstens gleich 3. Auf 
der rechten Seite von (1.8a) und (1.8b) kommen nämlich immer nur Ausdrücke der Gestalt 


Ode S% S, Co bzw. Di e t, t, de 


als Koeffizienten von x” vor (wobei die C,,. und D,,. bekannte Zahlenkoeffizienten sind 
und =, = 1 eingeführt ist). Dagegen waren die Gleichungen (1.2b) allein schon in den 
&, vom Grad M (wie sich ergibt, wenn man den Nenner der linken Seite als Faktor zu Z auf 
die rechte Seite nimmt), die Gleichungen (1.2c) in den &, mindestens vom Grad Men 
höchstens vom Grad 2M — 1°). 

Für das Weitere sei die Einschränkung angenommen, daß die R, und die S, ganze 
lineare Funktionen der a, sind. Dann gibt es außer den durch Koeffizientenvergleich aus 
(1.82); (1.8b) folgenden 2(M +1) Gleichungen 


Ra 52. u) = 2 Opre un He Fr Dusche > - = su. mar 
IS EIG Bar Did EI EIERN ED) 


noch (wenigstens)2 (M + 1) — n linear unabhängige Gleichungen zwischen den R, und den 8, 
(worunter auch Beziehungen zwischen den R, allein bzw. den S, allein sein können). 
Setzt man (1.9a) und (1.9b) in diese Gleichungen ein, so erhält man Beziehungen der 
Gestalt 
>= (Aue -L + Pxae) 550 + (Yrae L ie Ögie) 7) de =0...... (1.10), 


wobei die &; ß; y; ö bekannte Zahlenkoeffizienten sind. Damit hat man ein System von 
Gleichungen höchstens 4. Grades in den Unbekannten L; s,; tx; co; de. Das sind 


1+(n—1)+2(M— n+1)=2 M—n+2 


Unbekannte, für die ebensoviele Gleichungen zur Verfügung stehen. 
Falls die Anzahl n der Parameter nicht, wie bisher angenommen, ungerade, sondern 


gerade ist, gibt es nach (1.2b) im Inneren des Intervalls (A; B) 5 Argumentwerte &, mit 


geradem k, dagegen nur > — 1 Argumentwerte &, mit ungeradem k. Statt (1.7a) bzw. (1.7b) 


hätte man daher 


n_ 2 
DH We BlF 4 mot... 


2 
und 


2 


N 2 
Ua =[Ü+...0) art... 


Im übrigen bleibt der Gedankengang für die weitere Behandlung im wesentlichen derselbe 
wie bei ungeradem n, so daß sich eine genauere Durchführung hier erübrigt. 

Inwieweit das geschilderte Verfahren die numerische Behandlung der gestellten Aufgabe 
tatsächlich vereinfacht, hängt vor allem entscheidend von der Differenz M — n ab (je größer 
M bei festem n ist, um so mehr Gleichungen sind in (1.10) enthalten). Wesentlich für die 
praktische Brauchbarkeit der Gleichungen (1.10) ist, daß in vielen Fällen, wie auch aus den 
folgenden Beispielen ersichtlich ist, ein beträchtlicher Teil der Gleichungen, wenn nicht sogar 
alle, noch von niedrigerem als dem vierten Grad ist und daß sich beim Eliminieren einzelner 
Unbekannter noch weitere Vereinfachungen ergeben. 


5) Vgl. Fußnote 3. 


u 
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11. 
Wir betrachten die Aufgabe, 


2° +a,202 +0,20? +a,2? 


Max. le ae 
. +1 { 2 
möglichst klein zu machen ®) : 
Nach (1.1); (1.4a); (1.4b) und (1.4c) ist dann 
U=4%° +42 +74,° +4, —- ,2—1 
und a BRaN 
V == 4; % En 1 
(1.7a) und (1.7b) ergeben für das vorgelegte Beispiel die Identitäten 
U+l!V/ = I) +10 +s)0e..-. 2.2.20. (2.32) 
und 
U— IV =(c.—3) Ho +t)inde sohn. seien: 


Durch Koeffizientenvergleich unter Berücksichtigung von (2.2) folgt daraus, angefangen 
bei &°, zunächst 


9 EM EEE ER EL EN EDRELEL BEER SEEN NEF (2 4a) 
4,=c(2s,—1)=d(2t,—3), also 1 =sı +1 nach (24a) .. . . . (2.4b) 
L,=c +25 —-2s)=dl #2, —66) ale hei ts +} 

nach (2,42) (2:4b)am im. mild ei (2.4c) 

4,=c2 498 —- 93 —28)=dRh,hu 34H —6L), leo = —is1— 
nach (2.4a).; (2.4b) u.-(2.4c). . .... . (2.4d) 

83; t, und t, lassen sich also folgendermaßen durch s, ausdrücken: 

g=—3 1 — yaschis ieh de A (2.5). 


Durch Vergleich der linearen Glieder in (2.3a) und (2.3b) ergibt sich wegen (2.2) 


— 4,:(1—L)=c-.(&—2s, s,) 
In a Ve De ie ee era (2.62); 
— 4 (1 +L)=d.($ —6it, t5) 
Entsprechend folgt durch Vergleich der konstanten Glieder 
—(1—-D)=c-3 
und een Ars neg Einee Flair SEP TIIE EAEE (2.6b) 
—(1+D)=d-3% 
Wegen (2.4a) folgt aus (2.6) und (2.6b) schließlich 
t = L s ——_ 2 5189 82 
a Te (2.7) 
und daraus wegen (2.5) schließlich 
BASTI LE SO a 22a 
Von den Lösungen dieser Gleichung („= — 2 und ,=— !%) ist nur 3 = — 1 für 
die gegebene Aufgabe brauchbar. Damit erhält man unter Berücksichtigung von (2.5) 
a RE a TEE RR (2.8). 
Nach (2.7) folgt daraus 
A AO IT... - 22er (2,80) 
und daher aus (2.6b) 
BEN 093 294 10! ir Wesens eh. 


Dies ergibt zusammen mit (2.4a) bis (2.4d) 
A, = 0,0532 ; A,= — 05151714; 4,=1,92450; ,=—346410; L=—3 (2.9). 


6) Zur Anwendbarkeit des Gleichungssystems (1.2b); (1.2c) vgl. die Bemerkungen nach (2.10). 


a ne es a 
unter den Absolutstrichen stehenden i 
betrag, so müßte F— f in diesem Dat i 
Dies ist aber nicht möglich, weil nach (2.1) und (2.10) 


(ed astiheh EEE cal =: 
RE: Art) at) ER h 


weil also F- fi im Intervall A; 3) höchstens 4 Nullstellen haben Be Demi ist die Au 
(2.1) gelöst. i 
Speziell im Fall V =1 kann es sc kmähig sein, außer U-+ L auch noch die Ableitung 


U eB zu betrachten, da U’ dann sicher Bi Nr &,) als Faktor ee muß. Wir, 


dz 
betrachten dazu als Beispiel die Aufgabe, 


Max|* +0, 2,2 +, 82 us +9... ... ie 


0<r<s1 . 
möglichst klein zu machen. Mit 7E 
U= M=#+4, PAR +AR+ACtA: y_=1=.% 7212) 3 
läßt sich dann nach (1.7a) und (1.7b) setzen Re 
Url=e- ER ErLrsn er 
U-L=( 1) ® +h2 +6)” (0 — D. er Mana 21h 


(Der Koeffizient von x° muß wegen (2.12) auch rechts in (2.13a) und (2.13b) gleich 1 sein). 
Daraus folgt durch Koeffizientenvergleich, angefangen bei 2°, 


EN ER EN ER WE ln N 6) 
und daraus weiter 

A,=2:—-389 a ey De a ee ir 3 
A,=—-2,8— 23 = —24, —2# HAN — 204... 0.0. 0. (2146); 
L,=—41s +3 = 4, HÜH28- 2 14,6, — 2%... (21dc); 
A=-23u8 = — ul Altea Alıl,, aller lad) 

und schließlich , +L=0; 4, —L=d-t, also wegen (2.14) 
2 Ars 21 BET a me . (2.14e) 


(2.14a) bis (2.14d) ergeben 4 Gleichungen für die 4 Unbekannten s}; 83; tı; to. Wir 
nehmen dazu jetzt aber auch noch 


U=6(® +0 +8) HheH+tt)(e—b). . S1a.15, 


Unter Berücksichtigung (2.12) erhält man hieraus wieder durch en 
angefangen bei 29, zunächst 


Det, a re, 
und daraus weiter 
Ast la th - I —- sh]. ..2..02..0.(216a); 
4,=2 [66h — 1 — Sl, — Str)... vu Sr (2.16b); 
4=3 [se — (1 +) (set Hiesı)] » - > >.» (160); 
A067 I )sate 2 (2.16d). 


Damit stehen (nach Elimination der A,) zusammen mit (2. 14a) bis (2.14d) im ganzen 
sogar 8 Gleichungen für die Unbekannten s}; 83; {,; t, zur Verfügung, aus denen man nun die 
Gleichungen niedrigsten Grades bzw. einfachsten Aulbaues auswählen kann: 
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Aus (2.14a) und (2.16a), die beide in s, und t, linear sind, erhält man 


4, =33 36h +5, 476, —3 ....0.. 0.0. (2.17a) 
Ai a le ce Dat, A ti at re, 

Aus (2.14d) und (2.16d) folgt 
13 =- 3 +)sh= HH +t2Uw—Uls. - - - . (2.18) 


und daraus nach Division durch s, bzw. durch t, unter Benutzung von (2.17a) und (2.17b), 
wenn man jetzt statt s, bzw. t, einfach s bzw. i schreibt, einerseits 


ae — Wer at 3er 2er N 2.180) 
andererseits 


3 +10 st Hs — AR 16s5t +88 — 35 —4t=0 ... . (2.18b). 


. Das sind 2 Gleichungen 3. Grades in sund ti. Um sie zu lösen, kann man aus ihnen z.B. 
einerseits 2°, andrerseits s? eliminieren. Auf diese Weise bekommt man einerseits eine quadra- 
tische Gleichung für tin Abhängigkeit von s, andrerseits eine quadratische Gleichung für sin 
Abhängigkeit von t und kann dann den gesuchten Schnittpunkt der beiden zugehörigen 
Kurven bequem graphisch näherungsweise bestimmen. Dabei kommt von vornherein nur 
—2<s<0Q0und — 2<t<0in Betracht; denn die Extremstellen &,, deren negative Summe 
durch s bzw. t gegeben ist, müssen zwischen 0 und 1 liegen. Es ergibt sich 


SI = Sr — 1,27 ; t, SI er 0,77 sh rn Sa (2.19a) 
und damit aus (2.17a); (2.17b) 
s,= 0,33; 5.68 0.0682 8, 00.2) 3: ga Sera BR 


Mittels (2.13a) läßt sich daraus unter Benutzung vom (2.14) und (2.14e) die gesuchte 
Funktion U (= F) ermitteln. Es ergibt sich damit näherungsweise 


U(=F) = (a? — 1,27 x + 0,33)? («2 + 2,54) — 0,006. 


Daß diese Funktion (innerhalb der Grenzen der Rechengenauigkeit) die in (2.11) gestellte 
Aufgabe erfüllt, ergibt sich aus der Existenz von 6 Extremwerten gleichen Betrags und ab- 
wechselnden Vorzeichens im Intervall (0; 1) in ähnlicher Weise wie bei der Lösung (2.10) der 
vorausgegangenen Aufgabe (die Differenz zweier Funktionen der Gestalt 2 ta, Mt +...9 
kann höchstens viermal ihr Zeichen wechseln). 


Eingegangen am 6. Juli 1953. 


Die Ausbildung eines Elastizitäts-Störmoments beim 
schwingenden Pendel mit elastischer Schneide bzw. 
elastischer Unterlage 

Von H. Tesch in Bissendorf (Hannover) 


Verformie Oberflächenteile der elastischen Schneide oder Unterlage eines Pendels in der Nähe der Druck- 
flächenbegrenzung werden beim Schwingungsvorgang periodisch belastet und entlastel. Wegen der Trägheit 
des elastischen Mittels gegenüber Formänderungen entsteht durch Aufwalzung ein Blastizitäts-Störmoment, 
das die Schwingungszeit des Pendels verkürzt. 


Ifa pendulum ist oscillating, the deformed parts of ihe surface of ihe elastic knife edge or base next the 
plate of pressure are periodically loaded and unloaded. T'he rolling motion generates an elastic interference 
moment, that is due to the inertia of the elastic medium against deformation, and that has the effect of 
shoriening the period of oscillation of the pendulum. 


Des parties deformees de la surface du taillant &lastique ou dela base d’un pendule pres du plat de pression 
sont chargees et dechargees periodiquement pendant l’oscillation. En consequence de l’inertie du milieu 
&astique envers des deformations le mouvement roulant produit un moment de derangement de l’&lasticit£, 
lequel raccoureit la periode d’oscillation du pendule. 


edopMmupoBaHHzıe YACTH TeBEPXHOCTH YHPYTOTO HOMA HJIH ONTOPEI MAATHUKRA, HAXONAIIHECH 
BÖNHSH OTPAHHYeHHSI HAKHMHOH TOBEPXHOCTH, NOABePTaWTcH B IPONeCCE KOTeÖAaHHA TepHoNn- 
yeckoü Harpyske u pasrpyske. MHepıma yıpyroi cmpenst LO OTHONIeHMO AedopMmannli 
BEISBIBACT, BCIIEACTBUe PASBAJIBIOBAHHA, YIIPYyTut BO3MYINArOIMli MOMEHT, CORPAILSIOIIMÄ 
HPOANOIKHTENIBHOCTB KAYaHHÜ MAATHURA. 


Frühere Untersuchungen des Verfassers haben gezeigt, daß die Elastizität der Schneide oder 


der Unterlage eines Pendels eine Verkürzung der Schwingungszeit gegenüber derjenigen des 
starren Falles bewirkt. Dieser Sachverhalt mußte gedeutet werden als Auswirkung eines Elasti- 
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zitäts-Störmoments, das in bestimmter Weise von der Amplitude abhängig ist. Die vorliegende 
Arbeit handelt von der Ursache dieses Effekts. i 

Nach den hier wiedergegebenen Untersuchungen erweist sich die angegebene Störung in 
erster Linie als Folge der Trägheit des elastischen Mittels gegenüber Formänderungen. Das 
Material der elastischen Schneidenwalze bzw. der elastischen Unterlage erfährt in der Nähe der 
Druckflächenbegrenzung bestimmte Randverformungen, die im Rhythmus der Pendelschwin- 
gung periodisch belastet und entlastet werden. Da besonders in solchen Oberflächenzonen, die 
von der Druckfläche verhältnismäßig weit entfernt liegen, das elastische Mittel den kurzen Be- 
lastungseinwirkungen während einer Schwingung nicht augenblicklich folgen kann, kommt es 
beim Aufschwingen des Pendels vor der Walze bzw. in ihrem vorderen Teil zur Aufwalzung des 
bereits verformten Randmaterials. Die Folge ist die Ausbildung eines Elastizitäts-Störmoments. 
Bestimmte Annahmen über die Stärke des Normaldruckes in der Druckflächenbegrenzung und 
die näherungsweise Berechnung der Druckflächenvergrößerung beim Anschwingen des Pendels 
ergeben die Möglichkeit zur mathematischen Formulierung des Effekts. Ein Vergleich der Rechen- 
ergebnisse mit den Auswertungen derfrüheren Schwingungsversuche bestätigen im wesentlichen 
die Theorie. Statische Momentenversuche, die zur Klärung wichtiger Fragen beitragen, schließen 
die Arbeit ab. 


Einleitung 

Wird bei einem physischen Pendel die im Querschnitt kreis- oder keilförmige starre 
Schneide (Metall), oder die ebene starre Unterlage oder werden beide Teile durch gleichgeform- 
tes, elastisches Material (Kautschuk) ersetzt, so erfährt hierdurch die Schwingungszeit in je- 
dem Falle eine Verkürzung, die in bestimmter Weise von der Amplitude abhängig ist. 

Dieses Ergebnis hatten frühere, vom Verfasser auf Anregung vonM. Schulerin Göt- 
tingen durchgeführte Untersuchungen [1]. 

Quantitativ ergab sich der folgende Sachverhalt: 

Während die experimentell gewonnenen Schwingungszeit-Amplitudenfunktionen bei 


starren Verhältnissen, für keil und walzenförmige Schneiden hinreichend genau, durch die be- 


kannte Näherungsformel 
2 
Ta Tolı+%) ee 


wiedergegeben werden konnten, ergaben sich bei elastischen Verhältnissen in jedem Falle Funk- 
tionen von der Form 


2 
Da et nn +4:'9— B:p? ER (2). 
Das Verhältnis To, zu T(o) war gekennzeichnet durch die Beziehung 
To,< To. 

Bezeichnungen: 

p = Amplitude 

Tg) = Schwingungszeit des starren Systems bei der Amplitude 

Tg, = Schwingungszeit des elastischen Systems bei der Amplitude 

To) = Schwingungszeit des starren Systems für 9—>0 

To, = Schwingungszeit des elastischen Systems für 9— 0 


A, B = Positive Konstante. 


Ähnliche Abweichungen der Schwingungszeit-Amplitudenfunktionen vom theoretisch gefor- 
derten Verlauf der Gl. (1) sind häufig bei sehr genauen Pendelschweremessungen beobachtet 
worden. Sie haben den Anlaß zu den Elastizitäts-Untersuchungen gegeben. Bei der Klärung 
dieses Effekts führte die bezeichnete Arbeit des Verfassers zu dem Schluß, daß unter den ver- 
schiedenen Möglichkeiten zur Deutung der Abweichungen, diejenige der Ausbildung eines 
Elastizitäts-Störmoments M, die treffendste ist. Quantitativ ergab sich für dieses Moment bei 
kleinen Amplituden, in allen Fällen die Form: 


VASE Ar een 


wenn a die jeweilige Pendelauslenkung bedeutet und %,, k, positive Konstante sind. Aufgabe 
der vorliegenden Arbeit ist es, die nähere Ursache dieses Störmoments zu ermitteln. 


bar & 


ht - u La 
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2 Teil I. Theorie 
FallI. Das Pendel mit starrer Schneidenwalze auf ebener, 
elastischer Unterlage 


1.Die Verformungen der elastischen Unterlage beiruhendem 
Pendel. 


Bild 1 stellt das Querschnittbild der mit dem Pendelgewicht P belasteten starren Schnei- 
denwalze auf der verformten Unterlage dar. 


Gelten die Abkürzungen: 
P = Pendelgewicht, 
r = Radius der Schneidenwalze, 
! = Länge der Schneidenwalze, 
E = Elastizitätsmodul der Unterlage, 
m =Poissonsche Konstante, 


26°’—= Abstand AB = Drucklinie, 
so kann die letztere, bei Voraussetzung eines ebenen Spannungszustandes, nach Hert zdurch 


die Formel: 
il 
Bulıne) jur lets. in race) 


Ze l:n: BE 
berechnet werden. 
Bedeutet % den Abstand eines Punktes U der Berührungskurve (Bild 1) von der Mittel- 
linie MC und ist p, der Normaldruck im tiefsten Punkte der Walze (% = 0), so herrscht nach 


Hertzander Stelle Uein Normaldruck von der Größe 


Für die Drucklinienbegrenzungspunkte A, Bergibt sich hieraus miitau = +6, p =(. In Wirk- 
lichkeit bildet sich bei einem entsprechenden Belastungsversuch in der Nähe der Begrenzungs- 
punkte A, Beine Oberflächenverformung nach Bild 1, nicht aus. Durch den Belastungsvor- 
gang, mit dem mehr oder weniger ein Haften des elastischen Mittels an der starren Walze ver- 
bunden ist, werden die Oberflächenschichten gedehnt. Diese Schichten haben das Bestreben 


0, 0 


324 --- > 


Bild 1. Starre Schneidenwalze aufelast. Unterlage Bild 2. Starre Walze auf elast. Unterlage. 
Druckverhältnisse nach Hertz. ; Randverformungen 


sich zu verkürzen, was, bei entsprechenden Gegenspannungen mehr im Innern des elastischen 
Mittels, zu der in Bild 2 gezeichneten Oberflächenverformung führt. Wesentlich für die folgen- 
den Untersuchungen ist die Senkung der Punkte A’, B’. Sie können nicht, wie bei Bild 1 als 
druckfrei gelten. Gemäß Bild 2 mögen folgende Bezeichnungen festgelegt werden: 

9,6 — wirkliche Drucklinie = A’B’, Bogen A’ 0 B’ = Berührungskurve, Oberflächen- 
begrenzung O, 4’ C B’ O0, = Verformungskurve. 

Für die spätere Rechnung ist die Kenntnis des Normaldruckes in den Drucklinienbegren- 
zungspunkten 4’, B’ nötig. Erfahrungsgemäß ist die Hertz sche Formel für ö' nicht exakt 
richtig, ihre Anwendung ergibt aber gute Näherungswerte. Dies berechtigt zu der Annahme, 

27 


or Z. angew. Math. Mech. 
394 Tesch, Ausbild.eines Elastizitäts-Störmoments beimschwingenden Pendel 34.34 Nr. 10/11 Okt./Nov. 1954 
S BEE er ee 


daß für den Punkt A’ bzw. B’ der entsprechende Normaldruck p4 bzw. ps’ durch eine, der 
Gl. (5) entsprechende Formel 
PA pw = Vor — 3. ee fe et wie, aa ab aEe ee (6a) 


erfaßt werden kann, wobei unter ö* eine entsprechende fiktive Drucklinie zu verstehen ist. 
Druckversuche des Verf., nach denen anscheinend die wirkliche Drucklinie immer kleiner ist als 
die nach Gl. (4) berechnete, können in diesem Zusammenhang von Wichtigkeit sein. Bei einer 
Näherungsrechnung für p4- nach Gl. (6a) wird es genügen, den Normaldruck 7, aus der Forde- 
rung: Y 

: Pendelgewicht — Summierung aller Normaldrucke der Berührungsfläche in der Form 


"us BER FET, d 
P = 2 [1:72,28 |YorZaau=t:n[ Vila) +0%-aresin 
ö ö 


zu ermitteln. 
Für p, folgt daraus: 


7 
Pr ae : 
ı[-yı-(%) + 0% -aresin S,| 


Mit diesem Wert geht Gl. (6a) über in die Form 


2 
> Vila) 
AZ nn . . . u ET (6b). 


2 
10% — (5) + ö* -are sin =) 


ö . 
Genauere Untersuchungen dieser Gleichung zeigen, daß für ein Intervall 0,8 << lim 


Nenner der Formel (6b) das Verhältnis A gleich 1 gesetzt werden kann. Das Rechenergebnis 


wird hierdurch nicht wesentlich beeinflußt. 
Bei dieser Einschränkung für ö/ö* folgt aus (6b) für 94, bzw. pp’: 


aB / ö\? 
a ve De (6) 


2.Der Anrollvorgang, Ausbildung eines Elastizitäts-Rück- 
führmoments. 


Das in Bild 2 gezeichnete Querschnittbild wird bei dem Belastungsvorgang nicht mo- 
mentan in endgültiger Form ausgebildet. Wegen der Erscheinung der elastischen Nachwir- 
kung ist die Verformung des elastischen Mittels erst nach einigen Minuten Standzeit praktisch 
abgeschlossen. Wird die Walze um einen kleinen Winkel a aus der Ruhestellung ausgelenkt, so 
ist mit einer Rollbewegung des Walzenkörpers zu rechnen. Dieser Vorgang hat, wegen der Träg- 
heit des elastischen Mittels gegenüber Formänderungen, eine Stauung des bereits elastisch ver- 
formten Materials unmittelbar vor der Walze zur Folge, derart, daß eine Verlängerung des vor- 
deren Drucklinienteils erfolgt. Beim Verweilen der Walze in der ausgelenkten Stellung ist eine 
Verminderung dieses Effekts zu erwarten. Noch in einer anderen Weise kann sich der Anroll- 
vorgang auf die Gestalt der Verformungskurve vor der Walze auswirken. Beim Aufsetzen der- 
selben auf die elastische Unterlage werden die Oberflächenschichten je nach Stärke der Be- 
lastung gedehnt. Wegen der Haftung des elastischen Mittels am Walzenkörper sind die Deh- 
nungen innerhalb der Berührungskurve als umso stärker anzunehmen, je näher sich das be- 
trachtete Oberflächenelement an den Begrenzungspunkten 4’ B’ befindet. Beim Anrollen 
wechseln die verschieden gedehnten Oberflächenelemente ihre Lage innerhalb der Berührungs- 
kurve und verstärken den oben angegebenen Störeffekt durch Verlängerung des vorderen 
Drucklinienteiles. Beim Verweilen der Walze in der ausgelenkten Stellung ist eine Verminde- 
rung des zuletzt angegebenen Störanteiles, wenn vom Ausgleich durch Kriechen des elastischen 
Mittels abgesehen wird, nicht zu erwarten. 


Durch die geschilderten Verformungsverhältnisse wird ein Elastizitätsmoment ausge- 
bildet, das der Anrollbewegung entgegengerichtet ist, und das sich beim schwingenden Pendel 
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als zusätzliches Rückführmoment auswirkt. Die Verformungsverhältnisse hinter der Walze: 
lassen eine entsprechende Beeinflussung des Rollvorganges in gleicher Stärke nicht erwarten. 

Für die mathematische Formulierung des Störeffekts wird die jeweilige Verlängerung des 
vorderen Drucklinienteiles mit e bezeichnet. Während ö die halbe Drucklinienlänge bei ruhen- 
der Walze, also eine Konstante bedeutet, muß e als eine Funktion der Auslenkung a angesehen 
werden. Zu einem Punkt der vorderen Berührungskurve, der von der Vertikalen M’C’ den 
Abstand (ö+ e) hat (Bild 3), möge der Normaldruck p gehören. Es ist naheliegend, den letzteren 
näherungsweise in der Form 


p=pa—K:e 6: sole Puastornehälttaelt sta 


ansetzen, wobei 94 sich aus Gl. (6) ergibt und K eine positive Konstante bedeutet. Der Mo- 
mentenanteil, den ein Flächenelement } - de (| = Walzenlänge) mit dem Normaldruck 7 liefert, 
ist näherungsweise 


M|= 9:8 el ulcidei on en 7 Le Ra 
Für das volle Elastizitätsmoment der Aufwalzung folgt daraus: 


MI = fam=1-foHll)-@e—K-l)-de 2.2... @b) 


Bei Durchführung der Integration werden höhere als quadratische Glieder in e vernachlässigt. 
GI. (8b) geht über in die Form: 


IMl=1:|6- pa: lel+(pe 8: MS ee ee 


In dieser Gleichung stellt e = e(&) die Drucklinienverlängerung für die jeweilig Pendelaus- 
lenkung dar. Sie kann, wiespäter erörtert wird, für Fall I, Ilund IIInäherungsweise in der Form 
lel = A: al — B" o® 

angesetzt werden. 
‚Dabei bedeuten A, Bpositive Konstante und das Größenverhältnis beider Werte ist durch 
die Beziehung B>A gekennzeichnet. 

Nach dieser Feststellung beeinflußt der Ausdruck in e? der Gl. (8c) den in «linearen Anteil 
des Störmoments nicht, doch ist zu untersuchen, ob mit einer Einwirkung auf dasin a quadrati- 
sche Momentenglied gerechnet werden muß. 4 

Der Faktor K wurde durch die Beziehung (7) p = pa — K: |e| eingeführt. (p4 — K 6) 
bedeutet danach den Normaldruck für |e| = 6. 

Unter der Annahme, daß (pa — K - ö)klein ist gegen p4- selbst, und bei Berücksichtigung 
der Beziehung B>A kann das in e quadratische Glied der Momentenformel gegen das lineare 
vernachlässigt werden. Die Momentengleichung (8c) nimmt somit die vereinfachte Form an: 


IM|=1-ö.40:plel 
und nach Einsetzen des Wertes p4- aus Gl. (6): 


an = Ay)  , 


Die Auswertung dieser Formel setzt u.a. die Kenntnis der 
Funktion e = e(«) voraus. Zu einer Näherungsrechnung für 
sie führen die nachstehenden Betrachtungen: 

Der Bewegungsablauf wird als Rollvorgang der Walze 
auf einer, zur unverformten Unterlagsoberfläche parallelen 
Ebene angesehen. Wird die Walze um einen kleinen Winkel 
in Pfeilrichtung ausgelenkt (Bild3), so erfährt der ursprüng- 
liche Drucklinienteil A’C’ die Verlängerunge etwa bis A”. 
Der neue Begrenzungspunkt A’ beschreibt während des 
Rollens in einem raumfesten Koordinatensystem eine Pid#. Fre auf elaal.. Unferianp.ve 
Kurve, die im folgenden als ‚‚Spurkurve‘‘ bezeichnet werden 
soll. Da praktisch mit einfachen Mitteln eine Krümmung dieser Spur nicht zu beobachten ist, 
kann sie als sehr flach, also näherungsweise als Gerade angesehen werden. Über ihre Lage in 
dem raumfesten System ist zu sagen, daß sie jedenfalls durch den Drucklinienbegrenzungs- 
punkt A’ bei ruhendem Pendel verlaufen muß. Die Funktion e — € (&) läßt sich angeben, wenn 
der jeweilige Schnittpunkt der Spurgeraden mit dem Rollkreis bekannt ist. Dieser ist daher, 
bei Annahme eines Wertes für die Neigung m der Spurgeraden, zu berechnen. Die Näherungs- 

27% 
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rechnung, die wegen ihrer Weitläufigkeit nicht wiedergegeben werden soll, liefert, bei Be- 
schränkung auf quadratische Glieder in a für e = e(a) die Gleichung: 
str bare hal ne eo Er ee 
Yet Lmb 2/02 — be F me. bi) 
wobei b den Abstand bedeutet, den (bei Ruhestellung des Pendels) der Walzenmittelpunkt von 
dem Schnittpunkt der Spurgeraden mit der durch M verlaufenden Senkrechten hat. Nach 
Einsetzen des Ausdruckes (9a) in die abgeleitete Momentenformel (8) ergibt sich für das so be- 
rechnete Elastizitätsmoment eine formale Übereinstimmung mit dem in der Einleitung ange- 
gebenen Störmoment Gl. (3). Spezielle quantitative Angaben lassen sich noch nicht machen, 
unbekannt sind die Größen m bzw. b und ö/ö*. 

Vereinfachte Annahme über m. 

Nach Gl.(9a) zeichnen sich folgende zwei Fälle der Geradenneigung, bei Berechnung der 
Funktion e = e(a) als besonders einfach aus: 


je| = 


a ac ge ie) Mar re s 
mist sehr klein, die Spurgerade läuft nahezu horizontal. 
MEDIEN U Een en se 
Die Spurgerade geht genau oder annähernd durch den tiefsten Punkt der ruhenden Walze. In 
Gl. (9a) kann das Glied (r2 = 5) gegen m? - b2 vernachlässigt werden. Mit br gilt m =. N 
In diesem Falle geht Gl. (9a) über in die Form: 
ler kg er al. al, uch ai 


Eine Lage der Spurgeraden in dieser Weise läßt sich praktisch an einem sehr groben Modell 
beobachten. Die Art wie unter diesen Umständen die Neigung m von der Belastung, dem Ra- 
dius und der Länge der Walze abhängig wird, ist anschaulich annehmbar. Den weiteren Aus- 
führungen sollen die Verhältnisse dieses Falles zugrunde gelegt werden. Die Anwendung der 
Endformel für M, auf die angegebenen Versuchsergebnisse muß zeigen, wie weit die gemachten 
Annahmen zu Recht bestehen. 

Nach Einsetzen der Formel (9) in die Momentengleichung (8) ergibt sich für das Elastizi- 
täts-Störmoment beim Anrollen einer starren Kreiswalze auf ebener, elastischer Unterlage die 


Näherungsformel 
2 Datsscht H 3 | r 
M,| = IR Hua fan gebstenga Sshreie 
1. 7 „/! n ie am “ wu 


Diese Gleichung stimmt nicht nur formal, sondern auch zahlenmäßig mit dem in der Einleitung 
angegebenen Störmoment bei Schwingungsversuchen überein, wenn, bei Annahme eines be- 
stimmten Wertes für ö/ö*, die angegebenen Konstanten k,, k,sich darstellen lassen in der Form: 


Pe en ö\2 2P-12 2 ö 
ki = B3 a! (5) 5 Ka. ee, yı a mit oe . (10a). 


Der Quotient ö/ö* wird für verschiedene Belastungsverhältnisse nicht als konstant gelten 
können. Wegen der guten Näherungswerte der Hertzschen Drucklinienformel bei ent- 
sprechender Anwendung ist aber eine starke Abweichung der fikitiven Größe ö* von der wirk- 
lichen Drucklinie ö nicht zu erwarten. Im folgenden soll die entwickelte Momentenformel (10), 
durch Anwendung auf die in der Einleitung angegebenen Schwingungsversuche, nachgeprüft 


a: s 6) 
werden. Dabei wird probeweise gesetzt - 0,9. In der nachstehenden Tabelle ist das Ver- 


suchs- und Rechenergebnis zusammengestellt. In den ersten Spalten sind die Pendeldaten, da- 

neben die Konstanten %,, k, aus den Schwingungsversuchen eingetragen. Die eingeklammerten 

Werte wurden nach Formel (10a) berechnet. In der letzten Spalte sind die Quotienten kıfk 

der Versuche, darunter eingeklammert, die Quotienten (k,/k,) der Rechengrößen verzeichnet. 
Sonstige Daten. 


Elastizitätsmodul der elastischen Unterlage. 
Versuch ],, I E = 0,6 Kilopond/mm?2 
Versuch I],, I E=0,9 
Versuch IIJ,, In E=0)9 


LE} Le} 


A >] 
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ur 
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Außer den Versuchen bei endlichen Schneidenradien sind zum Vergleich unter III,, III, der 
Tabelle zwei Beobachtungen für den Fall ,‚Starre keilförmige Schneide“ (r = 0) auf ebener, 
elastischer Unterlage‘‘ wiedergegeben. In diesem Falle ist als Ursache des Störungsmoments ein 
rein-elastischer'Wangeneffekt, daneben, wie vorher, die Trägheit des elastischen Mittelsgegenüber 
Formänderungen, anzusehen. Ähnliche Verhältnisse bezüglich der Momentenbildung, wie sie 
für die spitze Schneide gelten, sind bei sehr kleinen Schneidenradien zu erwarten. 


Tabelle Il. Momentenkonstanten für den Fall „Starre Schneidenwalze auf 
ebener, elastischer Unterlage“. Klammerwerte berechnet 


Bezeichnung ae. Pendel- epalbe Gesamt- Konstanten des Störmoments 
des Wal, gewicht es Walzen- in Kilopond. Meter 
Versuchs alze inie länge 
r er 0) I kı k, kjlkz 
Ia 1 3,36 0,18 2,8 6,4 - 10-3 3,2 - 10-2 23.101 
cm Kilopond cm cm (10 - 10-3) &,5- 102) (1,8 - 10-1) 
1 1,8 0,15 2,8 4,3 10-3 2,4102 1,8 - 10-1 
cm Kilopond cm cm (5,4 10-3) (3,6 - 10-2) (1,5 - 10-2) 
0,6 4,74 0,1 2 12,7 - 10-3 7,7 - 10-2 1,7 - 10-1 
em Kilopond cm (7,5 - 10-8) (4,5 - 10-2) (1,7 - 10-1) 
0,6 2,13 0,07 3,5- 10-3 3.2.1022 1,1 - 10-1 
cm Kilopond cm (3,8 - 10-3) (3,3 - 10-2) (1,2 - 10-3) 
0 3,96 0,9 - 10-3 0,5 - 10? 1,8 - 10-1 
Spitze Kilopond 
. [Schneide 
III» 0 1,9 0,3 - 10-3 0,2. 102 1,5. 10-1 
Kilopond r 


Aus der Tabelle ergibt sich zum Teil eine recht gute Übereinstimmung der nach Formel 
(10a) berechneten Werte %,,k, mit den Versuchsergebnissen. Es kann daraus gefolgert werden, 
daß die angeführte Theorie den Elastizitätseffekt für Fall Tim wesentlichen erfaßt. Damit er- 
weist sich auch die Annahme über die Lage der Spurgeraden, über die Druckverhältnisse in der 
Nähe der Drucklinienbegrenzungspunkte und über das Verhältnis der wirklichen Drucklinien 
zur fiktiven, näherungsweise als zutreffend. 


Eine weitere Diskussion dieses Falles erfolgt im Teil Il dieser Arbeit. 


FallI. Das Pendel mitelastischer Schneidenwalze aufebener, 
starrer’Unterlage 


1. Druckverhältnisse beiruhendem Pendel. 


Bild4 zeigt das dem in Bild1 entsprechende Querschnittbild der belasteten elastischen 
Walze. Die für FallI festgelegten Bezeichnungen gelten hier sinngemäß. Die Drucklinie 
AB = 26’ kann nach Gl. (4), die Normaldruckverteilung über AB nach Gl. (5) berechnet wer- 
den. Für die Begrenzungspunkte A, Bgilt nach Hertz,p =. 


Abweıchnne des- Querschnittbildes von der in Bild 4 oe- 
zeichnete Form. 


Bei entsprechenden Druckversuchen wird ein Querschnittbild nach der in Bild 5 ge- 
zeichneten Art beobachtet. Der Belastungsvorgang hat zur Folge, daß auch seitlich der eigent- 
lichen Druckzone noch merkliche Formänderungen der Walze erfolgen. Wesentlich für die 
nachstehenden Untersuchungen sind die Ausbuchtungen außerhalb der Drucklinienbegren- 
zungspunkte A’ und B’. 


Die elastische Walze wird entweder von einem entsprechenden Lager gehalten (Bild 5) 
oder sie stellt den unteren Teil eines, sonst keilförmigen Schneidenkörpers dar. In beiden 
Fällen entsteht eine gewisse Symmetrie der Verformung in bezug auf die Winkelhalbierende 
des Winkels 4’ M D bzw. B’ ME. Die Oberflächenverformung läßt sich im Querschnittbild 
durch einen Kreisbogen annähern, dessen Mittelpunkt auf der bezeichneten Winkelhalbierenden 
liegt (M’), und der durch den Punkt 4’ bzw. B’ verläuft. Nach diesen Feststellungen muß im 
Begrenzungspunkt A’ bzw. B’ ein endlicher Normaldruck p4’ bzw. pr’ vorhanden sein. 
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Da im vorliegenden Falle, wie bei Fall I die Berechnung der Drucklinie nach Gl. Hanz 
Näherungswerte ergibt, ist auch hier mit einer halbkreisförmigen Druckverteilung nach Gl. 0) 
zu rechnen. Es ist daher statthaft, p4- bzw. pp’ durch das Verhältnis der wirklichen Druck- 
linie ö zu einer fiktiven ö’ festzulegen (Gl. 6). 


Bild 4. Elastische Walze auf starrer Unterlage, Bild 5. Elastische Walze auf starrer Unterlage, 
ohne Rand verformungen mit Randverformungen. Ruhestellung 


2.Die Verformungsverhältnisse beim Anrollen der Walze 
Ausbildungeines Elastizitäts-Rückführmoments. 


Gefragt wird nach der Änderung des in Bild 5 gezeichneten Querschnittbildes bei Drehung 
der Walze um einen kleinen Winkel «. Der Bewegungsablauf soll als Rollvorgang angesehen 
werden, wirksamer Rollradius sei r (Bild 5). Beim Anrollen des Walzenkörpers erfährt das 
elastische Mitte] eine Änderung der Beanspruchung. Diese kann von solcher Art und Größe 
sein, daß das ursprüngliche Querschnittbild lediglich eine Translationsbewegung erfährt, in 
seiner Form jedoch erhalten bleibt. Dem wirken besonders zwei Umstände entgegen, einmal die 
Trägheit des elastischen Mittels gegenüber Formänderungen, dann die obere Begrenzung der . 
Schneidenwalze durch eine Halterung oder die Erweiterung zur Keilform. Beide Einwir- 
kungen haben eine Stauung des bereits elastisch verformten Mittels in der vorderen Aus- 
buchtung zur Folge. Beider Bewegung kommt außerhalb der ursprünglichen Drucklinie, im 
vorderen Teil der Walze weiteres Material mit der Unterlage in Berührung, was zur Ausbildung 
eines Elastizitäts-Störmoments führt. Im hinteren Teil der Walze können, wie sich zeigt, Stö- 
rungen in gleicher Stärke nicht entstehen. 


Wird die Drucklinienverlängerung, wie früher mit e bezeichnet, und der Normaldruck » 
im Gebiete der Aufwalzung auch näherungsweise in der Form p = p4’ — K - |e| angesetzt, so 
liegen für die mathematische Formulierung des Störmoments die gleichen Verhältnisse vor wie 
bei Fall I. 


Die Berechnung führt daher ebenfalls auf die abgeleitete Momentengleichung (8). 


Die nächste Aufgabe besteht in der Ermittlung der Funktion e = e(a). Hierzu ist folgen- 
des zu beachten: 


Bei einer Drehbewegung der Walze in Pfeilrichtung wandert der linke Endpunkt 4” der 
verlängerten Drucklinie auf der horizontalen Unterlage. Dabei beschreibt er in einem Koordi- 
natensystem, das mit den unverformten Teilen der Walze fest verbunden ist, eine bestimmte 
„Spurkurve‘‘. Wegen der ursprünglich kreisbogenartigen Randverformung der Walze wird 
auch die Spurkurve näherungsweise als Kreisbogen, der zugehörige Mittelpunkt ungefähr auf 
der Winkelhalbierenden MD’, angenommen werden können. Wird das walzenfeste Koordi- 
natensystem im Walzenmittelpunkt so fixiert, daß die eine Achse bei ruhendem Pendel senk- 
recht, die andere Achse parallel zur Oberfläche der Unterlage verläuft, so kann der Spurkreis- 
mittelpunkt in diesem System durch die Koordinaten x, (in der Horizontalen) und y, (in der 
Vertikalen) festgelegt werden. Die Drucklinienverlängerung eist durch Berechnung des Schnitt- 
punktes der horizontalen Unterlage mit dem, mit der Walze bewegten Spurkreis, zu ermitteln. 
Auf die Wiedergabe der Rechnung soll verzichtet werden. Als Ergebnis derselben folgt bei 


unwesentlichen Vernachlässigungen und bei Beschränkung auf lineare und quadratische Glieder 
in afür e = e(a) die Funktion: 


a Due 
. lvo u. 2), 


10) 


el 


ii 


pr 
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Wird dieser Wert in die obige Momentengleichung eingesetzt, so ergibt sich für das Elastizitäts- 
Störmoment im Falle einer elastischen Schneidenwalze auf ebener, starrer Unterlage die 
Gleichung: 


e2PR1S o.\2 "© i Rs Pa 
Me rt] 02 
Dieser Ausdruck stimmt formal mit dem in der Einleitungangegebenen Störmoment (Gl. 3)über- 
Er Um die Gleichung quantitativ anwenden zu können, ist eine Annahme nötig über 6/6*, =, 
und %9- 

Nach der Festlegung des Spurkreismittelpunktes im Walzenkreis ist y, von x, abhängig. 
Beide Werte sind von gleicher Größenordnung, wenn die halbe Drucklinie klein ist gegen den 
Schneidenradius. Da &, mittelbar die Stärke der Aufwalzung bestimmt, ist es naheliegend, diese 
Größe als abhängig von der Drucklinie anzusehen. Für #, = ö/2 geht Gl. (12) annähernd über 
in die abgeleitete Momentenformel des Falles I. Wie die Versuche zeigen, sind für Fall Iund II 
die ermittelten Momentenkonstanten wirklich von gleicher Größenordnung. Der obige Sachver- 
halt kann also bedeuten, daß mit x, ungefähr gleich ö/2 die Verhältnisse der Aufwalzung und 
Momentenbildung näherungsweise getroffen werden. 


In Tabelle II sind die Konstanten k,, k, der Schwingungsversuche des Falles II wieder- 
gegeben. Die eingeklammerten Werte stellen Rechengrößen nach den Formeln 


TER, ö\2 a 
k we Ei | ee 
ide ! () Ba n). (128) 
(1 427 win . 


autark a u Te a a 
aha ya) ar 4 | 
dar. 


Bei Anwendung dieser Gleichungen erweist sich die Größe y, gegenüber dem Wert x, als 
von geringem Gewicht. Es besagt dies, daß die Annahme des Spurkreismittelpunktes auf der 
Winkelhalbierenden M’M keine besondere Einschränkung des Problems bedeutet. 

Bei der Berechnung der Konstanten wurde in allen Fällen von der Annahme x, = 0,65 ö 
ausgegangen. Durch diese Beziehung, die Größe ö selbst und den Radius ist nach GI. (12) die 


Tabelle II. Die Momentenkonstanten k,, k für den Fall: „Elastische Schnei- 
denwalze aufstarrer,ebener Unterlage“. Falll 


: Radius Pendel- Halbe Wälzen- M ven ken 
Bezeich- ende zen omentenkonstanten 
nase on gewicht ar länge in. Kilopond - Meter 
Versuchs 
r P ö l kı | k, kylla 
Ia 9 Wr; NEE 2,8 5-103 | 12-10 0,4 - 10-1 
cm Kilopond cm cm (3,8 - 10-3) (5,7102) (1- 10-1) 
mit  — 0,98 
In en.) "Ale 02 2,8 3.10. |Wg.10- 1. 10-1 
cm Kilopond cm cm (8,22110-3) (4,6 - 102) 05731022) 
mit; = 0,98 
rt ig 4,10 0,23 2 9,6 - 10-3 9.10* 1,1. 10-1 
cm Kilopond cm cm (9,6 - 10-3) (7,5 - 10-2) (1,3 - 10-2) 
Berechnet 
mit zw = 0,95 
JasH, 0,6 2,13 0,15 2 07 75.103 | 7,86. 10% 1. 10-1 
cm Kilopond cm cm (4,3 » 10°) (5,1 10-2) (0,8 - 10-1) 
Berechnet Ö 
mit Ur 0,95 
Illa 0 3,56 2 |. 38-10 0,5 - 10-2 8. 10-1 
Spitze Kilopond cm 
Schneide 
IIIv 0 1,96 2 3,4: 1028 0,3 - 10-2 12 - 10-1 
| Kilopond | cm | | 


Elastizitätsmodul des Schneidenmaterials: 
Versuch ah, b>E=03 Se 
Ian, InT>E=0,8 s 
Illa ‚, Ip > E = 0,8 >5 = 
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Charakteristik der Momentenfunktion bestimmt. In dem Wert k,/k, äußern sich diese Verhält- 
nisse zahlenmäßig. Wie die Gegenüberstellung der Versuchswerte kı/k, und der a 
(k,/k,) zeigt, stimmen sie in der Größenordnung überein. Versuch I, fällt offenbar heraus. Für 
die Größe der Momentenkonstanten selbst ist theoretisch das Verhältnis der wirklichen Druck- 
linie zu einer fiktiven, mit entscheidend. 6/6* wurde bei Versuch I zu 0,98, bei Versuch I zu » 
0,95 angenommen. Proportionalität der Konstanten k,, k, mit dem Pendelgewicht besteht bei 
diesen Versuchen offenbar nicht. Das kann Abhängigkeit des Quotienten ö/6* von der Belastung 
bedeuten, und zwar im Sinne einer Verkleinerung dieses Wertes bei Vergrößerung des Pendel- 
gewichts. Da 6/6* durch die Randverformung in der Nähe der Drucklinienbegrenzung bedingt 
ist, wäre der obige Sachverhalt anschaulich verständlich. Allein wegen der Kreisform des 
Schneidenkörpers sind bei geringer Belastung verhältnismäßig starke Deformationen in den 
Drucklinienbegrenzungspunkten zu erwarten. 

In den letzten Zeilen der Tabelle II sind die Elastizitätsmomente für den Fall „Spitze, 
d.h. keilförmige elastische Schneide auf ebener, starrer Unterlage“ verzeichnet. Da mit dem 
Belastungsvorgang auch in diesem Falle kreisbogenartige Verformungen in der Nähe der Druck- 
linienbegrenzung verbunden sind, liegt ein Sonderfall des besprochenen Falles 11 vor. Weitere 
Feststellungen über die Ausbildung eines Störmoments beim Pendel mit elastischer Schnei- 
denwalze werden in dem Abschnitt ‚Statische Momentenversuche‘ wiedergegeben. 


Fall III. Das Pendelmitelastischer Schneidenwalze auf ebener, 
elastischer Unterlage 


Nach den Untersuchungen des Falles I und II ergibt sich die Behandlung des Falles III 
zwangsläufig. Die Deformation der Walze erfolgt hier im Sinne des Falles II, die Verformung 
der Unterlage nach den Ausführungen des Falles I. Die frühere Feststellung bezüglich des Nor- 
maldruckes in der Druckflächenbegrenzung sowie die Ansicht über die Verlängerung des vor- 
deren Drucklinienteiles beim Anrollen der Walze ist hier zu übernehmen. Daraus folgt, daß auch 
in diesem Falle die Ableitung des Störmoments auf die allgemeine Momentengleichung des 
Falles I Gl. (8) führt. 

Die Funktion e = e(a) kann durch Aufsuchen des Schnittpunktes einer raumfesten Ge- 
raden (Spurgerade im Sinne des Falles I)mit einem, zum Walzenmittelpunkt exzentrischen, 
walzenfesten Kreise (Spurkreis im Sinne des Falles II) berechnet werden. Nach Einsetzen der 
so erhaltenen Drucklinienverlängerung in die Gl. (8) ergibt sich formal ein Störmoment wie es 
im Falle I und II berechnet und bei den Schwingungsversuchen des Falles III beobachtet 
wurde. Zur Anwendung eignet sich das Rechenergebnis wegen einiger, schwer feststellbarer 
Verformungsgrößen nicht. 

Ein günstigeres Ergebnis ist zu erhalten, wenn statt des angegebenen Spurkreises eine, 
zum Walzenmittelpunkt symmetrische, ‚‚Spurellipse‘‘ in Rechnung gesetzt wird. In diesem Falle 
ergibt sich für das Störmoment eine Gleichung, die nahezu mit der Momentenformel des Falles I 
übereinstimmt. Die Momentenkonstanten der Schwingungsversuche ergaben sich ebenfalls 
nahezu gleich den entsprechenden Werten des Falles I. Aus diesem Sachverhalt kann gefolgert 
werden, daß bei einem Pendel mit elastischer Schneidenwalze und elastischer Unterlage die 
Momentenbildung im wesentlichen durch die Verformungen der Unterlage erfolgt. 


Teil II. Statische Momentenversuche 


Die in Tabelle I und II wiedergegebenen Momentenkonstanten stellen die Glieder-eines 
Störmoments dar, das sich aus Messungen der Schwingungszeit des schwingenden Pendels er- 
gab. Der in Teil Ientwickelten Theorie der Momentenbildung lagen die Verhältnisse einer, aus 
der Ruhelage ausrollenden Walze zugrunde. Daß die hierbei gewonnenen Vorstellungen auf den 
Schwingungsvorgang des Pendels übertragen werden können, hat nach den folgenden Ver- 
suchen seinen Grund in der verhältnismäßig großen Trägheit des, von der Druckfläche ent- 
fernteren elastischen Mittels gegenüber Formänderungen. Gegen diese Verzugszeit ist die 
Dauer der im Rhythmus der Pendelschwingung periodisch wechselnden Belastung und Ent- 
lastung des Randmaterials gering. Nur ein Teil der Verformungsfläche steht beim schwingen- 
den Pendel dauernd unter nahezu gleichem, verhältnismäßig hohem Druck. Nach diesen Vor- 
stellungen muß das Elastizitätsmoment mehr oder weniger stark auch bei statischer Momenten- 
belastung zu beobachten sein, eine Vermutung, die durch die im folgenden beschriebenen Mo- 
mentenversuche bestätigt wird. 

Die Versuchsanordnung. 


An dem in Bild 6 im Querschnitt gezeichneten Pendel lassen sich die Versuchskörper in 
Form von festen oder elastischen Hohlwalzen konzentrisch zur Pendelachse anbringen. Sollen 
die Versuche mit vollen Schneidenwalzen vorgenommen werden, so dienen als konzentrische 
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Halterungen zwei, auf beiden Seiten der Achse aufgeschobene Rohre mit entsprechenden Aus- 
nehmungen an ihren Enden (Bild 6). Ein Faden, dessen eines Ende vermöge des Stiftes San der 
Pendelachse befestigt ist, wird in der gezeichneten Weise oben über die Pendelachse geführt. Er 
trägt an seinem freien Ende die Gewichtsschale @. Bei Belastung der Schale erfährt das Pendel 
eine Auslenkung a, die mit Hilfe eines Lichtzeigers an der Skala Sk abgelesen werden kann 
(Spiegel Sp am Pendelkörper). AB stellt eine hebelartig gelagerte Unterlage dar (Drehpunkt 
4), mit deren Hilfe eine allmähliche und 
ruhige Senkung der belasteten Schale @ 
ermöglicht wird. 


Ist das Hebelende B soweit gesenkt, 
daß die Belastungsschale frei schwebt, so 
herrscht Gleichgewicht zwischen dem 
„Schalenmoment‘‘ und dem Rückführ- 
moment des Pendels. Bedeutet R den Ra- 
dius der Pendelachse, @, das Gewicht der 
belasteten Schale, P das Pendelgewicht 
und a den Abstand des Pendelschwerpunkts 
vom Walzenmittelpunkt M, so gilt in die- 
sem Falle bei starren Verhältnissen 


GER = Pe,.a=sin.a: 


Wird die Schale oder die Unterlage, oder 
werden beide Teile durch gleichgeformtes 
elastisches Material ersetzt, so ist zur Aus- 
lenkung des Pendels auf den gleichen 
Winkel & ein größeres Schalenmoment 
@,: Rerforderlich. In diesem Falle lautet 
die Gleichgewichtsbedingung: 


G,:R=M,+P-a:sino, 


Bild 6. Versuchsanordnung für Momenten-Versuche 


wobei unter M, das störende Elastizitäts- 
moment zu verstehen ist. Für die Pendelauslenkung a ergibt sich dieses Moment aus den 
obigen Gleichungen zu: i 
M,=(&,—@)-R. 

Zur Bestimmung der Momentenfunktion M, = M(«a) muß die Funktion R-G= f(«a) einmal für 
starre, ein andermal für elastische Verhältnisse ermittelt werden. Dabei ist es vorteilhaft, den 
Schwerpunktsabstand « des Pendels klein zu wählen. Entsprechende Versuche wurden vom 
Verfasser für Schneidenwalzen und Unterlagen durchgeführt, die z. T. in der Form nnd den 
elastischen Eigenschaften den Versuchskörpern der in Abschnitt I angegebenen Schwingungs- 
versuche ähnelten. Alle diese Versuche ergaben Elastizitätsmomente in gleicher Art und 
Größenordnung wie sie bei den entsprechenden Schwingungsversuchen auftraten. Als Beispiel 
für die Wirkungsweise des elastischen Mittels wird im folgenden ein Versuch für den Fall 
„‚Starre Schneidenwalze auf ebener, elastischer Unterlage‘‘ wiedergegeben. 


Pendeldaten: 
Radius der Schneidenwalze #==1,25cm 
Pendelgewicht P = 3,5 Kilopond 


Elastizitätsmodul der Unterlage EZ = 0,6 Kilopond/mm!?. 


In dem Diagramm in Bild 7 sind auf der Abszisse die Schalenmomente @ - Rin Kilopondimeter, 
auf der Ordinate die Ausschlagswinkel « in Graden angegeben. Die verschiedenen Kurven stellen 
die Abhängigkeit der Pendelauslenkung « von dem jeweiligen Schalenmoment @. R unter 
folgenden Versuchsbedingungen dar: 


K uvvenl, Gersueht. 
Starrer Fall, die Schneidenwalze besteht aus Stahl, die Unterlage aus Glas. 


KaursviesLEHN\V ezisiuichhz2; 
Elastischer Fall. Die Schneidenwalze besteht aus Stahl, die Unterlage aus Hartgummi. 
Walzenlänge = 2: 0,8 cm. 
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Das Pendel schwingt nach jeder Einzelbelastung frei aus, dann erst wird bei entsprechen- 
der Momentenvergrößerung die neue Messung vorgenommen. Dieses Verfahren hat den Zweck, 
einseitige Verformungen des elastischen Mittels während der Versuchszeit zu vermeiden. Die 
Winkelablesung erfolgt unmittelbar im. Anschluß an die Freigabe der Momentenschale. Bei 
einer längeren Standzeit des unausgelenkten Pendels vor der Durchführung der Versuche ergibt 
sich eine ähnliche Kurve wie oben, der Elastizitätseffekt erscheint jedoch vergrößert. 


I ® 
738 I/ 
12° Momenten-Versuch (statisch 
= Pendeldaten: 
u Starre Walze auf elastischer Unterlage 2 
Pendelgewicht P=3,5 Kilopond 
10° Walzenradius r»125cm 


länge der Einfachwalze #£,=2x0,8cm=16cm 

länge der Dreifachwalze £43-6x08cm+4öcm 

Elastizitatsmodul e & 
der Unterlage. E=06 Kilopond Im m? 


Auslenkung a in Graden 
> = 


S 
° 
LU} 


(® 
6° ° o a er 
Zeichenerklärung: 

5a I Starrer fall Te 

I Elastische Unterlage | tn Di a 

A L-15 cm, Einzelbelastung 
z IT Elastische Unterlage 
£,=16cm, Dauerbelastung. | —————— 
3° Je Messung 2Min-Standzeit 
E ® IV Elastische Unterlage oo 

2° D £3=48cm,Einzelbelastung 


10 15 20 25 30 35:70" 
Schalenmoment 6.R.in 10°” Kilopond- meter 


Bild 7. Abhängigkeit der Pendelauslenkung vom äußeren Moment bei verschiedenen Versuchsbedingungen 


Kurvelll, Versuch3. 


Elastischer Fall. Pendeldaten und Art der Unterlage wie bei Versuch 2. Die Momenten- 
belastung wird fortlaufend, ohne Zwischenentlastung durchgeführt. Die Winkelablesung er- 
folgt bei jeder Messung erst nach 2 Minuten Standzeit des ausgelenkten Pendels. Nach einer 
Höchstbeanspruchung durch das Schalenmoment 3,8 - 10-3 Kilopond.meter wird eine schritt- 
weise Entlastung vorgenommen. 


Dieser Versuch zeigt deutlich den Einfluß der elastischen Nachwirkung. Bei längerer 
Standzeit des ausgelenkten Pendels wird das Elastizitätsmoment geringer. Die Spannungen des 
elastischen Mittels gleichen sich aus. Bei der Entlastung tritt bereits bei 7° Auslenkung ein 
Gegenmoment in Erscheinung. 


Kurve IVveN\ersu ch 


Elastischer Fall. Pendeldaten und Unterlage wie bei Versuch 2 und 3, doch wird die 
Länge der Schneidenwalze auf das dreifache vergrößert. Länge der Schneidenwalze = 3.1,6 cm. 


Die Durchführung der Momentenbelastung geschieht im Sinne des Versuchs 2. (Ent- 
lastung nach jeder Messung). 


Ein Vergleich der Kurve IV mit Kurve Il zeigt wesentliche Eigenschaften der Momenten- 
bildung. Bei kleinen Winkeln, im vorliegenden Falle etwa bis 4° Pendelauslenkng. weichen 
beide Kurven nur wenig voneinander ab, bei größeren Winkeln nähert sich Kurve IV mehrund 
mehr der Geraden I. Zur Ermittlung eines Störmoments aus Bild 7 wird durch die entspre- 
chende Winkelgradzahl eine Horizontale gezogen. Der Abstand des Schnittpunktes dieser Ge- 
raden mit Kurve I von dem Schnittpunkt der Horizontalen mit Kurve II bzw. IV ergibt das 
Elastizitätsmoment für den betreffenden Winkel. 


ai 
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In Bild 8 sind die auf diese Weise ermittelten Störmomente in Abhängigkeit vom Aus- 
schlagswinkel & für Versuch 2 und 4 dargestellt. Die Momentenwerte des Versuches 2 lassen 
sich durch die Funktion 
M, —=:10,8 + 10? le| a37:10°"Kilopond: meter ne > 

2 =70,8:10"3a-3.1-10°2a ? 
— 3,1-10-2- 02 Kilopond.meter, at Zertilmmen h 
die Momentenwerte des Versuches 4 » 
durch die Gleichung: 


M, = 10,8:10-3 a ar wo 
— 4,8 - 10=2 - a? Kilopond.meter - ee 
wiedergeben. z 


Beide Kurven sind im Diagramm 
von Bild 8 eingezeichnet. 
Die entsprechenden Rechenwerte 


nach der MomentengleichungdesFallesI » 


2 / 2 e 
Er DP"Fr6 1 ö %, FETT ER a u RE TE SE 
e >= Me ö* —— j* N — En Auslenkung in Graden 
Bild 8. Elastizitätsmomente bei den Versuch 2 und 4 
lauten: 


Für Versuch 2: M, = 10,8 : 10°. — 4,1 - 10-2 - a2 Kilopond.meter. 
Für Versuch 4: M, = 10,8 : 1073 - x — 6,4 : 10-2 a2 Kilopond.meter, 


wobei das Verhältnis ö/ö* wie bei den Versuchen des Abschnittes I gleich 0,9 gesetzt wurde und 
die Drucklinie bei Versuch 2 gleich 0,33 cm, bei Versuch 4 gleich 0,2 cm war. 


Dieses Ergebnis zeigt, daß die vereinfachenden Annahmen, die zur Enwicklung der 
Näherungsformeln gemacht werden mußten, im großen und ganzen bei bestimmter Ausschlags- 
begrenzung zu Recht bestehen. Es zeigt mithin, daß die Störung bei den angegebenen Ver- 
suchsbedingungen im wesentlichen durch die beschriebene Momentenwirkung des verformten 
Randmaterials hervorgerufen wird, und nicht etwa durch die ebenfalls mögliche Ausbildung 
einer wirksamen Abrollkurve. 


Zusammenfassung 


Bei Schweremessungen mit den gebräuchlichen Schneidenpendeln ergab sich nach Be- 
richten verschiedener Beobachter eine Verkürzung der Schwingungszeit bei sinkender Ampli- 
tude, auch nach Berücksichtigung der üblichen Korrekturen. Der Verfasser des vorliegenden 
Berichts konnte in einer früheren Arbeit nachweisen, daß eine Verkürzung der Schwingungszeit 
immer eintritt, wenn die starre Schneidenwalze bzw. Unterlage des Pendels durch gleichge- 
formtes, elastisches Material ersetzt wird. Versuche am Pendel mit Gummischneiden und 
Gummiunterlagen zeigten dies. 

Die entsprechenden theoretischen Untersuchungen führten zu dem Ergebnis, daß durch 
ideal elastisches Verhalten des Schneidenkörpers bzw. der Unterlage ein derartiger Eifekt nicht 
zu erklären ist. Die genannten Abweichungen konnten jedoch als Folge eines Elastizitäts-Stör- 
moments gedeutet werden, das sich für alle Versuche von gleicher Form und Größenordung 
ergab. Die Klärung der Ausbildung eines derartigen Moments, die mathematische und einfache 
experimentelle Erfassung desselben war Aufgabe der vorliegenden Arbeit. 


Nach den hier wiedergegebenen Untersuchungen erweist sich der angegebene Effekt in 
erster Linie als Folge der Trägheit des elastischen Mittels gegenüber Formänderungen. Das 
Material der elastischen Schneidenwalze bzw. der elastischen Unterlage erfährt in der Nähe der 
Druckflächenbegrenzung bestimmte Randverformungen, die im Rhythmus der Pendelschwin- 
gung periodisch belastet und entlastet werden. Da besonders in solchen Oberflächenzonen, die 
von der Druckfläche weiter entfernt liegen, das elastische Mittel den verhältnismäßig kurzen 
Belastungseinwirkungen während einer Schwingung des Pendels nicht momentan folgen kann, 
kommt es beim Aufschwingen vor der Walze bzw. in ihrem vorderen Teil zu Stauungen des 
bereits elastisch verformten Randmaterials. Die Folge ist eine geringe Aufwalzung, Vergröße- 
rung des vorderen Druckflächenteiles und somit die Ausbildung eines Elastizitäts-Störmoments. 
Bestimmte Feststellungen über die Normaldruckverteilung in der Begrenzungszone der Druck- 
fläche sowie die näherungsweise Berechnung der vorderen Druckflächenvergrößerung in Ab- 
hängigkeit von der Pendelauslenkung geben die Möglichkeit zur mathematischen Formulie- 
rung des Problems. 
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Als Näherungsformel für das Störmoment M, im Falle einer starren Schneide auf einer 
elastischen Unterlage, oder im Falle einer elastischen Schneidenwalze auf elastischer Unterlage 
ergibt sich die Momentengleichung 


u nee), 


wobei die einzelnen Zeichen die folgende Bedeutung haben: 


— Pendelgewicht 
r = Radius der Schneidenwalze 
ö = wirkliche Drucklinie | 
= — 0,9, ein Verhältnis, das bestimmte Eigenschaften des elastischen Mittels kennzeichnet. 


Durch eine geeignete Versuchsanordnung konnte das Störmoment in obiger Form und Größe 
bei Schneidenwalzen und Unterlagen aus Kautschuk durch statische Momentenversuche ge- 
messen und die Funktion M, = M(o) ermittelt werden. 
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Über Iterationsverfahren für Gleichungen 
und Gleichungssysteme. (Teil II)” 
Von Rudolf Ludwig 


In gleicher Weise wie in Teil I werden hier Iterationsverfahren für Gleichungssysteme unter dem @e- 
sichispunkt der Ordnung des Verfahrens behandelt. 


Following the methods of part I, iteration processes for systems of equations are trealed from the point of 
view of the order of the process. 


Dela m&öme manitre que dans Part I des methodes d’iteration pour des systemes d’equalions sont traitees 
du point de vue de l’ordre du procede. 


B Hacroaımei paöore, cIeAyA MeTonaM YacTu ], TPAKTYOTCH HTePALMOHHLIeE METOABI Als 
cHCTeM ypaBHeHHÜ C TOYEH CpeHuA TOPAAKA MeToNa. 


2. Iterationsverfahren für Gleichungssysteme 
2,1.’ Die Ordnung eines Iterationsverftahrens bei Systemen 
Ein System algebraischer oder transzendenter Gleichungen 
EEE er EEE 
läßt sich als Vektorgleichung schreiben 


(Fe) 
ie nn ni! 
F&) 
Ein JV 2) kann man ähnlich wie für eine Gleichung mit einer Unbekannten in der Form 
Ka a) (kan ne Eee ER 


ansetzen mit den Vektoren 


oh) e vr 
u = | ) und Mia [ .) 
u R n(£r) 


ade, *) Teil I: Z. angew. Math. u. Mech. 34 (1954) 8. 210-225. Formeln aus Teil I werden durch vorgesetzte / 
!) Vektoren werden durch klein, Matrizen durch 1 
leine, große deutsche Buchstaben b : 
Spaltenvektoren erhalten einen FRPSTER: Zeilenvektoren einen oberen Index. Zur Vereinfachung a a 
Fila... Sl) und CR] Ba 
x %r —_— Oo 
°?) JV = Iterationsverfahren, JVr£ = JV k-ter Ordnung. 


nenn emer. 1954 Ludwig, Iterationsverfahren für Gleichungen u. Gleichungssysteme (Teil II) 405 


Um die Ordnung eines JV zu definieren [2], geht man aus von der Abweichung t; des Vektors 
t. von einem Lösungsvektor TE: 
197 >= 7 — T Phrase 8 ie en a sdeikrdeh, re (1,3). 


Entwickelt man tx;ı in eine Ta ylorreihe nach tr, so ist bei Einführung des vektoriellen 
Differentialoperators 


0 
= 
ie Tal, hol Al HaPBRm plaıs AA 
0 
Om | 
PORE COEE JE Er BELGRE 3 Er BrRE-T RE) 277 Br SEEREEEREEE € Yu 


Ein JV wird nun von p-ter Ordnung genannt, wenn die Taylorreihe erst mit dem p-ten Glied 
beginnt, also 


Er WE DITE) Her muncden Harman win lb) 


Man definiert daher: 


Ein JV für ein Gleichungssystem I, = f(f;) heißt von p-ter Ordnung, 
wenn alle partiellen Ableitungen der Funktionen f; bis zur (p— 1)-ten Ordnung für einen 
Lösungsvektor £ verschwinden, aber wenigstens eine der p-ten partiellen Ableitungen 
für £ von Null verschieden ist. 


Speziell müssen also bei einem J Vırfür ein Gleichungssystem alle n?ersten partiellen Ab- 
leitungen der n Funktionen f; für den Lösungsvektor £ verschwinden. Bei einem JVrr müssen 
außerdem noch die n? (n + 1)/2 zweiten partiellen Ableitungen für £ verschwinden. 


22 'Iterationsvekfahren zweiterund.dritter Ordnung 


Für das Gleichungssystem (1,1) f(x) = 0 werde zunächst vorausgesetzt, daß sich der 
Vektor f (r) mit einer Matrix 9 (r) und einem Vektor g (2) auf die Form 


fr) = HH): — g(ı) = 0 ee ee 


bringen läßt. Wird für einen Lösungsvektor £ des Systems weiterhin vorausgesetzt, daß 
det H(£) = H(r) #0, so kann man (2,1) auf die Form eines JV bringen: 


rad lg). - = 0 nee rn 22), 


Da aber (2,2) im allgemeinen nur von erster Ordnung sein wird, soll nun ähnlich wie in Teil I 
der Vektor f(x) umgeformt werden durch Einführung zweier regulärer Matrizen B und Q, 
deren Elemente 9;; bzw. g;z zunächst als willkürliche Funktionen des Vektors r vorausgesetzt 
werden: 


FOPDITSTTFRIT - ana er, (2,3). 

Führt man ein 
SSR derer - - - - 2.4), 

so hat (2,3) wieder die Form (2,1), und es folgt daraus ein JV 

Ba) =) - - 2000: 2). 
Aus (2,4) folgt = 

8-5—-%9 
und ” £ 
g=Pg+9re—- P9Ee=-9r—- Pi. 

Also ist £ Ei 

Mer N PITN - te +: 6): 


Damit nimmt das JV die allgemeinste Form eines solchen an 


uau=t + fe), 
bei der die &, eine Folge regulärer Matrizen bedeuten. 
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Damit (2,5) nun aber ein JV von mindestens zweiter Ordnung ist, müssen für einen 
Lösungsvektor £ die Bedingungen erfüllt sein: 


ji up DE En ee 7, 


Wegen f(£) = 0 und unter der Voraussetzung 


det SEA) EN een ne (2,8) 


ist die erste Bedingung für jede reguläre Matrix Bund Derfüllt. Für die n zweiten Bedingun- 
gen folgt: 


Dw_a 0-5 - 


0%, 0%, = 0%, 
EN EN 
Da = — e, der v-te Einheitsvektor ist, gilt 
na 
5 0 
oder 


22-50... W=1,...,n.. (29). 


Faßt man diese n Vektorgleichungen zu einer Matrizengleichung zusammen unter Einführung 
des Matrizen-Differentialoperators 


ö 0 
BU (,:32) a (2,10), 


der auf einen Vektor f angewandt die Funktionalmatrix liefert, so gilt 


EUEDEETEIERDPIT] een al). 


Statt (2,9) kann man auch schreiben 


AL OR [OS ne . (&12) 


[SIR 1: = 5®] Aesken a ent Mr  WANSE 


Die Matrizen ® und DO sollen nun so bestimmt werden, daß (2,11) bzw. (2,13) erfüllt werden. 
Dazu gibt es — ebenso wie in I— zahlreiche Möglichkeiten, von denen hier einige angegeben 
werden sollen, die besonders vorteilhaft sind. 


a) Fordert man, daß die Matrizengleichung (2,11) identisch für alle zerfüllt sein soll, so ist 


5 BADEN a: - - Free 


bzw. 


und es folgt daraus nach (2,6) 


fe) =E— (RE) DEEDT Pa) fo): 
Die Matrix ® fällt also heraus, und das JV lautet: 


| zur = 20 — Dr tk] en | en, 


unter der Voraussetzung DE) #0, (IV). 


Dies ist aber gerade das Newtonsche Verfahren für Systeme (s.z.B.[1]). Aus (2,15) 
folgt für die k-te Verbesserung 


= rn — u = — Drlfle)] fire) 


nit die für die praktische Berechnung günstigste Form als lineares Gleichungssystem 
ür die d;: 


Del, Um). ...... (2,15%), 
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b) Wählt man in (2,11) von vornherein P(r) =€& (= Eistee so muß 
DE) = DLfe)] — HE) 


erfüllt werden. Diese Bedingung kann man näherungsweise dadurch erfüllen, daß 
man mit einem Näherungsvektor x, ansetzt: 


De) = Dlfleo)] — 9): 


Ho) = Hd) + Dr) = Dff(z,)] = const , 
und dam lautet die rechte Seite des JV: 


| 9 -r- ta |. 79 I PEE (2,16), 


und für die Verbesserungen bj; lautet das lineare Gleichungssystem 


| D [flto)] % + () = 0 |» 


Dies ist analog I (3,16) die bekannte Vereinfachung des Ne wtonschen Verfahren mit kon- 
stanter Funktionalmatrix [1], die sich durch wesentlich verringerten Rechenaufwand aus- 
zeichnet. Das Verfahren ist aber nur von 1.Ordnung, da (2,11) nur näherungsweise erfüllt 
worden ist. 


c) Geht man von (2,13) aus und fordert die identische Erfüllung: 


DIR) KO] = PR) HE) + DM), 
‚so lautet nach (2,6) die rechte Seite des JV 


| fe) = 2 — DU I) KO Bo) fa |, 


mit n? willkürlichen Funktionen. Dies Verfahren ist das Analogon zu 1 (3,8). Die Verbesserun- 
gen berechnet man aus dem linearen Gleichungssystem: 


| DIBE) F(Er)] dr + Bl) I) = 0 | I Pie 3 2 et 


d) Einen Spezialfall von (2,17) erhält man durch die Wahl von 
n 


Dann ist 


(J 149) Pe (2,16*). 


Eu ae (2,17) 


B— Fk Ss 
Dann ist 
0 a BR rR of 
Re ei: 
Mit der Matrix $* = (f,..., f) und der Diagonalmatrix 8 = (k;) ist 


[io =:- en+ Bil, 070 222.218), 


und für die Verbesserungen db; gilt: 


| (D[fE)] + Tr) KR) 0: + Für) = 0 * (IV)... (2,18%). 


Auch dieses Verfahren hat in I (3,9) ein Analogon. 
) Aus (2,11) ließe sich auch folgern 
De, a Keinen, N). 
Die Funktionen ?;, und g;; könnte man leicht aus den F,; aufbauen. Diese Festlegung wider- 
spricht aber der Voraussetzung (2,8), führt alsozu keinem JV. 
f) Die Bedingungsgleichung (2,11) läßt sich auch erfüllen durch 


Be) (DEI — IE) = Ar) — &) 


NR nl 2. (2,20): 


unter der Voraussetzung 


Dann ist 


S=BH+D=BOfI+E. 
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Setzt man noch BP = &», so lautet die rechte Seite des JV: 
id =r-(OI + ID | YIm). . .-- 21), 


unter der Voraussetzung DIE] +0, EE)=0. Wegen c;;(£) = 0 kann man die Matrix 
durch Linearkombination der F, erhalten, also z. B. 


a ee er (2,22) 


mit den Matrizen &,; = (0%), (,=1,...,n). Das Verfahren enthält also n® willkürliche 
Konstanten (bzw. Funktionen). Für die Verbesserungen b; lautet das lineare Gleichungssystem 


| (DLflz)] + Cxr)) dr + f(tr) = 0 I (U) ee 


g) Formt man nach (2,21) um: 
DD) +CL=- DeE+ INN, 
(DU +9t=(E+ DO)" DIN] 
und wendet formal die Neumannsche Reihe an: 
E+ N OT-E- DICH... 
beschränkt sich aber auf die ersten beiden Glieder, so ist 
fe) =r— (E- DO) DIflf 
=r- ff IN EMIT 


unter der Voraussetzung &(£)=0. Da & D![f] = — C* ebenfalls eine Matrix ist, die der 
Voraussetzung (2,20)entspricht, so kann man ansetzen — statt C*ist wieder & gesetzt worden — 
de. Arm)... (2,23). 


Für die Verbesserungen gilt: 


| Di) v,+ (E+ Ex) Fler) = 0 E (FF) 2 


h) Durch die frei wählbare Matrix C in (2,23) muß es nun möglich sein, das J Vır zu einem 


fe) 


JVınmzuerweitern. Es müssen zusätzlich noch die Bedingungen 3 = 0 Ian en 


%,0%, 
erfüllt werden. Man erhält nach kurzer Zwischenrechnung: B 
- je) _ aaa ee of 
ae, wo dx, dg, BLU 0x, 0x, 
o& of o& of 
no nn gemer.in. | 
+®% m(ze 32. + d2, “) ee (9. 
Es gelten nun die folgenden Beziehungen: 
Durch Differentiation von D-[f] D[f] = & ergibt sich 
SD] _ Zr 2 n- 
urT = D if] Br D if] ET u en (a). 
Weiterhin gilt: 
rg Of 
)) MN, — e, RES I N en Eu (B), 
E31 Org Ak graph en nn 
Ai 0x, 0x,’ RER uST DR, OT, m BE) SE (r) 
Setzt man («), (ß) und (y)in (*)ein, erhält man 
SEE) fe) , ICE) HE) _ FE) 
02,30, ä: 0, 0 ö7, 0%, Ws Gieschungns Sl, 


als weitere Bedingung für die Matrix ©. 


®) ® fällt bei der Rechnung heraus, wenn man setzt C = Ber, 


ss CI), 
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Wegen &(£) = 0 kann man die Elemente von & durch Linearkombination der Funktionen 
F;erhalten. Zu diesem Zwecke seienn symmetrische Matrizen &;=(d,)= (d,) ange- 
nommen mit den Zeilenvektoren e?: 


Re rei RE N ze 


Für die Matrix E gilt nun: 


a) &(£) = ist nach (2,27) evident. 
ß) Der letzte Summand von (2,23) lautet, wenn f’ der zu ftransponierte Vektor ist: 


6f- > Gr — (6) 


r&i 
also 
== 2 FF, „et hlr Aut. (2,28). 


Die p-te Komponente von f*ist alsoeinequadratisch e Form inden F, mit der Koeffi- 
zientenmatrix Q,. 


y) Es ist Kö ir: : 
f 


0%, 0%, 


Für die linke Seite von (2,25) ergibt sich somit: 
= or, of or, _ © 
a) = 2 0%, e 0%, % > Or, er 
Wie unter ß) folgt aus der Symmetrie von &;: 


EL EN 


,(E) = Di a lan, "on, 


of 
Er 


„a 
0%, 
Soll (2,25) identisch für alle g erfüllt sein, so folgt zur Bestimmung der Matrix @,die Gleichung: 


turen [289 


Die zweite partielle Ableitung von F, nach &, und x, ist also gleich der Bilinearformin 
den ersten partiellen Ableitungen 9f/0x, und 9f/dx, mit der Matrix &,. Die n (n + 1)/2 Ele- 
mente von @,lassen sich dann gerade aus der gleichen Anzahl der in den c?, linearen Gleichungen 


für ,»=1, :--, n bestimmen. Die c?, sollen zusammengefaßt werden zu einem Vektor c, der 
Dimension n (n + 1)/2. Aus (2,29) folgt also für p = const und „v =1,::-,n 


SS fo 1, a Ve) 
mit der von punabhängigen Systemmatrix D,[f] und den Vektoren: 
Be | @ 1 Bel ji) — N 
; ch, 2 nam | [Fa] 
also ist Br 
0,= Dr I 
und Er 
fs == CH r® Be Be ee Je (2, }: 


28 
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Das Lösungsschema für die n Gleichungssysteme mit der gemeinsamen Systemmatrix D;[f] ist, 


bei Anwendung des Matrizenverfahrens zur Lösung linearer Gleichungssysteme (Zurmühl 


[A)): 


2 2 Da 


Die Verbesserungen v; errechnet man mit den hiernach errechneten Werten für j* = Ef nach 
einem ähnlichen Schema aus: 
DEE] 9 + Fe) + TE) SO rn 33). 


i) Eine Vereinfachung des letzten Verfahrens kann man in der Weise erreichen, daß man 
die Gl. (2,29) nur einmal für den Anfangsvektor t, löst, also 


EP IE DERSR F u 28 a0, 


po 
0%, 0%, 2 02,0%, 


und diese festen Koeffizienten in (2,23) einsetzt, also: 


| DLfeEr)] dr + Fler) + Co) I) =0 |, IV) » » . » (34%). 


Dies ist aber nur ein JVı, da (2,29) für £ = t, nur näherungsweise erfüllt ist. 


Zusammenfassung: 


Für ein gegebenes Gleichungssystem f(z) = 0 werden JV von der Form tx }ı = f(t) 
aufgestellt. Die behandelten Verfahren werden in Tabelle I, Teil II (S. 415) mit An- 
gaben über den erforderlichen Rechenaufwand zusammengestellt. 

Vergleicht man die Verfahren für einfache Gleichungen und für Gleichungssysteme, 
so fällt die Gleichartigkeit sofort ins Auge, wenn man die Verfahren für die Verbesse- 
rungen v bzw. b angibt: 


Gleichungen: F(x)=0 Gleichungssysteme:f(rJ)=0 
JVn: (a) Fo+F=0 Q[f]p + = 0 (Newtonsches Verfahren) 
(b) (pFyv+pF=0 DB fID + Pf 0 
a SE 
TC))V® =, a 
(e) For Fri) PR = 0 Stflo De A Vor Al 
ar: (l) Fo+r +4, m=0 flo +f+ft=0 
mit p =f &,i . 


LER BE 
AT Ba 
Beispielzudem JVmm(f): Zu lösen ist das Gleichungssystem: 
F,=i1+2-235=0 


Fr, = =I=0. 


und 


Dann ist 
43 44% 2 
m 7), Vi-nm Mit ma 
2 ı Dr A 5 2 ) 
| ı[2 1 
m=-2AMF, M=>[0, e=->11|; 
Br 2 * [0 
4,5 0,75 0,5625 
2=(3 1%) -( ). Ma)=| 0 
0 
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Nach den Rechenschemata (2,32) und für das Ne wto.n sche Verfahren (JVır) (a) ist: 


[4,601 4 en td: 
JYm: n=(19s5 97) ZVe: &=(19538)- 


Wählt man als neuen Ausgangswert:}, = Ba ‚„soist: 


nachNewton(JVnm): u = He: 24 ee ; 
- __[4,601 594 920 
Bars lem J/Fır: &2 =(1.955 843 Ar 
nach dem vereinfachten z __[4,601 594 912 
Verfahren (2,34): & = (1,955 843 610) ° 
Der exakte Wert ist: r= Be 28 ns . 


2,3. Konvergenzbedingungen 


Es sollen hier Konvergenzbedingungen für JV angegeben werden, die einen Vergleich für 
Verfahren verschiedener Ordnung zulassen. Für spezielle Verfahren existieren teilweise schär- 
fere Abschätzungen, z.B. für das Newto.nsche [1]. Ohne auf die besondere Form des 
Verfahrens einzugehen, werden hier Konvergenzbedingungen hergeleitet, wenn man nur die 
Ordnung des Verfahrens kennt. 


Geht man aus von dem JV, wie (1,2), 


Deal Job nsyaemiill maters sb au (3,1) 


mit einer Lösung t, für die also gilt er = i@), so folgt daraus die Ta ylorentwicklung (1,4) und 
speziell für ein JV, (1,5): 


ler ni (ee (3,2), 


oder, wenn die % (k=1,2,---)und £in einem Bereich ® liegen und ein Vektor x diesem Be- 
reich angehört, so kann man dafür schreiben: 


ang =] be) ed. ee (3,3). 
Die g-ten Komponenten der p-ten Ableitungen seien beschränkt in ®, also 
ea Mer ee nn N. 
Dann kann man nach (3,3) abschätzen: 
a Si, | B3 a — 2, j re rt 7 (3,5). 
Setzt man noch T 
K, = 21:0 — 3, ie pe (3,6), 


so ist: 
(1 = p 
|a7’ 2%) S My. RK: 


a ae ei Dh ieh, * 


(k) | al RE | ® PET pETt Li p24n 
%q — | = Mo za — | Ss Meg K, Mo. 
—=1 = 


= M,,. vr 2 M Er = M,,. K > M ‚| 2. 
A=1 m 


28* 
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Damit diese Folge konvergiert, ist also hinreichend: 


Be Sem, ; Kcal. a. 
)=1 


Daraus folgen die speziellen Konvergenzbedingungen für ein 


2 « 
IVı: | Salt. ala 
JVn: K, - SM;; ZIERT. IN 20 EINER. Weyer, ER (3,9), 
JVm: K2 I. RR > (3,10). 


Mit Ausnahme des JVı kann man also nach (3,7) einen Anfangsvektor x, stets so wählen, 
daß das Verfahren konvergiert. Ein JVı konvergiert also dann, wenn es in der Umgebung von 
reinen Bereich ®* gibt, für den (3,8) erfüllt ist. 


2,4. Verbesserung einer Iterationsfolge 

Genau in der gleichen Weise wie für JV für einzelne Gleichungen (I; 1,5) kann man für 
Systeme auch einen verbesserten Vektor £ angeben. Die Herleitung soll hier nicht angegeben 
werden ®). Das Ergebnis ist: 

Aus den ersten Differenzen der Iterationsfolge der k-ten Komponenten werden die Vek- 
toren gebildet: 


[ 2 — ar) IEAR den ar an 1 | 
vn BO eretke (ME) ap selinnek 7 Bden JaLL N ae 
| zn — a | dy n 0 | 1 
Dann gilt: 
> det (d,,- ++, d,,09) 
re 2 En ER ı E 
| a det (d4,. --„Du, ©) (42) 


Die Berechnung dieses Determinanten- Quotienten ist sehr einfach bei Anwendung des Matrizen- 
verfahrens [4], da die Determinanten sich nur in der letzten Spalte unterscheiden. 


2,5. Iterationsverfahren höherer Ordnung durch Linearkombi- 
nation von Iterationsverfahren niederer Ordnung 

In ähnlicher Weise wie in (I; 2,5) soll gezeigt werden, wie man durchLinearkombination von 

JV 1. bzw. 2. Ordnung neue Verfahren zweiter und dritter Ordnung herleiten kann. Es seien p 


JV gegeben; 
Hiahe), Veh po. een 
mit der Matrix der rechten Seiten 
1 
= =(.,) a A (5,2). 


Mit Hilfe der Matrix der » : 9 ,, Gewichtsfunktionen‘“‘ 
= (4) =, all, Weil, +95 ke} an) 8 . (9,3) 


soll das neue JV 1 = f(xx) gebildet werden, dessen rechte Seite der Diagonalvektor des 
Matrizenproduktes 2% ist: 


: hf 
0 =(;.) Se 


*) Vortrag auf der GaMM-Tagung 1952 in Braunschweig, [3]. 
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Für einen Lösungsvektor £ muß gelten: 


Deu, re Dil ee Fans) 
und daher 
ee (ar) 
En In(E) & 


also müssen die Gleichungen 
Lo), =1 er; (5,6) 


erfüllt sein. D. h. das mit (5,4) gebildete JV ist dann wieder mindestens von 1. Ordnung, wenn 
die Spaltensumme der Matrix 2 für die Stelle r gleich 1 ist. Für ein JVı können die Elemente 
von Xauch Konstante sein mit der Spaltensumme 1. 


Damit das neue Verfahren von zweiter Ordnung ist, müssen alle ersten partiellen Ablei- 
tungen von (5,4) an der Stelle£ verschwinden. Für die k-te Komponente f; =!;(t) f*(z) muß also 
gelten: 


oc) _ A = j une) a 
et ee 
oder wegen (5,5) 


ed. 1 - De a en 


Fordert man, daß (5,6) identisch für alle r erfüllt ist, also 


LEN Eee ia ze 
so gilt auch 
ols(E) , 021,(2) } 
= —e = I RE 
027, ,=0; 0, 02, Az00 ;% (5,8%) 


Also wird wegen (5,8) und (5,8*) aus (5,7): 
= ME 
BD Ed 2 near 


Sind die Verfahren (5,1) von zweiter Ordnung, so ist (5,9) erfüllt, sind sie jedoch von 1.Ord- 
nung, so muß man auch hier wieder fordern, daß (5,9) identisch erfüllt ist, also 


1a) ED = 0 el a ER (5,10). 


Für ein Verfahren 3. Ordnung müssen außerdem noch alle zweiten partiellen Ableitungen an 
der Stelle T verschwinden: 


fe) _ LG) ee) A E) 8°) 


oz, 00,1) 0%, 0%, Or, 0%, 


ALLE) a),  _ 


0%, 0%, 0%, 00, 


En 


Aus (5,10) folgt: 


alle) fra) oe) _ 
ei 


olx(£) of*(z) Ban) ef (%) =) 


aber auch | dm, 0% 


ED, Oluoeh 0, 6. 0%, 


Wenn man die Identität fordert, so müssen erfüllt sein: 


0? f*(£) —(, a n) RE (5,11). 


(2) 02,00, 
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Es gilt somit der Satz: 


Aus denp JV Ları, = filte), = 1, , p) läßt sich mit der Matrix der Gewichts- 
funktionen 2= (A;,(£)) ein JV bilden tr} = Ru dessen rechte Seite der Diagonalvektor 
des Matrizenproduktes 2% ist, also fx(£) =1x(£) f"(k)- 

a) Für gegebene JVr ist un neue JV von 2. Ordnung, wenn die A;, die folgen- 
den Gleichungssysteme (k = 1, - : -, n) erfüllen: 

Ko)e,=1, l.(£) —— 2 = 0, Wels m) ra. 


Es sind also p= n + 1 JVınotwendig. i 
b) Für gegebene JVı ist das neue JV von 3. Ordnung, wenn außerdem 
die A;, noch die Gleichungen 


ae =0, wr=l,,n). . 618) 


erfüllen. Es sind also p = (n +3n+2)/2 JVınotwendig. 
c) Für gegebene JVır ist das neue JV von 3. Ordnung, wenn die A;; die Glei- 
chungssysteme 


=, = rl). (614) 


erfüllen. Es sind also p = (n® + n + 2)/2 JVır notwendig. 


Beispiel: (wiein 3,4): Zu den beiden gegebenen, Gleichungen des Systems 
F="—,=0, R=3ı, — 3=0 
fügt man eine dritte hinzu: 
FI, =. —-F,=o2(,—3)—0. 


Daraus wird die Matrix der JVı gefolgert: 


[ ae See ae 
a 2 
Ei % % 
= 
ee 
% © 
mit den Ableitungen: 
03 23 ae | 5 I 3° x, «8 
| en a 
beallirgigaan nie alla Bl {an ala DEE Fa 
| U, 2 | | X, | 
Es müssen für die A;z die Bedingungen (5,12) (k=1,2;v =1, 2) erfüllt sein: 
AıtAatAı el Aa t/a tie = 
9% Aa _ = 632 
Bar » 12 ut, IE =. 
Ayı 2 Agı 
2 t3 "An u 3 = 
Daraus folgt für die rechten Seiten des JV: 
= 15 2 = 7 


= Id 637 h=ar2 


| [eu rigen De 
A; = %— De 8) mn“ Ra) ie "J 


5 | 41 = 


—(® Ne) + a) N 


ee Ge Be en un a le en u E = er 


Lg — DEI (FE) + Cry) fe) 
Vor.: &(g) = 0 
Gr %— I Lie) (Hey) > 3)) 
mit =ißi, 


WA L_er 
Oxu Po = Oxu 0x, 


DI [HE] (Hr) + Er fley)) 


HE) ME) _ 


E- id 0° F(xo) 
da 0x 


2 Oxu 0% 


m |__________ DL en en Ba un en ee ee nn 


k 
. 82,89) 
det (dı ye..y dns sp) 


geg. JVI: Di — ii = EM 
&rrı ” if) mit F(t) = L,® f* (x) 


ll II aus: |x(X) fg 1 
ED _ e 
ke, = 0:7 aM an 
ee ehli)si=ls...,p 
| Errı — ee 
mE, p=ynm +3n +2) 
ojk = 021% (1) 
12 = u NEL 
II |kW)e=]1, 1k(ı 0, 14(£ uhr Fr 
b) Jun, p= m +n+2 
. ee. 


ı) F, = Zahl der Rechenschritte zur Berechnung der Komponenten von ®. 
2) C; k= Zahl der Rechenschritte zur Berechnung der Elemente von ©. 


et) I + Felt) Yin) ar en 3 Far EN ah | 


ER ee Eu 


i=1 


ei! 


=r zn (ms + 3n—1) ?) 


Euer 


+ m +60) 


Er+ Eir+t & Fan 


LT k= 


zn mt DM tn + Im 
+18n +8) 


ne +3n—]) 
1. Schritt: wie 8 
2. u. folg. Schritte: 
Rn 
Tre nm +1) + zn +3n—) 
i=1 


Zn dm +9n +09) 


B3 


Bi formelmäßig ausgerechnet! 


n 


ZF; 


il 


meist formelmäßig ausgerechnet! 


dsgl. 


?) Zahl der Rechenschritte zur Berechnung der Verbesserung, nicht der Iterationen. 


5,12 


5,4 


5,13 


u  Ee 
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Geht man aus von dem Anfangsvektor £,= ei so erhält man f, = (1'928) undßg = 1,9 3313) 


mit einem Fehler „= 1, — i = — 10° (5 : 


Es kommt bei dieser Methode darauf an, durch geschickte Umformungen die JVr so 
aufzustellen, daß sie möglichst einfach sind und sich auch leicht differenzieren lassen. 


2,6. Rechenaufwand für die behandelten Iterationsverfahren 
— Tafelder behandelten Iterationsverfahren— 


In gleicher Weise wie in Abschnitt I; 1,8 sollen hier die behandelten J ” zusammengestellt 
werden mit Angaben über die Zahl der erforderlichen Rechenschritte. Dabei wird die Zahl der 
Rechenschritte zur Berechnung eines Funktionswertes F,(g) mit F; und zur Berechnung einer 
Ableitung AF,/dx, bzw. 02F,/0x, 0x, mit F;, bzw. F;,., bezeichnet. Zur Lösung eines linearen 
Gleichungssystems mit » Unbekannten soll, wie bereits an gewissen Rechenschemata gezeigt 
worden ist, vom Matrizenverfahren [4] ausgegangen werden. Man berechnet leicht, daß dazu 


m=anm+3n— 1) 


Multiplikationen und Divisionen nötig sind (ohne Rechenkontrollen). Hat man außerdem 
noch das System mit derselben Matrix aber anderer rechter Seite zu lösen, so sind dazu nur 
noch n? Rechenschritte nötig. 


Zusammenfassung 


. Man kann also auch für Systeme zeigen, daß es sehr verschiedene Möglichkeiten gibt, JV 
bestimmter Ordnung aufzustellen. In weit stärkerem Maße als bei Gleichungen mit einer Un- 
bekannten machtsich bei Systemen allerdings das Anwachsen des Rechenaufwandes bei Erhöhung 
der Ordnung bemerkbar. 
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Konforme Abbildungen durch elliptische Funktionen 
Von E. Fick in München 


Es werden in der folgenden Arbeit die konformen Abbildungen dargestellt, die die Funktionen &, 


C, 0, 0,, sn, cn, dn, am, E, zn, d, und ®, bei rechteckigem Gitter liefern. Die Kurven wurden 
numerisch berechnet. 


The author represents in this paper the conformal mappings that are given by the functions 0,6, os, on, 
en, dn, am, E, zn, ®,, and ®, in the case of a rectangular lattice. The curves are computed numerically. 


Dans Varticle suivant les representations conformes fournissant les fonctions $, £, 0, 03, sn, cn, dn, am, E 
zn, d, und d, ü grille rectangulaire sont presentees. Les courbes furent calculees numeriquement. 


’ 


B Hacroameh paÖbore IPeACTaBıAmTcH KOHdOPMHEIe OTOÖPasKeHHa, Narınne ÖyHkoun 
p, & 0, 05, sn, cn, dn, am, E, zn, d, uU d, B cıyyae IPAMOYTONIBHOH pemerku. Pacyer KPHBBIX 


ITDOH3BOAAJICH YHCJIEHHBIM IIYTeM. 
1. Einleitung 


Die durch elliptische und damit verwandte Funktionen vermittelten konformen Abbil- 
dungen sollen im folgenden wiedergegeben werden. Die funktionentheoretischen Sätze und 
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Formeln dieser Funktionen können als bekannt vorausgesetzt werden [3], [6]. Sie bilden den 
Ausgangspunkt zur Ermittlung der Abbildungseigenschaften. 


Es wird ein rechteckiges Gitter zugrunde gelegt. Die Halbperioden sind ® = 1,0 und 
o&' = 1,5 bei allen Funktionen, mit Ausnahme der zn- und d-Funktionen, bei welchen sich 
zeichnerisch ungünstige Werte ergeben würden?!). Aus ® und w’ lassen sich die übrigen Para- 


' 


meter wie e,, n, k, K, usw. über 9, (0 = 2) bestimmen [3]. 


Im allgemeinen wurde zunächst der Verlauf der Funktionen längs der reellen und imagi- 
nären Achse berechnet ?). Mittels der Additionstheoreme läßt sich dann jeder weitere Wert der 
Umrandung und des Inneren des Periodenrechteckes bestimmen. 


Der Werteverlauf längs des Randes wird im folgenden dadurch wiedergegeben, daß dieser 
in genau angegebener Weise längs einer Achse aufgewickelt wird und die zugehörigen Funk- 
tionswerte als Ordinaten auf diese Achse aufgetragen werden (vgl. z. B. Bild1). Die Form 
dieser so entstehenden Kurven läßt sich leicht aus den bekannten Eigenschaften der Funktionen, 
wie Lage der Pole, Nullstellen, Verzweigungs- bzw. Kreuzungspunkten (= Horizontalstellen 
dieser Kurven), sowie aus Stetigkeit und Monotonie erschließen. 


Bei der Darstellung der konformen Abbildungen sollten insbesondere auch die Zusammen- 
hangsverhältnisse der Anschauung zugänglich gemacht werden. Um dies zu erreichen, wurden 
einzelne Blätter der Riemannschen Fläche nebeneinander gezeichnet. Sie sind längs der 
Verzweigungsschnitte aneinandergeheftet zu denken. Die beiden Ufer der Verzweigungs- 
schnitte, sowie ihre Bilder in der zugeordneten Fläche sind getrennt wiedergegeben. Die Zu- 
ordnung der aufeinander abgebildeten Flächen wurde sowohl durch verschiedene Strichstärken, 
als auch durch Pfeile mit Numerierung deutlich gemacht. 


Um die Zusammenhänge der einzelnen Funktionen untereinander wiederzugeben, wurde 
ein Viertel des Periodenrechteckes durch einen Schraffur-Saum hervorgehoben, und der 
entsprechende Teil in jeder anderen Ebene ebenso charakterisiert. 


Real- und Imaginärteil werden — mit Ausnahme von o und d# — durch die Indizes 1 
und 2 dargestellt: 
E= fiseile: 


2. Die Weierstraß sehe p-Funktion 


Mittag-Leffler sche Partialbruchdarstellung: 


1 , 1 1 
p (u) = >] Si a; 2] ’ 


wobei: 
oo=2no +2mw 
DGL.: 
3 
= 46% 99 Gg=4 II (6 —e). 
ü=1l 
Umkehrung: 
3) 
U= j Se 
a Me®—n9—-9 
(Jg nach Weierstraß). Rechtechiges Gitter. 
w reell (w = 1,000) 
Additionstheorem: @' rein ımaginär (@'=1,500i) 
1 ll — 6 (po) 12 e, reell (@-+1,648 
bu +n=-- ed — pa) - pß). : 2.0647 
&--1001 
Zusammenhänge: 9,9, reell (4-+8272 
F= Ve, —&U 44269) 


4>0 


p(u)= 8 + (& — Es) ns? F. 


Bildi 


1) In diesen Fällen wurde das Periodenverhältnis = 0,75 i angenommen. 
2) Es dienten dazu z.T. die Tafeln vonMilne-Thomson[4],in welchen die Ja co bi schen Funktionen 
n, cn und dn für reelles Argument tabelliert sind. 
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Bild 2 


Bild 3 


Bild 4 


ee En 
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3. Die Weierstraß sehe £-Funktion tem 


U 
zw | |pwzla vin= sm 
0 


d 
a | 
49-96 —-9 | 
(Jırnach Weierstraß). / 
| 
Additionstheorem: 
I1pW-— = (v) 
Zusammenhänge: 
= EF—iK,h V jr -,—W—o)+n Rechteckiges Gitter. 
Fe 1, Bw) . n reell (n=0,821) 
77) v (v, T) $ n‘ rein ımagınär (n'’=-0,340i) 
mit 
ee. u er w 
ee ne 


Bild 6 


Bild 7 
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4. Die Weierstraßsehe o-Funktion 


77 
[kw-4lau 


(u) = we? ; wi. 
Additionstheorem: / | 
ou Fv)o(u —v)= 0%(u) 0!) [pl) — pw]. 
Zusammenhänge: 
n 
32 "d,(w, 7) 
ou) = 2we 09 
Bild 8 
Bild 9 


I 


a | 


Bild 10 


Mech 


1954 
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&(u) 


d. Die o,-Funktionen 


(wa + u) ze N4% 


ox(u) = o(u)p (u) — ex = 


(0x) 
u ,?T), 
a Ra 8. x 1 
; Den (0, A ‚ 228) | 
(Nach [6]; in [3] sind die ox anders numeriert.) d 


In den Bildern 11—13 ist die Funktion o, (u) wiedergegeben. 


— _clu) 


0, (u) = ln nit F=Yye — u 


Bild 12 


Bild 13 
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ee seen 


Kreelle#,(K=1628) 
Krreell2 £, (K=2442) 
OSk27,  (k=0366) 


S 


%: ds 
aa 


‚BER VE VE Den Be see Il.) 


nn un un 


@ > 


eiı_ı 


6. Die elliptischen Funktionen von Jacobi 
s=sn(F, k) f 
ds 


Ss 
Umkehrung I (s, k) = ITmermern 
1— 82) (1 — %? 8?) 
zii ( 


(Jrnach Legendre). 
Additionstbeorem: 


dnG GenFdnF 
nee sn F cn Gdn@G + sn | 


1 — k2 sn? F sn? @ 
sn( F,k) =isc(F,k’) 
Vollständiges Integral: 
F6=1,h)=K(k)=Yea -,0= 5 920,7). 


Zusammenhänge: 
ılaa _ Ye) _ 1 dor), 
lie ing: ee Van Y%9,0% 
sn"'F=cenF-dnF; 
_ 9%(0, 7) K 


k X — VI; K'(k)=K(W)=—iVe —e0; T=i cc 


0,7)’ 


a Ds Bild 15. 


% (Die gestrichelten horizontalen Geraden 
ergeben in den s-Ebenen Kreise; 
vgl. [1, 8. 278]) 


Bild 16 
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== en (F, k) 


Umkehrung: 1 
ae er: 
YA) @”+ Re) 
Additionstheorem: 


en(i F,k) = nc (F, k’) 
Zusammenhänge: 


nF, = VER a _ Aw) _ 1 Ho 
pw) —e 0,(u) k d,(w, ?) 
en" = —snFdnF; ®ftsnFr=1. 


Bild 17 


a 
I WI- RZ cd) 
Gyar:} 1° fa 
! 
Peer DB 
Ta) a5’ 


LOBPFEP-SPHR 


Bild 18. 


(Anstelle des zum primitiven Periodenpaar 

K,2(K +iK’) gehörigen Perioden- 

parallelogrammes (dicke Linien) kann man 

das Periodenrechteck mit den nichtprimi- 

tiven Perioden 4K, 4iK’ verwenden. — 

Weitere konforme Eee vgl. [1, S. 281 
u. 


d=dn(F, k) 
Umkehrung: 
1 
dd . 
F(d, k) — (rn 
d 
Additionstheorem: 
— k? 
Se FL ET I 


1 — k2sn? F sn? @ 
dn(i F,k) = de(F, k’) 
Zusammenhänge: 
hm) -& aW) wm) 
dn"F=—k2snFenP; dnF +l2sn?F =1. 


dn (F, k) = DW) —e; 93 (u) Z d (v, 7) 


Bild 19 
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Bild 20. (Vgl. auch [1], 8. 284—285). 


7. Die Jaeobische Amplitude 


o=am(F,k); sinp=sn (F,k). 
Umkehrung: 


Y 
dp 
F (9, k)= == 
| Y1 — k2sin? 
) 
(Jı nach Legendre). 
Zi Zeillridt Zusammenhänge: 


on(F,k)=cosp; dan (F,k)=Y1 — k?sin? p; 


‘ am'’F=dnF. 
1 
| 
ıLL® Bild 21 
ORFr2i” F 
Fri Lt 
& d 
F= P m 
-/ m k’sin’p nn 
[0] 


Bild 22 
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Bild 23 - 


8. Das elliptische Integral 2. Gattung nach Legendre 


F 
5 dn?(F,k)dF, 


E(F,k) mit = sn(F,k)=sinp 
Q 
— k? 2 une Lu n 
E(s, k) = | -- —-ds, E(p,k) = fr — k?sin2pdp. 
0 
Additionstheorem: 
E(F +@G)=E(F) + E(G) — k2sn Fsn@Gsn(F +6) 


B(iF,k) = i[— EiF,k%) +F + sc(F, k’) dn(F, k’)]. 
Vollständiges Integral: 


an um me, ao tn 
Eß=1,k) =E(k) -42 77; 6%) =El) ae IT 
e) —& Ve, — & 
Zusammenhänge: 
1%?) su trtwto)— 
A a u 
) 2K9,W, 7) en 
Bild 24 
Bild 26 SE OTT Tr 
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p=amf 

19 

s 

441 

N 

Ka 2 

Ye ee see mn 

Bi | mw ı Il 

Iol_1P IL 1 19 

BE 2 es EEE u 

1 | ot B. ja I mi E) 

mn jan, i | 

Kette eng 
FHrteogı! 


r A H | {) 


„Rechteckiges Gitter: 
K,KE,E' reell 

(K = 1628) 

(K' = 2442) 

(E = 1517) 

(E'= 1131) 


© © 
Era 
® 
Ne 7 \L- #) 
\ L 
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Bild 26 


Rechleckiges Gilter: 13 


i er 3 Se 9. Die Jacobische Zetafunktion 


znF IT O8: E(k) 


aM =ERN- mE: 


r Additionstheorem: 
zn(F +G0)=znF+zn@ —k?snFsnGsn(F +0), 


mi P,k) = i[sc(F, k’) dn(F, k') — N — zu (F, Ir)] 
- Zusammenhänge: 
1 %,?) 
N ee 
FETT den 


(In den bisherigen Bildern 1—26 wurde 7 = 1,5 i angenommen; in den 
Ki folgenden Bildern 27”—34 wählen wir = 0,75 i.) 


Bild 28 


be a u a 
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10. Die Thetafunktionen 
(vr) = 2 ga > (-YPARrFtD sin Ant) mv. 
nn n=0 
int ‚ 


mit q=e 


letz). 


Bild 30 


© 
Sr EIER 
Mi 
U 


ar 
2 


Bild 81 


| 

I 
) 
run. 
{una \ol 
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Ka T=071 


Bild 29 


u 


3 5 " h. 
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dw.) =1+2L% (— 1)"g" cos (2anv) 
— n—1 


Bild 32 


Bild 33 


a 


| | | Ä 
i | je | | Bild 34 
| DR BER Te a2 2 Mi 


ricom i: Eiliptische Fe | un. 1048. 
 Eing gongon am 12.11.1958. 0000 DE, 


= KLEINE MITTEILUNGEN 


Ein Satz vom Kreis. - strecken und Halbmessern gleich der Anzahl der Teil- - 
N, Wird der Umfang eines Kreises vom Halbmesser 4 Punkte, also (Bild 1) 
j in n gleiche Teile geteilt und ein Teilpunkt mit den üb- (n—1) u, 
rigen (a — 1) Teilpunkten durch Gerade verbunden, UI =- Pa 
so ist das Produkt der Quotienten aus den Sehnen- ne gr” ’ 


Für A= lwird 


Nach Bild list 


Be _ 
Eee _ 24- = “ 
= ey sin; . 
“ 0% IB 
. 83 2 A -Birl ne 2 Asin en, (2). 
nm —24-sin 1924 u 
Unter Berücksichtigung von (2) geht Gl. (la) über in 
(n—1) 
j 201, I] eu =n 
< (= 
Bild 1 bzw. De (ih): 
(n—1) 


Bild 2 
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Zu diesem Ergebnis gelangte der Verfasser zufällig 
auf empirischem Wege bei der Bearbeitung seiner 1927 
der Technischen Hochschule Berlin vorgelegten Disser- 
tation [1]. In dieser Arbeit wurde u.a. für den Sonder- 
fall der Grundrißanordnung einer Anzahl r gleich weiter 
und tiefer Absenkungsbrunnen, in gleichen Abständen 
auf einem Kreisumfang verteilt, eine schnelle Lösung 
der vonPh.ForchheimergefundenenAbsenkungs- 
gleichung gesucht. Für beliebige Verteilung der Brunnen 
im Grundriß lautet die allgemeine Gleichung von 
Forchheimer ([2], [3]), unter Berücksichtigung 
der Reichweite Rsowie der Spiegelhöhe 4 des ruhenden 
Grundwassers über der grundwassertragenden dichten 
SohlenachKyrieleis-Sichardt[#], 


mg UmrInz,2 2...m) 8). 


| 
PTFTTTTTTTTR 


En 
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Für das unter Berücksichtigung von Gl. (2) um- 
geformte Produkt . i 
e zoRr . (n-D)r 
81 82... 901 = Wr) sin Ze ae 
n—1 2 
— 2n—1 sin — 
4 c=1 3 
wurde für verschiedene Werte n das Ergebnis errechnet 
und in einer Tabelle zusammengestellt, wobei ersicht- 
lich wurde, daß die Lösung gleich nist [4]. Leicht ana- 
lytisch zu verfolgende Beispiele (n = 2, 3, 4, 6, 8 usw.)- 
führen zum gleichen Ergebnis. 
Man erhält so aus Gl. (4) die viel bequemere Um- 
formung [#]: 
N. e 1 A 
B-M--(mR-mA+— In.) 
ck , Rn 


er DR (4a). 


Halbmesser r 


n-1 

Hierin bedeuten gemäß Bild 2: 

H = Spiegelhöhe des ruhenden Grundwassers über 
der grundwassertragenden Sohle, (m), 

Z = Standrohrspiegel über jener Sohle in einem 
Beobachtungsbrunnen, (m), 

n= Anzahl der Entnahmebrunnen, 

= Br je Brunnen in der Zeiteinheit 
m/s, 

k = Bodendurchlässigkeitswert nach Darcy, m?js, 

R= Reichweite der Absenkung, m 

% %9,%g....%n = Abstände der einzelnen Brunnen 
vom Beobachtungsbrunnen. 


Für den erwähnten Sonderfall (Bild 3) wird der Be- 
obachtungsbrunnen am Außenmantel eines der Ent- 
nahmebrunnen stehend angenommen ; Gl. (2) geht dann 
überin 


3 


nq 1 
H? — = —ı In R Zn 81, 89, 83 ...8n-1 + Inr| (4). 


Bild 3 


Infolge starker Inanspruchnahme durch 
Baupraxis und Bearbeitung technisch- 
wissenschaftlicher Veröffentlichungen kam 
der Verfasser bisher nicht dazu, die mathe- 
matische Fachwelt auf die durch die Pro- 
duktenreihe Gl. (1b) ersichtliche Eigen- 
heit des Kreises aufmerksam zu machen. 
Ein ihm von einem Mathematiker der 
Siemens-Werke zur Verfügung gestellter strenger Beweis 
für Gl. (lb)ging leider durch Kriegseinwirkung verloren. 
Prof. Sza b 6 (Ord. f. techn. Mechanik a. d. Techn. 
Univ. Berlin) gab den Beweis für die Richtigkeit der 
Gl. (1b) an, der nachfolgend zusammengefaßt wird !). 
Zu beweisen ist, daß 


m 2 CH n 
N 7 sin—— —— [vel. Gl. (1b)]. 
n AR— 
c=1 
Man erhält mit 
ion ion‘ 
eH 1 EI 
Eine lu a 
n 23 
in ion 
M N 
ee ei e TER er ih zn] 
CR 1 1 Bi 
sin — = — (tt — 17) = — Eike 
n 2% 2; ze 


t) Der Verfasser dankt Herrn Prof. Sza b öfür diese Hilfe,ebenso 
Herrn „Prof. Dr.-Ing. Lorenz für seine freundliche Vermittlung. 


Aa nd 
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nach einigen Umformungen 


.n—l en 1 n—1 
11 a 1] ME), 
c=1 c=Ii 
wobei 
Face 
=e " 
ist 


Ferner ergibt sich mit 
un—1]= ei2nc; 
n_ an te 
w= w.—ye2rc=e "— \e *) =a82e, 


wn — 1 = (w — wo) (w — wı).... (w — wn-.,) 
(w — 1) (w — 2?) (w— 2°)... (w — x2(n-D)) 


l 


und 
un—] 
DT wiltw24t..+w+l 


= (w— 2?)-(w— @).... (w — 22: (n—l)), 
Mit w= 1 wirderhalten 
= (1—-22)(1 — x)... (1 — 22(n-h) 


oder 
n—1 
Te 293 =n...:.:. an 
c=1 
(II)in (T)eingesetzt liefert schließlich 
n—1 ee F 
II sin Pe [vgl. Gl. (1b)]. 
t= 


Es sei noch vermerkt, daß Gl. [1b] mit Bezug auf den 
Kreis besagt: 

Wird der Umfang eines Kreises vom Halbmesser 
„Eins“ in rn gleiche Teile geteilt und ein Teilpunkt 
mit den übrigen »n — 1 Teilpunkten durch Gerade 
verbunden, so ist das Produkt der Sehnenstrecken 
gleich der Zahl der Teilpunkte. 


Bremen. W. Sichardt }. 
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Zum Wärmeaustausch zwischen einer geheiz- 
ten Wand und einer Konvektionsströmung. 

In [1], [2] sind für verschiedene Randbedingungen 
Lösungen der partiellen elliptischen Differentialglei- 
chung 


En (Tex En Toy) = Y . 2 a (1) 


bzw. 


mit 
— (7 2)" ER „= (5 al di | (2), 


angegeben, die unter gewissen vereinfachenden An- 
nahmen für das Temperaturfeld über einer durch eine 
geradlinigeCouetteströmung gekühlte Wand gilt (T Tem- 
peratur, A Wärmeleitfähigkeit desströmenden Mediums, 
v seine kinematische Zähigkeit, o seine Dichte, c»seine 
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spezifische Wärme bei konstantem: Druck, 7 Wand- 
schubspannung). Der Couetteströmung entspricht die 
zur Wand parallele Geschwindigkeitsverteilung 


er: 5 0 Yy . 
Mittels einer einfachen Transformation sind die dort er- 
haltenen Lösungen auch anwendbar für das Temperatur- 
feld über einer geheizten Wand, die durch eine Strömung 
mit der Stromfunktion 


y=0,79 . DU5 .y-05..42. y2 (Yyy Zu, — Yo) (3) 


gekühlt wird. (3) stellt dieerste Näherung für die Strom- 
funktion eines Grenzschichtprofiles zur Potential- 
strömung U, = D.x°/? dar. Dieses Profil tritt z. B. in 
der Anordnung von Bild lauf [3]. } 


675° 
X 


Bild 1 


Mit (3) erhält man bei homogenen isotropen Stoff- 
werten und Vernachlässigung von Reibungs- und Kom- 
pressionswärme im zweidimensionalen Fall für das 
Temperaturfeld die Differentialgleichung 


C-(Tıa+Ty)=y-22-Ta—y?-x-Ty (4) 


mit 
rent, I ET (5) 
(4) geht mittels der Transformation 
in I (6) 
über in die Differentialgleichung 
ee 


die mit (2)identisch ist. 

Esist jedoch zu beachten, daß die in [2] formulierten 
Randbedingungen für y=0Ound z<0 und y=(, 
x <{ 0, sich auf die negative und positive x-Achse be- 
ziehen, die entsprechenden Randbedingungen bei (7) 
für 6=0,0<0 und i=0,0>>0 dagegen auf die 
positive y- und positive x-Achse. Es werden also in [2] 
Lösungen für die obere Halbebene, bei (7) jedoch 
Lösungen im 1. Quadranten untersucht. 

Abschließend sei an einem einfachen Beispiel demon- 
striert, daßtrotz der mathematischen Analogie zwischen 
(2) und (7) bei den physikalischen Ergebnissen wesent- 
liche Unterschiede existieren können: 

Im Fall der isothermen Wand ergibt sich bei der 
Couette-Strömung nach [2], Formel (39), wenn 


N 
Ba 
für den Temperaturgradienten an der Wandoberfläche 


in » T\L/8 
IT yo = — 0,538 4er = ‚aus... (8). 


.y 


F Er 
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Path of a vortex round the rectangular bend 


of a channel with a uniform flow. 


1.Vortex pathsin the L-plane. 

Let the channel with rectangular bend have the 
opposite sides parallel, that is Au B|| Au D and 
BC„ || DC» and uniform breadth h throughout;let U 
denote the velocity of the undisturbed uniform flow. 

BythetheoremofSchwarzand Christoffel 
the interior of the channel can be transformed into the 
half plane above the real axis of the ö-plane by the trans- 
formation 

dz _ 

de 
where k, A are constants to be determined; the corre- 
spondencesin the two planes are 


Ar ir re ae Sie a ec Te 
Bot=h; An >t=to. 

The uniform flow in the z-plane may be supposed to 
be due to a source of output Uhat A„ anda sink of 
intake Uh at C„. Hence, in the Z-plane, at the origin 
there is a sink ofstrength Uh/n (since the sinktakesin 


the volume Uh per unit time over an angle n). This 
flow is given by the complex potential 


Uh 
“Rs REBEL.) 


whence 
dr _ do de _ Un KA 
a dk nn ke-+tisn' 
d 
But, #4, CE = + ®) zZ =U,if the real axis 


(of z-plane) be parallelto BA... 


eh whence k = 
sch 


al=, 


er Re). 


d 
Againat Cx(E=0) = iU,if BO be taken as 


imaginary axis. 


1 Uh Uh 
ZN NEE SEE 
£ ck "mr Ben 


ker + Ya dyar 2.) 


= 


whengn er Fe Ar . (1.4). Aso D 


Therefore, (1.1) becomes 


a_h +48 nr 
En tem m TO 
ee (um 
where f”(£) ee ni 


Byintegration 


to rc, where f(&) = cosh-!&+sec-1L, 


use 2 
3A -SEH., j 


ah 


"PLANE 
Fig.1 


S 


-—————|um.. 


...__- —-. 


[7 
Eh 


© PLANE 


Fig. 2 


- .-. -..-. 


Therefore, if uz, vg denote the components of velocity 
of the vortexin the Z-plane, 


u + in 


KENN 
4n 


dw, en? 
sehe Tat 


Gt, 1 


22060 — &o) a=i] 
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Writing &, = oei?, and then equating real and imagi- 


nary parts . 
Uh cos 
a 0 ‚r 
| 
4rlosin @ Di 
ı. Uh sing 
vu = ne e 4 
re o?cos2p — 1 | 
4nl eo ot — 202cos2p +1! 


To obtain the paths in the {-plane we further calculate 
the two velocity components 


= urcosp+tmvwsingp 

2_0Uh,T [cosp (0 +1)sin@ | 
70 a ee Tr Sl 17 
0P = — using +vrcosp f A 

_T,__(e@—1)eosp 

An gt —2p%cos2p+1' 


from which by division and simplification one obtains 
cos 1 d 2 +-2ocospo-t]|1 
Amer erreg 
This on integration gives 
e+2gcosp+l] 
0 —2ocosp-+1! 
—4Ukhp=cons.. .... (1.82). 
If 0. 9, be the initial values of oand prespectively, 


then the above equation can also be thrown into the 
form 


= „ Er, (ei — 29005 + ai 
0? — 20cosp-+1l 08 +20,c0osp, +1 
P—Pe 


ish A 
sin 99 


—4Uh=0. 


T.|iog sinp+ log 


here «= ——.. 

where & Un 

These give the entire family of the paths of the vortex 
in the Z-plane. 


2. Vortex paths in theL-plane by 
graphical method. 


As the inverse transformation of (1.6) in the form 
E = Fe), or oeip = F(r ei®)is too complicated, a direct 
transformation of the vortex paths in Z-plane to paths 
in z-plane is not so simple. We can, however, draw the 
paths in £-plane, and then easily infer the nature of the 
paths in z-plane. 
" Now, from (1.8a) if U =0, then the path is given by 


5 1 0?+2ocosp-+tl 
T..|iog sing-+ 7 log 3 Spt cE 1| = const. 

Again,if =0(U=0), then the paths are 9= const 
i.e.the paths are radial. 

Before going into the actual graphical representation 
let us see how the line B.Din z-plane is transformed in 
the £-plane and also whatline in the z plane corresponds 
to the imaginary axis C„ O n of £-plane (Figs. 1, 2). 

As zchanges on BD, amp. zas wellas amp. dzremains 
fixed each being equal to n/4. Consequently amp. f(£) 
andamp.d{f(£)} must remain fixed on the corresponding 
path on Z-plane, and each must be equal to n/4. 

Hence, on the curve on £-plane corresponding to BD 


Z = amp. d{f(0} = amp. dt —amp.? 


— 2 amp. (€) amp. (+1), 
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which implies that B.Dis transformed on £-plane to the 
unitsemi-circle BO.D Fig. 2). 

Also with alittle analysis, one can easily see thät the 
imaginary axis 0. On corresponds to a line 0 „0A, in 
2-plane (where BO—=0O* 78 h nearly), the first quadrant of 
the Z-plane being transformed into the part of the chan- 
nel below the line CO. O0 A. and the second quadrant 
into the part above this line (Fig.1). This line is of 
great importance. For, it divides apparently the chan- 
nel into two parts such that the vortex originally in one 
part cannot, during its motion, cross into the other. 
This can be easilyseen from the nature of 0@. 

n 
2 
result that on the unitsemi-circle (0 = 1) and theimagi- 


Nowep=0gives o=1, orp= — leading to the 


nary axis (o = >) the vortex can have no crossradial 


velocity. t 

Again, putting 6 = 0 in the first equation of (1.7), 
one finds that the vortex will be stationary on the unit 
semi-circle at the value of pgiven by 


(2cosp9—1)—2sing=0!) . ... (21), 


En on the imaginary axis (p = en at the points given 
N 


oa © 02 1 een 
It may be noted that for real positive roots of (2.2) 


a must < — 2,andalsothatforra > — 2,0 <oO 
9= 


wa 


T 
ı)%& (- im) is to be taken positive if T'be*) and negative 


if Tbe/. It should be noted that the main channel flow indicates 
negative Totation oftheliquid in thefigures drawn. 


Fig. 3b 
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Fig. 4b 


always. And from the same equation it follows that if 
the vortex be stationary at a point e,(< 1)on the imagi- 
naryaxis,it will be also stationary at the inverse point 
1/o, outside the semi-circle. 

We can now draw the paths of the vortex. Writing 
equation (1.83)in the form 


0 +2ocoso-+]1l 
0? — 2ocosp-+] 


9=Fa 


(2.3) 
and determining for a particular value of &the values of 
F, at certain points on a number of radial lines, one 
finds byinterpolation the points on them corresponding 
to a fixed value of Fa; the join of these points gives a 
possible path of the vortex for that value of &. Any 
point on this path may be the initial position of the 
vortex. Similarly other possible paths for that value of 
a are drawn. Thus, giving to & different values the 
corresponding paths can be drawn. These paths on the 
&-plane are, however, drawn separately for each qua- 
drant. Figs. 3a, 4a, 5a show these paths, from which 
figs. 3b, 4b, 5b showing the nature of the actual paths 
in the channel can be sketched. 


3-Pathsin the Channel. 


The paths fall into three categories determined by 
the conditions: 


1 
%-|logsing + 7 log 


() a>o0, 
(ü)0>a >—2, 
MG) a<—2. 
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Category (i): : . 

In the first quadrant of the £-plane all the paths in a 
region adjoining the imaginary axis are closed, one of 
them reducing to the stationary point P on the semi- 
circle and the imaginary axis itself being the returning 
branch of a closed path, and in the remaining region 
close to the real axis the paths are all open running 
from C’„ towards inf. (Fig. 3a). (As a goes on decreasing 
remaining always positive the stationary point P will 
continuously approach the point B, and consequently 
the region close to the real axis in which the paths are 
open will become narrower and finally disappear with a). 
In the second quadrant all the paths are open running 
from inf. towards C,. Hence, we infer that in the part 
As 0 0% BA» (Fig. 3b lower) of the channel the vor- 
tex describes closed paths except in a narrow space close 
to the straight boundary where it continuously moves 
up the stream, whereas in the other part Au 0 Co DA 
(Fig. 3b upper) the vortex moves on continuously down 
the stream. 


Category (ü). 

In this category all the paths in the first quadrant of 
t-plane are open running from infinity towards the 
point C'. In the second quadrant the paths are of two 
types, namely (a) paths which cross the semi-circle and 
(b) paths which do not cross the semi-eircle. Under type 
(&) all the paths are open, some of which run from inf. 
towards CO, (these disappear only when = —2i.e. 
when the stationary point 8 coineides with O), while 
othersrun from C'„ towards inf. (these disappear when 
&= 0,8 then coinciding with D). The paths of type (b) 
exist only when «+0. Under this there are again two 
kinds of paths separated by the semi-circle, the kind 
outside the semi-cirele being just the image of that 
inside; both ends of each path of the kind outside the 
semi-cirele reach infinity, while the ends of those inside 
the semicircle converge to C„. Hence, we infer the 
following. Inthe part Au0 0. BA. (Fig. 4b lower) of 
the channel the vortex moves on continuously down the 
stream, whereas in the other part A. 00%. DA 
(Fig. 4b upper) the vortex can never cross the line BD 
if it be initially in any one of the two dotted regions 
(One in each leg) whose extension depends on & (#0) 
and which disappear for & = 0, and within each of 
these regions the vortex paths turn back. IT a=— 2, 
there appears between the central line and the above- 
mentioned regions a strip of area within which the vor- 
tex moves on continuously down the stream; and in the 
portion between the above-mentioned regions and the 
boundary Au DC. the vortex (&=# 0) crosses the line 
BDand moves continuously against the stream. 


Category (fü): 

In this category there are both open and closed paths 
in the first quadrant of the Z-plane, the open paths 
running from inf. towards 0’. (Fig. 5a); in the second 
quadrant the paths are, as in category (ii) of two types. 
In one type all the paths are open running from (0. 
towards inf., and the other type is somewhat similar to 
the type (b) Category (ii). It may be noted that in this 
case there are three possible stationary positions of the 
vortex, two T, T’ ]ying always on the imaginary axis 
(2.2) and the other Qlying on the semi-circle (Fig. 5a). 
All these three points however coincide with O when 
= — 2. Wenowinfer thatin the part Au O0 C„ BAs 
of the channel (Fig. 5b lower) the vortex moves on 
continuously down the stream exceptin a certain region 
depending on &% (which disappears when «= —2, 
points T, T’ and Q then coinciding with O) and within 
this region it describes closed paths round the stationary 
point Qon BO. In the other part Au O0 Cu DA, of the 
channel (Fig. 5b upper) the vortex will move against 
the stream crossing the line OD except when it lies in 
any one of the two dotted regions in which the paths 
turn back asin category (ti). 
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4. General case. 


If the channel be bent at any angle, 


problem can still be tackled. 
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nn say, the 


With the notations as before the transformation in 
this case er 


HK ae Da a ib a 


2. 
so that we may write 
2 fee ee (42) 
where 
ee 4 a Se 


Proceeding in the same way as before, the equation 
determining the path of the vortexis found to be given 


by 


aan a 
0? — 20cosp +1! 
= Uhop=const . .. . . (4.83). 


T'\logsin Br "10g 


The drawing of the paths is now very much more 
laborious, but the essential features of the types des- 
cribed before are expected to be found here also. 

I take the opportunity to express my indebtedness 
to Prof. N. R. Sen for his kind helpin the preparation 
of this paper. 
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Caleutta. K.R. Singh. 


Die Methode der konformen Abbildung ebener 
Spannungszustände von L. Föppl, ange- 
wandt auf Probleme des Parallelstreifens. 
Für die Untersuchung des durch konzentrierte Rand- 
kräfte belasteten Parallelstreifens stehen die von 
Filon!) und Seewald?) angegebenen strengen 
Lösungen zur Verfügung. Ihre numerische Auswer- 
tung stößt jedoch auf erhebliche Schwierigkeiten, 
wie z. B. die Arbeitvon G.Mesmer?) zeigt, der die 


 Filonsche Lösung für das Grundproblem des 


Parallelstreifens (zwei gegenüberliegende Lasten) zah- 
lenmäßig eingehend untersuchte. Mit zunehmendem 
Abstand vom Lastquerschnitt versagte bald infolge 
der immer schlechter werdenden Konvergenz die gra- 
phische Auswertung der Fourier- Integrale, deren 
Integrandenkurven sich schon in unmittelbarer Nähe 
des Lastquerschnittes über so weite Gebiete des Para- 
meters erstreckten, daßsie durch rechnerisch integrier- 
bare Näherungskurven fortgesetzt werden mußten. 


Der Zweck der vorliegenden Mitteilung ist, zu zeigen, 
daßsich diese Schwierigkeiten vermeiden lassen, wenn 
man zur Herleitung der Lösung die von L.Föppl®) 
angegebene Methode der konformen Abbildung ebener 
Spannungszustände benutzt. Dieseneue Lösung enthält 
zwar auch ein Fourier-Integral, das aber außer- 
ordentlich schnell konvergiert, da es im Gegensatz zur 
Filonschen Lösung die Last-Singularitäten nicht 
enthält und nur zur Aufhebung einer Schubspannungs- 
verteilung längs der Ränder benötigt wird, während 
Randnormalspannungen bei der Abbildung nicht auf- 
treten. 


Aus der Lösung des Grundproblems lassen sich leicht 
auch andere, technisch wichtige Spannungszustände 
ableiten, wie z.B. derjenige der Balkenbiegung, 
dessen Spannungsfunktion für den Fall des frei auf- 
liegenden, in der Mitte durch eine Einzelkraft belasteten 
Balkens mit überstehenden Enden angegeben wird. 


1) 1.N.S.Filon ‚Phil.Trans. A. Bd. 201 (1903), S. 106. 

2) F. Seewald, Abh. a.d. Aerodyn.Inst. d. T. H. Aachen, 
Berlin 1927. 

®2) G.Mesmer, Techn. Mechanik und Thermodynamik, Berlin 
1930, 8.85. 

“%) L.Föppl, Z. angew. Math. Mech. 11 (1931), S. 81. 


(1) 
cos yu+lojyo ’ 
2 Sinyv 

IT osyu +Eoiyv 


2. Konforme Abbildung der bekannten Spannungs- 
funktion der nach Bild 2 belasteten Walze (von der 


Stärke 1): 
2y an: 9) 


| 
In sin y u ‚ | 


+y—l1l 2 
auf den Parallelstreifen von der Stärke 1. 


z-Ebene 


N 26 arctg 


Bild 2 


Nach L. Föpplwird F,in die folgende Form um- 
geschrieben: 


2y 


Rn — lost nz 0) En 6,82 


a 4?) 
(3). 


Für die Spannungen ind len 
schrieben: 2 


2 Fu zu 


are A 
3. [en -( 30 ) 2° 

Da essich zeigt, daßsich mit den drei Gliedern des Aus- 
drucks (4) die Randbedingungen nichtidentisch erfüllen 
lassen, werden die beiden mit C, und c, behafteten 
Glieder, die von einem unwesentlichen hydrostatischen 
Spannungszustand in der z-Ebene herrühren und daher 
mit der Singularität des Problems nichts zu tun haben, 


weggelassen ((,=0;c,=0). Mit C, re aus der 


dritten Bedingung wird: 
12 Sinyv 
ne nen . (5). 
Für v= +b (sin yu= +1; cos yu= (0) wird 
er 
(ou)ı = 9»: > 0 


für jedes v mit Ausnahme von v = 
wird. 

Auf den Rändern verbleibt — außer den beiden Ein- 
zelkräften P — nur noch eine Schubspannungsver- 
teilung: 


_ 0®F, 
(Tuv)ı we, = AR d 


pP vKojyv — Sin 
-r2,( In u. ). 
yv 


die zum Verschwinden gebracht wird durch Überlage- 
rung der Spannungsfunktion: 


0, wo (au = =» 


Fi f (Ans it + Bio) = Sin Soon da 
0 


b 
(6). 
2F, 
5 sb. ei wird: A(&) = — B(a).0-Tga (7), 
OU 
= 0v r 
ergibt: 


Bi; ( yv&ojyv — Sin yo ) 
m’ Sin?yv 


f a\? [9 )) 
= I@) 30 ( Sina+ — x in u da 
0 


ei h 
Ou Ovu=ın 
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oder nach ® integriert, mit der Abkürzung, 


Die) = > 30) (ein + u) (8). 


x 


&v P yv 
D ee El DA 
| DE nz Eimyv 
0 


Mit Hilfe des Fourier-Theorems findet man hieraus: 


oo 
Pr vb 
D(e) = — ee ER 
n Eu b b 7 

2b 


{1} 


wenn mit J das Integral bezeichnet wird. Um den Wert 
von J zu erhalten, wird zunächst J,=f J da gebildet: 


° &v 
sin —— R 2sin —- 
Br ER) (5) - a% il 
er_. 3 
0 
2a 2& 
et —e 
= = 290°) 
er OLE f 
Damit wird 
dJ, 1 


und 


Schließlich folgt aus (8): 


2b 1 
een aCoja(Sin2a+2a) ' 


womit nach (7) auch A(&) bekannt wird. 


s) siehe Bierens de Haan: Nouvelles Tables, 8. 388. 
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Die Spannungsfunktion des Parallelstreifens zwi- 
schen 2 Einzellasten lautet somit: 
‚Einyv. 

cos yu 


F=-R+B= ZU E aro tg 


ke} 


j ein —% Ro 
2 ee N Be | 
Eoja (Sin 2% +20) 


x cos — “e) 


mit 


Bild 3 


Aus (10) folgen die Spannungen: 


an az | 2ycosyubßojyo 
“Tom m loosyu +Eintyo 
= en Sinne 
) > ( cos?yu + Sin?yv 


2 cos yu&oj?yv Sin yv ) 


E; (cos? yu + Ein? yo)? | 


("ein [a e 
fe & Sin 7 TgoaE&ei 7 


ar 
Ei a (Sin2a +20) 


a) 
AV® 
x c08S Z da & 


er cos yu Ein yv 
ou? m! cos?y u+ Sin?yv 


2 Sin yvsin?ywcosyu 
(cos? yu + Sin? yv)? 


a ES ur ein .- — a Tg a Coj 2 


„2P bu,imbs 
ab Ciao (Sin 2a +2) 


&V 
x cos do, 
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Sinyvsinyu 
EZ Mo 7’ \cos2yu+ Sinayv 


x sin? da. 


(cos? yu + Sin? yo)? 
2P | x 
Unb 
b 


Bild 3 zeigt für die Schubspannung 7,, an der Stelle 
u= 0,65; v = b die graphisch zu integrierende Kurve 
sowie zum Vergleich diejenige nach der Lösung von 
Filon (gestrichelt). 

Auch in größeren Abständen vom Lastquerschnitt 
klingen die Integrandenkurven der Lösung (10) schnell 
ab. So wird z. B. der Wert des zu (o4),_, gehörigen 

v=2b 
= 4 schon verschwindend klein. 


‚ AU AU au 3 au 
(ein +60 —arıa Sin “) 
Coja (sin2a +20) 


Integranden für & 


Bild 4 


Aus dem Grundproblem läßt sich auch leicht die 
Spannungsfunktion für den durch Einzellasten ge- 


Buchbesprechungen 


Z. angew. Math. Mech. 
Bd. 34 Nr. 10/11 ‘Okt. me 1954 


bogenen Balken herleiten, die z. B. für den in Bild 4 


dargestellten Fall folgendermaßen lautet: 


he A 2 ie. 
ma Gyr + cos yu — sin yu 
Sinyw+Hi) 

= = +I) are tg Eofyw+D-+Fcosyutsin yu yu 

I „ne zu, 

. w—) are tg Cojy(o —N+teosyu+tsinyu 


= 


Av 


= / [46 Ci + D(a) — = ein &#|(cos = y 
4 ze: 

x | +5 008 [»° zu da 

5 &i- &# (co = 


a B El Goos =)“. 


I cos E 


+| [36 Sin = +C(e) — 
f) 


2 b 2 b 
a. RL 

7 ja (Sin2a +2a)’ 

Pb 1 


 % Sofa (Sin2a+2o)’ 
1- 3 
4 (® Sin2a+a Sin 2a Sin?a 
x2Cin?20 (Sin2a — 20) 
2 (Coja — 1) (Sina — Eoja) 
2 Sin2a(Cin2a—- 20) |) 


ee Coja (Sina +a)—a 
n 03GSin2a(Gin2a — 20)’ 
Karlsruhe. Rudolf Sonntag. 


BUCHBESPRECHUNGEN 


Dr. A. Hess (ehem. Prof. am Kantonalen Technikum 
in Winterthur) Praktische Mathematik. Für 
Studierende der Technik und zum Selbststudium. 116 S. 
m. 127 Abb. Zürich 1947. Rascher Verlag. Preis 
geh. 8.90 DM; 

Das vorliegende Bändchen enthält als Ergänzung zu 
den Lehrbüchern der Planimetrie, Trigonometrie und 
Analytischen Geometrie des gleichen Verfassers rd. 
150 Aufgaben mit Lösung aus diesen Gebieten. Außer- 
dem gibt es eine kurze Einführung in einige Kapitel der 


praktischen Analysis (Funktionspapiere, genäherte 
Integration und Nomographie). 
Dresden. Willers. 


Georg Aumann (o. Prof. a. d. Univ. München), 
Reelle Funktionen. (Grundlehren der 
mathematischen Wissenschaften. Bd. LXVIII.) 
VIII +416S. m. 22 Abb. Berlin/Göttingen/Heidel- 
berg 1954. Springer-Verlag. Preis geh. 56,— DM, 
geb. 59,60 DM. 

Das lang erwartete Buch des bekannten Verfassers 
liegt nun dem mathematisch interessierten Publikum 
vor. Der über 400 Seiten starke Band gibt dem Leser 
einen vorzüglichen Überblick über den augenblick- 
lichen Stand der Theorie der reellen Funktionen. Auf 
gewisse Lücken, es wird z. B. nicht auf die Theorie 
der Distributionen eingegangen, weist der Verfasser 
selbst hin. An Kenntnissen setzt das Buch verhältnis- 
mäßig wenig voraus, seine Lektüre erfordert aller- 
dings, schon wegen der konsequenten Verwendung 


logistischer Zeichen, eine gewisse „Bereitschaft zur 
Abstraktion‘. Die ersten vier Abschnitte behandeln 
abstrakte Mengen, Ordnungen (teilweise, linear ge- 
ordnete und wohlgeordnete Mengen, Zornsches Lem- 
ma), Verbände und Räume. Die nächsten zwei Ab- 
schnitte gehen auf Punktfunktionen ein — man fin- 
det hier auch den Weierstraß-Stoneschen Approxi- 
mationssatz und den Brouwerschen Fixpunktsatz — 
und auf Funktionen in Produkträumen. Der 7. Ab- 
schnitt beschäftigt sich eingehend mit den reellen 
Funktionen einer reellen Veränderlichen. Im Zu- 
sammenhang mit dem sehr allgemein gehaltenen Auf- 
bau einer Differential- und Integralrechnung wird 
hier u.a. das Perron-Integral und das Lebesguesche 
Integral behandelt. Der 8. Abschnitt befaßt sich aus- 
führlich mit der Maßtheorie, während der letzte Ab- 
schnitt sehr elegant (nach M.H. Stone) eine all- 
gemeine Integrationstheorie wiedergibt, die verhältnis- 
mäßig einfach (und ohne Maßtheorie) zu den grund- 
legenden Sätzen gelangt. Der Abschnitt schließt mit 
der Behandlung iterierter Integrale und einer neu- 
artigen Fassung des Fubinischen Satzes. Am Ende 
des Buches befindet sich ein sorgfältig zusammen- 
gestelltes Literaturverzeichnis, das besonders auf ein 
weiteres Eindringen in die Theorie zugeschnitten ist, 
und ein ausführliches Namen- und Sachregister. Ein 
näheres Eingehen auf das inhaltsreiche und vom 
Springer-Verlag in gewohnter Weisse hervorragend 
ausgestattete Buch ist hier nicht möglich. 


Dresden. M. Landsberg. 
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 A.H. Shapiro, (Prof. am Massachusetts Institute 

of Technology). The Dynamics and Ther- 
modynamics of Compressible Fluid 
Flow BdI. XIII -+647 Seiten m. zahlr. Abb. 
New York 1953. The Ronald Press Company. Preis 
geb. 16 $. 

Der Verfasser hat es verstanden, theoretische Über- 
legungen, physikalische Beschreibungen und experi- 
mentelle Erkenntnisse auf allen Gebieten der Gas- 
dynamik zu verbinden. 


Neben der für Rohrströmungen wichtigen ein-. 


dimensionalen Behandlung wird auch ausführlich die 
zweidimensionale Strömung im Unter- und Über- 
schallbereich dargestellt. Die modernen experimen- 
tellen und theoretischen Methoden zur Analyse der 
Strömungsfelder nahe der Machzahl 1 sind nachzu- 
lesen. Der Ingenieur findet Durchrechnungen der 
Gitterströmung. 

Das Buch ist als Nachschlagewerk, aber auch für 
das Selbststudium zu empfehlen. Der Lernende findet 
am Ende eines jeden Kapitels eine Anzahl von Auf- 
gabenstellungen. 


Dresden. Albring. 


Dr.-Ing. Bruno Eck, Technische Strö- 
mungslehre. Vierte verb. Aufl. X + 4228. m. 
407 Abb. Berlin/Göttingen/Heidelberg 1954. Sprin- 
ger-Verlag. Preis geb. 29,40 DM. 

Das Buch von Eck, das 1940 erstmalig erschien 
und nun in der 4. Auflage vorliegt, wird von Tech- 
nikern gern gelesen; denn der Verfasser hat den Ver- 
such gemacht, die technisch wichtigen wissenschaft- 
lichen Erkenntnisse der modernen Strömungslehre 
aus den Spezialwerken und den Veröffentlichungen 
der Fachzeitungen herauszulösen und sie in einer dem 
Praktiker verständlichen Form darzustellen. Dabei 
kommt dem Verfasser besonders sein guter Überblick 
über die modernen Forschungsarbeiten zu statten. 
Allerdings bleiben noch viele Wünsche offen. Das 
bezieht sich keineswegs auf den Umfang der Dar- 
stellung — von der Staubabscheidung, von den Ver- 
brennungsvorgängen bis zur Gasdynamik ist kaum 
ein Anwendungsgebiet ausgelassen — sondern auf die 
Art der Darstellung. Dazu einige Beispiele: Das 
Gebiet der Grenzschichten ist für die moderne Technik 
von allergrößter Bedeutung. Die Abhandlung von 
Eck hat in methodischer Hinsicht manche Mängel, 
auch sachliche Korrekturen sind hier und dort not- 
wendig. So trifft der Hinweis, daß das Diagramm 


A—=f —— für Kreisrohre bei Umrechnung auf den 
v 


hydraulischen Radius nur noch für die- turbulente 
Strömungsform benutzt werden kann, nicht zu, wie 
schon die Messungen von Nikuradse und 
Schiller (Handbuch der Experimentalphysik 
Band IV, Bild 96 und 97) zeigten. Diese Auftragung 
hat natürlich auch für den laminaren Bereich Gültig- 
keit. Das Kapitel über den Umschlag laminar — tur- 
bulent liefert keine Auskunft über die Berechnungs- 
methoden, die z.B. von Schlichting in einer 
auch für den Praktiker auswertbaren Form (kritische 
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Reynoldszahl über dem Pohlhausenparameter) zu- 
sammengestellt sind. Es bleibt auch unklar, wann 
die Umschlagkriterien von Tollmien und wann 
die einen ganz anderen physikalischen Vorgang zu- 
grunde liegenden Kriterien von Taylor zur An- 
wendung kommen. Im Kapitel Gasdynamik vermisse 
ich das Arbeiten mit der Stoßpolaren nach Buse- 
mann. Dieses Diagramm ist für den Ingenieur 
außerordentlich nützlich und ermöglicht die Ab- 
lesung von Stoßwinkel, Stoßverlust, Ablösung der 
Strömung abhängig vom Keilwinkel und der Mach- 
zahl. Über den A-Stoß (auch als gegabelter Ver- 
diehtungsstoß bekannt) wird erstmalig in dieser Auf- 
lage referiert. Von seiner Wirkung — der Strömungs- 
ablösung — wird aber nicht gesprochen. Jedoch 
sollen diese Hinweise nicht den grundsätzlichen Wert 
eines Buches bestreiten, das dem Praktiker Hilfe und 
Stütze sein soll, aber ein Werk von diesem Umfang 
muß eine dauernde Weiterentwicklung und Ver- 
besserung erfahren. 
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Gestützt auf drei eigene frühere Arbeiten zur Be- 
rechnung turbulenter Grenzschichten, entwickelt der 
Verfasser hier zwei Berechnungsverfahren für in- 
kompressible Strömungen bei glatten und bei rauhen 
Wänden. Das erste (E-Verfahren) benützt nach dem 
Vorgang von Gruschwitz neben dem Impuls- 
satz der Grenzschichttheorie noch eine zweite, aus 
dem Energiesatz hergeleitete Differentialgleichung 
für die Änderung der Profilform. Das zweite (B-Ver- 
fahren) setzt wie bei Buri eine unmittelbare Ab- 
hängigkeit zwischen Profilparameter und Druck- 


“ gradient voraus. Hier ist die Integration nur einer 


gewöhnlichen Differentialgleichung erforderlich. Als 
grundsätzliche Annahmen dienen: Aufbau der Grenz- 
schicht aus wandnahem und äußerem Teil, Einpara- 
metrigkeit der äußeren Profilteile, eindeutige Bezie- 
hung zwischen Schwankungsverteilungen und Ge- 
schwindigkeitsprofilen. Mit beiden Verfahren werden 
berechnet Impulsverlustdicke, Verdrängungsdicke und 
örtlicher Reibungsbeiwert. Durchgerechnete Beispiele 
ergeben für die glatte Platte bei beiden Verfahren be- 
friedigende Übereinstimmung mit Meßwerten. Für 
Grenzschicht mit starkem Druckanstieg ergibt sich 
bei glatter und rauher Oberfläche mit dem E-Ver- 
fahren gute Übereinstimmung mit Meßwerten von 
Schubauer, ebenso für das FlügelprofilNACA 65 
(226) — 222 (Approx). In Gebieten fallenden Drucks 
oder mäßigen Druckanstiegs sind die Differenzen 
zwischen beiden Verfahren nur unbedeutend. 


Weilburg. L.Schiller. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 
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NACHRICHTEN 


Freiburg i.Br.: Die von der Gesellschaft für 
Angewandte Mathematik und Mechanik (GAMM) 
ins Leben gerufene Stellenvermittlung 
zur Förderung des wissenschaftlichen Nachwuchses 
auf den Gebieten der Mathematik und Mechanik ist 
nach einem auf der Münchener Jahrestagung 1954 
dieser Gesellschaft vollzogenen Beschluß an den 
Interessenverband deutscher Mathematiker (IVM) 
zur Weiterführung übergeben worden. Bewerbungen 
sind daher in Zukunft an die Anschrift dieses Ver- 
bandes, Stuttgart-Degerloch, Felix-Dahn-Straße 77, 


zu richten, Die Geschäftsstelle der GAMM stellt 
sich jedoch für Fragen der Studien- und Berufs- 
beratung der Diplom-Mathematiker nach wie vor 
zur Verfügung. 

Die auf der Münchener Jahrestagung 1954 der ° 
GAMM gehaltenen Referate über Ausbildung und 
Stellenvermittlung der Diplom-Mathematiker mit 
Auszügen aus der anschließenden Diskussion er- 
scheinen Ende Juli 1954 als kleine Broschüre im 
Physik Verlag, Mosbach. 

H. Görtler. 
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